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A discussion is given of the different possibilities which 
would lead to the assumption of the presence of mesotrons 
in extensive atmospheric showers, or the so-called Auger 
showers. The difficulties involved in the problem of the 
experimental evidence for the existence of these mesotrons 
are shown to be, in general, considerable. Three kinds of 
experiments are described in which mesotrons associated 
simultaneously with other particles are found. These 
particles are present even at large distances (~20 meters) 
from the observed mesotrons. However, they are much 
more infrequent than the particles which are found in 
the Auger showers, investigated under conditions similar 
to those of the present experiments. Qualitatively, the 
observed phenomena could be due, not to mesotrons 


associated with an Auger shower, but to single mesotrons 
accompanied by collision electrons which are projected by 
these mesotrons along their path in the air. There exists, 
however, a definite quantitative disagreement between the 
observed values and the theoretical values which are 
calculated for the considered collision processes. Therefore, 
at least part of the observed mesotrons could actually 
belong to an Auger shower. But, on the other hand, 
explanations are formulated which might, at least partly, 
account for the mentioned quantitative disagreement. 
Thus, only after these explanations are tested by experi- 
ment will it be possible to decide more definitely as to the 
presence of mesotrons in Auger showers. 





PRELIMINARY CONSIDERATIONS 


HE problem of the presence of mesotrons in 

extensive atmospheric cosmic-ray showers, 
or Auger showers (here referred to as A showers) 
has been investigated since the first experiments 
were carried out on these showers.' The par- 
ticular interest of this problem arises from its 
possible relation to different phenomena, which 
are mainly: (a) the development of a cascade 
shower originating from a particle of very high 
energy, up to 10'* ev, and (b) the mechanism 
itself of the generation of mesotrons in cosmic 
rays. 





'P. Auger, R. Maze, and T. Grivet-Meyer, Comptes 
rendus 206, 1721 (1938); W. Kolhérster, J]. Matthes, and 
E. Weber, Naturwiss. 26, 576 (1938); P. Auger, P. Ehren- 
fest, R. Maze, J. Daudin, Robley, and A. Freon, Rev. 
Mod. Phys. 11, 288 (1939); H. Euler and H. Wergeland, 
— 27, 484 (1939); N. Hilberry, Phys. Rev. 60, 1 

1). 
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We shall briefly examine here the different 
possibilities which would lead to the assumption 
of the presence of mesotrons in A showers, and 
then discuss the problem of their experimental 
evidence. 

An analogous discussion can be applied, in 
part, for heavy penetrating particles, other than 
mesotrons, which might also be present in A 
showers. 


A. Showers Originating from Primary 
Soft Particles 


In case the primary soft particles (electron or 
photon) multiply only by processes described in 
the theory of cascade showers and if one admits 
that the cascade theory is valid for particles 
having energies up to 10'* ev, then one should 
not expect to find mesotrons in A showers since 
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this theory does not predict at all any creation 
of mesotrons in the multiplication process leading 
to the building up of the shower. 


B. Showers Originating Indirectly from 
High Energy Mesotrons 


An A shower can originate from a mesotron 
by a close collision process with an electron,” * by 
a bremsstrahlung process providing a high energy 
photon,’ or possibly by a decay of a short life 
mesotron of very high energy. The main differ- 
ence between cases A and B is that the high 
energy parent particle of the A shower is a 
secondary instead of being a primary. In the 
present case, however, the A shower starts at 
any altitude in the atmosphere wherever the 
high energy secondary particle is created, and 
not at the top of the atmosphere, as in the case 
of the primary parent particle. We will leave 
aside the question of how far the origin of the A 
showers could be attributed to such mechanisms 
since our intention is to limit the present dis- 
cussion only to the problem of the presence of 
mesotrons in A showers. 

It is expedient to remark here that in order to 
obtain a widespread shower, the energy transfer 
from the mesotron to the high energy soft 
particle should have taken place in the upper 
atmosphere, otherwise only a burst of particles 
of a small extension could be observed. In any 
case, the probability of detecting in the whole 
shower the only present mesotron—which is the 
source of that shower and is assumed to have 
lost only part of its energy in the above-de- 
scribed processes—will be very small. On the 
other hand, the probability of detection of this 
mesotron should decrease at least in proportion 
to the frequency of its occurrence in the differ- 
ential energy spectrum of the mesotron radiation. 


C. Secondary Production of Mesotrons 


1. The problem of detecting mesotrons in the 
A showers would look less. hopeless if there 
should be a creation of a larger number of 
mesotrons within these showers. In case the 
created mesotrons are of high energy, they must 
have originated in the core of the A shower. 
~ #H. J. Bhabha, Proc. Roy. Soc. 164, 257 (1938). 


*R. F. Christy and S. Kusaka, Phys. Rev. 405, 411 
(1941). 


They will travel together with the high energy 
particles of the central part (core) of that shower 
and remain close to its axis. 

In order to detect these mesotrons one will 
have to eliminate from the core of the A shower 
the soft high energy particles which are mixed 
together with the mesotrons. 

2. From the experimental point of view, the 
identification of the high energy mesotrons in 
the core of the A shower will present great diffi- 
culties. Indeed, in counter experiments, the only 
reliable criterion for the identification of the 
mesotron is its small probability of creating a 
shower after having traversed a larger number of 
radiation units of any material. Generally, it is 
sufficient to place below this material a group 
of counters, the function of which is to distinguish 
between a single particle or a shower emerging 
below the material. If, however, there are present 
many mesotrons of high energy, they could dis- 
charge this group of counters and would act then 
in the same way as a cascade shower. The above 
criterion will fail in this case. It will therefore 
be necessary to increase still more the thickness 
of material in order to eliminate eventually the 
less energetic mesotrons until one single mesotron 
remains. Such a procedure could lead to the use 
of prohibitively large thicknesses of absorbing 
material. However, if one could rely on the 
strict validity of the cascade theory for ex- 
tremely high energies (10'*-10'*® ev) of the soft 
particles which are involved in the A showers, 
then all particles, whatever their number might 
be, emerging below, for example, a lead absorber 
of a thickness of, say, } meter, should consist of 
penetrating particles (mesotrons), or at any rate 
of particles originating from a penetrating par- 
ticle. This follows from the fact that any soft 
particle of an A shower of energy up to 10" ev 
should be stopped in the indicated lead absorber 
since the cascade theory shows that for high 
energy soft particles the average cross section for 
the pair production and the fractional energy 
loss by radiation processes are practically inde- 
pendent of the energy of these particles and are, 
for one radiation unit of length of any material 
(0.5 cm for Pb), equal to ~1. In any case, there 
will always be a possibility that the mesotrons 
finally identified were, originally, not present at 
all in the A shower, but were created in the 
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absorber—for example, a lead absorber—where 
the cross section for such a process might be very 
different from that in air. Therefore, it would be 
best to use, as absorbers, materials having about 
the same atomic number as air, e.g., water. 
Such experiments could be carried out by im- 
mersing the detecting apparatus in a lake. 

3. In cases where one could demonstrate that 
the particles emerging below the thick absorber 
did not belong primarily to a dense part of an A 
shower, in other words, if one could prove that 
the density of particles which are then present 
above the absorber is relatively low, the particles 
below the absorber could be attributed to a 
mesotron with a considerably higher probability 
than in the previous case. Indeed, in the present 
case, the particles below the absorber could not 
have originated from a high energy electron or 
photon, since such particles would almost always 
be accompanied either by other particles of the 
core of the A shower, or by a cascade shower 
originated along their path in air, or even more 
frequently by both. It is evident that these cases 
in which the mesotrons are present in the low 
density region (wing) of the A shower are the 
most favorable ones for their detection. However, 
as will be seen more in detail in the experimental 
part of this paper, the question may then arise 
whether the observed mesotrons belong actually 
toan A shower or whether they can be considered 
as single mesotrons ejecting electrons along their 
path in air by collision processes. 

4. Finally, it should be pointed out that in 
case the mesotrons created within the core of 
the A showers are of low energy, of the order of 
magnitude of their rest mass energy, the proba- 
bility of their detection would again be small 
since they would very likely be stopped in an 
absorber of a thickness smaller than required to 
filter out all the high energy soft particles present 
in the core. They would naturally be stopped 
already in the air, in case they are created high 
in the atmosphere. 


D. Delayed Particles in A Showers 


As a direct consequence of the possible presence 
of low energy mesotrons in an A shower there 
should exist a certain number of electrons, re- 
sulting from the decay of the mesotrons, and 
occurring with delays—of the order of magnitude 
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of the mean lifetime (~ 2X 10~* sec.) of the meso- 
tron—with respect to all other particles in the 
same parent A shower. Inversely, a demonstra- 
tion of the existence of such delayed particles in 
an A shower would provide probably the best 
evidence of the presence of mesotrons in A 
showers. For that purpose, an apparatus was 
built which makes it possible to measure delays 
occurring between a group of counters selecting 
A showers and another counter placed in the 
neighborhood of this group of counters. This 
outfit records only those delays which occur 
within ~.10X10-* sec. A more detailed account 
of this method will be given at a later date. 


EXPERIMENTS 


Several’ types of counter-tube experiments 
were performed in order to obtain direct evidence 
of the possible presence of mesotrons in A 
showers. These investigations were carried out at 
sea level (Chicago, altitude 610 feet) and at a 
higher altitude* (Echo Lake, Colorado, altitude 
10,618 feet). In all experiments here described 
the master group method was used. The general 
principle of this method was described in a 
previous paper,® to which we shall refer here as I. 
All counters used in the present experiments are 
filled with a mixture of argon and alcohol.® 
They have an effective area of 5X50=250 cm? 
and operate at 1200 volts. The resolving time of 
the coincidence circuits is about 210~* sec. for 
a separation between counters less than 2 meters 
and about 4X10~-® sec. for a separation of 22 
meters. 

Throughout this paper, any counter which 
does not belong to the master group will be 
called an analyzing or distant counter. When 
speaking about discharges or counting rates in 
such a counter, it should be understood, unless 
specified otherwise, that the master group was 
simultaneously discharged. 

The discharges of the analyzing counters were 
recorded on a photographic film by means of 
neon lamp flashes. The photographic recorder 
was of a type similar to that used by Schein, 


‘ During the Mount Evans Expedition of the University 
of Chicago in August-September, 1943. 

5 A. Rogozinski, Phys. Rev. 65, 207 (1944). 

6 These counter tubes were designed by Dr. V. H. 
Regener. 
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Jesse, and Wollan’ in their stratosphere ex- 
periments. 


Experiment I 


In the first experiment,’ which was carried out 
at sea level and was of a preliminary character, 


© © 
© 


© 
© 


—— vp to 1S meters 


Fic. 1. Counter arrangement in Experiment I. 


the master group consisted of two counters M, 
and M, connected in twofold coincidence and 
separated by a horizontal distance up to 15 
meters. The role of this master group was to 
select A showers (Fig. 1). A set of 4 analyzing 
counters F, G, G’, H, placed below M, could be 
separated from that counter by lead absorbers 
up to 12.5 cm of thickness. The counters M2, F, 
G, G’, and H were located inside a Celotex box 
which was placed in open air and could be heated. 
The counter M, was placed indoors below a thin 
tile roof. 

It was found, for different thicknesses of the 
lead absorber, that a certain fraction of the four- 
fold coincidences M,M2GG’ was not accompanied 
by a discharge in any of the counters F and H. 
Therefore, it is assumed that these events could 
be attributed to single non-shower producing 
particles present in that part of the A shower 
which strikes the lead absorber. These particles 
which emerge below a lead thickness of several 


radiation units without producing an abundant . 


shower should be heavy particles of the pene- 
trating type, such as mesotrons. 


Experiment II 


1. In order to determine the density of par- 
ticles in the region of the shower in which the 
penetrating particle is located, experiments were 
performed (at Echo Lake) with the same counter 
arrangement as described in I. In this experiment 
the master group consisted of two counters M, 
and M:2 connected in twofold coincidence and 
placed one above the other. In order to increase 

7M. Schein, W. P. Jesse, and E. O. Wollan, Phys. Rev. 


59, 615 (1941). 
’ P, Auger and A. Rogozinski, Phys. Rev. 65, 64 (1944). 


the counting rate, two counters connected in 
parallel were used for M,; and for Mz (Fig. 2), 
At a horizontal distance of 22 meters from the 
master group there was placed the set of seven 
analyzing counters K,, Ke, K3, F, G, G’, H, of 
which the upper set of counters, K;, Ke, K,, 
could be separated from the lower set, F, G, 
G’, H, by variable thicknesses of a lead absorber. 
The four counters of the master group were put 
outdoors in a small wooden box: which had a 
wall thickness of § in. and which could be heated. 
The set of analyzing counters was located in- 
doors under a thin wooden roof. 

2. The results obtained in Experiment |] 
dealing with those cases where a shower emerges 
below a lead absorber (coincidences M,M2FGG'H) 
were already described in I. Here we shall in- 
vestigate in detail such events, which could be 
attributed to penetrating particles, like meso- 
trons, emerging below the lead absorber without 
being accompanied by a shower. These events are 
recorded as coincidences M,M2GG’, which are 
not accompanied simultaneously by a discharge 
in any of the counters F and H. Such an “‘anti- 
coincidence”’ will be represented by the symbol 
M,M.GG'— FH. . | 

The results of these experiments are sum- 
marized in Table I. However, for comparison 
with those’ cases in which a shower emerges 
below the lead absorber, the corresponding data 
are also given in this table. The first column of 
Table I gives the total time of recording for a 
certain thickness of lead, which is indicated in 
column II. The columns from III to XI repre- 
sent average counting rates, expressed in num- 
ber of coincidences per 10 hours. Column II] 
gives the frequencies of the anticoincidences 


OO 
8 
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Fic. 2. Counter arrangement in Experiment II. 


M,M.GG’' — FH or the coincidences M,M,FGG'H, 
regardless of discharges which might have oc- 
curred simultaneously in the upper set of counters 
K,, Ke, Ks3. Columns IV to X give the fre- 
quencies of those events M,M.,GG’—FH or 
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TABLE I. Experiment II. Coincidences, per 10 hours, occurring between the master group M,Mz¢ 
and the different analyzing counters F, G, G’, H, Ki, Ko, Ks. 








I IV Vv 


III 
MiMGG’ —FH (1) 
Hours of or 
MiM2FGG'H (2) 


recording 


Ki Ke 


VI VII Vill IX x 
MiMGG’ —FH or MiM:FGG’H 
in coincidence with 


Ks KiK2 KikKs; KK; KiK2K; 





(1) 75 
(2) 66 3 
(1) 56 
(2) 64.8 
(1) 27.8 
(2) 26 


8.4 
3.7 
3.4 
0.5 


7.5 6 3 0 3 

4.5 9 3 12 22.5 
1.7 3 
3.7 6 33 
1.5 1 
20 


6.7 1.3 8.4 


3.4 6.4 
4.4 1.5 1 


0.5 1 1 0.5 








M,M:2FGG'H which are accompanied by simul- 
taneous discharges in the 7 possible combinations 
of counters of the upper set K,, Ke, Ks. The last 
column A represents the differences between the 
frequency of M,M,GG’—FH or M,M2FGG'H 
and the sum of all the frequencies of the cor- 
responding 7 types of coincidences listed in 
columns IV to X. These A values represent 
the frequency of events M,M.GG’—FH or 
M,M:FGG'H which are not accompanied by a 
simultaneous discharge in any of the upper 
counters K;, Ke, K3. The average frequencies 
(counts per 10 hours) of coincidences between the 
master group and just one, or just two, or all 
three counters K,, Ke, K3, were, respectively, 
320, 80, and 170, regardless of any discharge 
which might have occurred in the lower set of 
counters. It results from these figures that the 
average frequency of those A showers, which 
discharge at least one of the counters K,, Ko, 
or K; in coincidence with the master group, is 
320+2 X 80+ 170=650 in 10 hours. 

3. As we have mentioned in I, the chance co- 
incidences could be neglected for the higher 
order coincidences reported there. In the present 
case, however, this factor has to be considered, 
mainly for the single-particle events M,M.GG’ 
-—FH, where larger thicknesses of lead were 
used. Although such an event is essentially a 
fourfold coincidence, the computation of the. 
chance coincidences in the actual case involves 
coincidence rates of different orders, the most 
important contribution of which comes from the 
twofold coincidence rates between (a) the proper 
twofold coincidence rate M,M, of the master 
group and (b) the proper twofold coincidence 


rate of the telescope GG’. The frequencies of the 
chance coincidences (per 10 hours), thus calcu- 
lated, for M,M,GG' — FH are, respectively, 19.2, 
18, and 16 for lead thicknesses of 2.5, 5.0, and 
10.0 cm. This shows that the corresponding A 
values for 5.0 and for 10.0 cm of lead could be 
due entirely to chance coincidences. For 10.0 cm 
of lead, the A value of 10.6, which is appreciably 
lower than the value of 16 which represents the 
calculated chance coincidence frequency, could 
be due to a fluctuation in the average number of 
the chance coincidences; in addition, it could be 
caused by a small contribution of chance coin- 
cidences arising from the fivefold coincidences 
M,M.GG'K,, M,M.GG' Kz, and M,M.GG'K;3. 
Incidentally, the above-mentioned reasons— 
apart from the fact that the geometrical arrange- 
ment of the counter telescope GG’ presents a 
certain “‘leak,’’ regarding the solid angle of GG’ 
which is not entirely covered by the upper set 
of counters—makes uncertain, on the basis of 
the present data, any attempt to correlate with 
a non-ionizing particle the passage of a single 
penetrating particle through the lead absorber. 
4. In order to obtain a clearer picture of the 
meaning of the experimental data of Table | 
concerning the density of particles which strike 
the lead absorber in case of events M,M,GG’ 
— FH, a comparison was made in Table II be- 
tween these events and the shower events which 
are respectively coincident with a discharge 
occurring in one or more of the upper counters 
K,, Ke, K3. For this purpose, the counting rates 
for M,M,GG' — FH and M,M,FGG'H were each 
reduced to a value of 100, and the counting rates 
of the corresponding higher order coincidences 
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were calculated accordingly. The equivalent 
three counting rates for coincidences between 
M,M.GG'—FH or M,M:2FGG'H on one hand, 
and K,, Ke, or Kz; on the other hand, were re- 
placed by their sum. The same was done for the 
three counting rates involving K,K2, K,K3, or 
K2Kz3, instead of Ki, Ke, or K3. In computing 
the values of Table II, the calculated chance 
coincidence rates of 19.2 and 18 were taken into 
account for 2.5 and 5.0 cm of Pb, respectively ; 
the chance coincidence rate for 10.0 cm of Pb 
was taken equal to 10.6 in order to make the A 
value which corresponds to this thickness of lead 
equal to zero. 


Discussion of Experiment II 


Let us consider, in Table II, the data referring 
to single or shower events which occur below 
the lead and are associated with simultaneous 
discharges in the upper set of counters K;, Ko, Ks. 
One can notice immediately the apparent fact 
that for all lead thicknesses, the single-particle 
events M,M.2GG’ — FH are much more frequently 
accompanied by discharges in a single upper 
counter than the shower events M,M.,FGG'H. 
These, on the contrary, are more frequently 
accompanied by discharges in two and mainly 
in three upper counters than the single-particle 
events. In particular, for 10 cm of Pb, the ratio 
between the percentages of events M,M,.GG’ 
—FH or M,M:FGG'H, each accompanied by 


TABLE II. Experiment II. Relative values of counting 
rates for coincidences occurring between (1) M,M2GG' 
— FH, or (2) M:M:FGG'H, and the upper counters 
Ki, Ks, K3. 








M.M:2GG’ —FH or MiM2FGG'H 
in coincidence with 


M.1M:GG’ —FH (1) 





or (Kit+Ke (KiK2+Kiks 
cm Pb M:1M2FGG’H (2) +Ks3) +K2K:3) KiK2K;3 A 
(1) 100 54 16 5.5 24.5 
2.5 
(2) 100 16 36.5 34 13.5 
(1) 100 60 16 22 2 
5.0 
(2) 100 13.5 25 51 10.5 
(1) 100 71 23 6 0 
10.0 
(2) 100 7.7 9.6 77 5.7 








discharges in the upper set of counters, changes 
from 71:7.7-~~9 when a discharge occurs in just 
one of the upper counters, to 6:77~1/13 when a 
discharge occurs in all three upper counters 
K;, Ke, K3. This variation is of the order of mag- 





ROGOZINSKI 


nitude of 9:1/13~100. The obvious consequence 
which follows from these figures is that the mean 
density of particles which are associated, above the 
lead, with single-particle events M@,M.GG’ — FH 
is very small in comparison with the correspond- 
ing density for the shower events /@,M2FGG'H. 

This fact that a single-particle event M,M.GG’ 
— FH, which is presumably caused by a meso- 
tron, is usually associated with a very small 
number of ionizing particles in its immediate 
neighborhood, could be interpreted in two differ- 
ent ways: Either the mesotron belongs to a low 
density region (wing) of an A shower and, hence, 
is present at a relatively large distance from the 
core of the A shower (A shower hypothesis), or 
the mesotron is located in a shower of low par- 
ticle density, other than an A shower. Such a 
shower could arise, for example, from collisions 
of the mesotron with electrons along its path in 
the air (collision-electron hypothesis). 

Thus, if the A shower hypothesis is correct, 
the mean density of particles which cover a 
certain relatively large area around the mesotron 
should be approximately constant all over this 
area. In other words, one should expect that the 
considered particle density at any point of this 
area would be approximately independent of the 
distance (a few meters) of this point from the 
mesotron trajectory since in the present case the 
mesotron would be located far from the high 
density region of its parent A shower. 

In case the collision-electron hypothesis is 
correct, the just considered density of particles 
should be an essentially decreasing function of 
the distance from the mesotron trajectory. This 
follows from the fact that, on one hand, electrons 
arising from collision processes of mesotrons in 
air should have an increasing energy in order to 
reach points located at increasing distances from 
the parent. mesotron, and, on the other hand, 
the probability of ejecting a collision electron of 
a certain energy decreases, in first approximation, 
as the square of that energy. A theoretical com- 
putation of this effect is given in the paper (II]) 


which follows the present one. 


Experiment III 


As a consequence of these considerations 
further attempts were made to determine which 
one of the considered hypotheses is correct. 
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1. In the experiments I and II the master 
group could be discharged by both soft and 
penetrating particles of energies higher than the 
minimum energy of a few Mev which is necessary 
to penetrate the counter walls. In the experi- 

















(a) 














(b) 


Fic. 3. Counter arrangement in Experiment III. 


ments which are to be described, the master 
group selected only particles which are able to 
penetrate 10 cm of lead, and, in particular, those 
particles which emerge below the lead without 
being accompanied by a shower and are, there- 
fore, classified as mesotrons. One can say that 
in the latter case the master group worked as a 
kind of a mesotron selector. 

Thus, the master group consisted of a vertical 
telescope in twofold coincidence M,Mz and, in 
addition, of two side counters 7; and 7, which 
could be connected in anticoincidence to M,M, 
(Fig. 3). 7, and 7: were in most cases connected 
in parallel. The whole master group was sur- 
rounded by a lead shield 10 cm thick. According 
to what was explained above, the anticoincidence 
events M,M.—T;\,2 were attributed to pene- 
trating non-shower-producing particles, such as 
mesotrons. In some very rare cases, however, 
which will be discussed below, high energy soft 
particles might also produce such an anticoin- 
cidence event. 

2. Two series of experiments were carried out. 
In the first, several analyzing counters were 
placed at different distances from the master 
group (Fig. 3a). In that way it was possible to 
obtain, from the same record of coincidences, the 
counting rate in an analyzing counter as a func- 


tion of its distance from the mesotron selector. 
In the second series of experiments, a set of 7 
analyzing counters, similar to the set described in 
experiment II, was used in order to investigate 
the penetrating power of the particles associated 
with those particles which discharge the master 
group (Fig. 3b). Since all the data obtained in 
these experiments, which were carried out in 
Chicago and at Echo Lake, are not yet analyzed 
and since similar experiments are still in progress, 
only part of the results will be reported here and 
summarized in Table III. 

In Table III the first two columns I and II are 
self-evident; column III (1, 2) gives the proper 
counting rate (counts per 10 hours) of the master 
group for the two cases where only mesotrons 
(M,M:z—T);,2) or both mesotrons and high energy 
soft particles (44, M¢2) are registered ; column IV 
gives the counting rate (also per 10 hours) of 
discharges in an analyzing counter for different 
distances D of that counter from the master 
group; column V represents the ratio »v between 
the counting rates which figure in columns IV 
and III for events M,M.—T),2. The last column 
VI gives the theoretical values for », which 


TABLE III. Experiment III. Coincidences, per 10 hours, 
occurring between (1) the mesotron selector M,M2—T)},2, 
or (2) the master group M,Mz2, and an analyzing counter K. 








I II Ill IV Vv VI 
Hours M:iM: (MiM?2 
of re- —T1,2 (1) —Ti,2)K (1) 
cord- Din (MiM:—-Ti,2)K 


or or 
ing meters MiM: (2) M.iM2K (2) oe ee —tee rth 





Chicago; Altitude: 610 ft. . 


50 (1) 2.7 x108 19.7 7.1 X1074 8.6 X10 
34 aa (2) 3.4X10* 48,3 

50 (1) 2.7 X108 8.6 3.2 x10" 4.2 x1075 
34 ” (2) 34X10 22 

50 (1) 2.7 106 5.2 1.91074 2.0 X10-5 
34 " (2) 3.4X104 13.3 

Echo Lake; Altitude: 10,618 ft. 

51 (1) 5.0 X108 15.3 3.1 X10" 1.4107 
26.5 - (2) 6.210 67.4 

76 (1) 5.0106 6.7 1.3 X10" 0.5 X107 
42 - (2) 62X10 43.8 








represent the ratio », between (a) the frequency 
of discharges in a distant counter, arising from 
collision electrons ejected in air by those meso- 
trons which traverse the selector, and (b) the 
proper counting rate of the selector. The formula 

















for this ratio is derived in paper III and is 
given by 
wy, = (BS/xD)Omax(D, EF). 


Here, B is a constant factor which is equal to 
10~ for standard air if the effective area S of 
the counter, placed at a distance D from the 
selector, is expressed in m? and D in meters; 
Omax(D, E;) is a function of D and of the minimum 
energy £; which is necessary for an electron to 
traverse the wall of the counter. Some of the 
calculated values for @max(D, Ei) are given in 
Table IV. 

In Table III the values of », were computed 
for S=2.5X10-*? m? and E,; =2 Mev. The altitude 
effect on », was taken into account (see Altitude 
Effect, in ITI). 


Discussion of Results of Experiment III 


It was found that: (a) the frequency of dis- 
charges occurring in an analyzing counter in co- 
incidence with a discharge in the mesotron 
selector is, at least for the first few meters, a 
strongly decreasing function of the distance D of 
this counter from the selector; (b) most of the 
discharges in any particular analyzing counter 
are only very occasionally accompanied by a 
discharge in other analyzing counters; and (c) the 
counting rates in these analyzing counters are 
considerably lower than those which are usually 
observed for A showers using unshielded counters 
separated by the same distance D. 

A proof that the observed phenomena, in 
which only the passage of a mesotron through 
the selector M,M2—T)\,2 is involved, are of a 
different nature from those which involve soft 
particles, is given by the experiments in which 
the side counters 7; and 7, are removed. In these 
experiments, the master group M,M, could 
record, in addition to mesotrons, also the high 
energy soft particles which emerge with a large 
shower below a lead thickness of 10 cm. Although 
in this case the proper counting rate of the master 
group increased only by ~25 percent, the corre- 
sponding increase in the counting rate of the 
analyzing counters was of the order of ~300 
percent. Moreover, the proportion of cases where 
multiple coincidences between the different ana- 
lyzing counters occurred increased at the same 
time by an even more considerable factor. 
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Because of the strong dependence of the 
counting rate of an analyzing counter upon its 
distance from the selector, the observed phe- 
nomena seem so far to be in qualitative agreement 
with the collision-electron hypothesis. But, on 
the other hand, it is seen from Table III that 
there exists a considerable discrepancy between 
the experimental and theoretical values calcu- 
lated on the basis of the collision-electron hy- 
pothesis. Therefore, it would be legitimate to 
conclude that at least part of the observed 
mesotrons belong actually to A showers. How- 
ever, there is a possible explanation, which could 
account, at least partly, for the above-mentioned 
discrepancy. 

As we have seen, the role of the mesotron 
selector is to register only those events in which 
the vertical telescope M,M, is discharged, but 
not any of the side counters 7; and 7». It is 
well known that the great majority of these 
events can be attributed to single mesotrons 

TABLE IV. Values of @max(D, E:). 








\. D meters 
E: Mev’, 1 2 5 10 
2 0.53 0.50 0.41 0.30 
4 0.42 0.40 0.34 0.27 





emerging below the lead absorber which sur- 
rounds the selector. In two cases, however, the 
selector will fail in its function as mesotron 
selector.® First, the selector will miss all meso- 
trons which emerge below the absorber with a 
collision electron discharging one of the two side 
counters 7, or 72. This effect, however, reduces 
the counting rate of the distant counters by only 
a few percent. Second, a high energy soft particle 
emerging, for example, with a narrow shower 
accompanied by a few side particles, could dis- 
charge the telescope M,M, and very occasionally 
miss both" side counters 7, and 7». Such an 
event will be recorded and considered naturally 
as being produced by a single mesotron. The 
influence of this effect on the counting rate of 
the distant counter will, on the average, be much 
more important than in the first case since a 
high energy soft particle has a much higher 


® See also Section C2 of the preliminary considerations. 

10Counter inefficiencies are too small to contribute 
appreciably to such a double missing. On the other hand, 
the anticoincidence circuit is practically 100 percent 
efficient. 
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probability than a mesotron to be accompanied 
by a shower; hence, the probability of discharg- 
ing the distant counter will, in first approxima- 
tion, also increase in proportion to this effect. 
It was emphasized before that in connecting the 
anticoincidence counters into the circuit of the 
master group, one reduces the proper counting 
rate of the master group by ~20 percent, while 
the corresponding reduction of the counting rate 
in the distant counter is then ~75 percent. 
Therefore, it is easy to understand that in case 
several percent, out of about 80 percent of the 
discharges of the master group which are at- 
tributed to mesotrons, were actually due to high 
energy soft particles, the corresponding reduction 
in the counting rate in the distant counter could 
be appreciably higher than the observed ~75 
percent. Consequently, the frequency of dis- 
charges in the distant counter, occurring simul- 
taneously with a passage of a mesotron through 
the master group, might be actually much 
smaller. It should be pointed out, however, that 
the hitherto admitted occurrence of that par- 
ticular narrow shower below a lead thickness of 
10 cm is probably a very infrequent event and, 
hence, could account only partly for the observed 
large discrepancy. 

It would be possible to obtain a more precise 
idea of the importance of the just-described 
effect—and to reduce it eventually—in using, for 
example, counters of a smaller size for the tele- 
scope M,M, and of a larger size for the side 
counters 7; and 7;. Considerable difficulty would, 
however, arise then from a very strong reduction 
of the observed counting rates which are already 
very low in the present investigations. 

This is the characteristic difficulty in counter 
experiments since it is not known whether a 
discharge in a counter was produced by a single 
particle or by several particles. Therefore, a few 
cloud-chamber photographs" of the particles 
which discharge the mesotron selector would be 
of great value in analyzing this effect. 

The altitude effect, which should be different 
for each of the two hypotheses, could throw some 
light on the problem considered here. In case of 
the collision-electron hypothesis, the observed 


a ~ 1 Valuable hotographs of tracks in a cloud chamber 
controlled by inal counters were obtained, at an altitude 
of 3200 meters, by P. Auger and J. Daudin, Comptes rendus 
209, 481 (1939). 
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effect should increase parallel to the increase in 
the intensity of the mesotron component of the 
cosmic radiation. In the alternate case, this effect 
should follow the increase in the frequency 
of the A showers. We found that between 
Chicago and Echo Lake the increase in the in- 
tensity of the mesotron component is 1.9, where- 
as the increase of the frequency of A showers of 
an average extension is ~7. The observed in- 
crease in the counting rate of a distant counter is 
~3 and lies between these two values. But 
before a significant comparison regarding the 
altitude effect can be made, it will be necessary 
to carry out further measurements in order to 
extend the common range of the distance D for 
the two altitudes and to prolong the time of 
observation. 


CONCLUSIONS 


By way of conclusion, we shall point out that: 
(a) Counter experiments are practically unable 
to give unambiguous information regarding the 
presence of mesotrons in the dense part of A 
showers. (b) Mesotrons were found to be associ- 
ated with particles of very low density, which 
could be considered either as belonging to A 
showers or as collision electrons ejected by these 
mesotrons along their path in the air. (c) Al- 
though the hypothesis of mesotrons accompanied 
by collision electrons seems to agree qualitatively 
with the experimental results so far obtained, 
the experimental values, however, are consider- 
ably higher than those given by a theoretical 
calculation for this case. Therefore, it would be 
legitimate to conclude that at least part of the 
observed mesotrons actually belong to A showers. 
But on the other hand, explanations are indi- 
cated which might account, at least partly, for 
the above-mentioned qualitative disagreement. 
Thus, only after these explanations are tested 
by experiment will it be possible to decide more 
definitely as to the presence of mesotrons in A 
showers. 

It is a great pleasure for me to express my 
gratitude to Dean A. H. Compton for placing at 
my disposal the facilities for this work, and to 
the Rockefeller Foundation for its generous sup- 
port. | am also very thankful to Professor Pierre 
Auger and to Professor Marcel Schein for their 
continued interest in the present work. 
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It was found in the preceding paper that mesotrons are 
associated with particles present at several meters distance 
from the mesotron trajectories. At least part of these 
particles are collision electrons ejected by these mesotrons 
along their path in the air. The theoretical results obtained 
in the present paper show that: (a) Only mesotrons of 
momenta p> po~10* ev /c can eject in air collision electrons 
sufficiently energetic to reach and to traverse a counter 
tube (C) placed at a distance D (a few meters) from the 
mesotron trajectory. (b) The number » of collision electrons 
arising from a single mesotron and which are able to dis- 
charge (C) is practically independent of the mesotron 
momentum ? if the influence of the spin of the mesotron 
is neglected and is given, in first approximation, by 


vo=(BS/xD) X Onax(D, E1). In this formula, if the effective 
area of (C) is expressed in m? and D in meters, then B = 10-2 
for standard air. @max is a function of D and of the minimum 
energy £; necessary for an electron to penetrate the wal] 
of the counter. For E;=2 Mev, for example, @nax decreases 
from 0.50 for D=2 m, to 0.30 for D=10 m. (€) When the 
mesotron spin (0 or 4) is taken into account, v is about 
equal to $0 for the lowest mesotron momenta, but ap- 
proaches asymptotically vo as p increases. (d) A fraction of 
these collision electrons can give, by their subsequent 
cascade multiplication in air, a certain contribution to the 
counting rate in (C), in excess of ». However, when aver- 
aged adequately over the mesotron spectrum, this shower 
effect appears as a small correction only. 





HE preceding paper,' here referred to as II, 

describes experiments (III) which show 

that mesotrons are associated with particles 

present at several meters from the mesotron 

trajectories. At least a part of these particles are 

collision electrons ejected by these mesotrons 
along their path in the air. 

The method used in these experiments consists 
of measuring the frequency of discharges oc- 
curring simultaneously (a) in a group of counter 
tubes, responding in principle only to mesotrons, 
and (b) in a counter tube located at a variable 
distance from that group. This group of counters, 
which will be called here mesotron selector or, in 
short, selector, is surrounded by a lead shield 
10 cm thick ; it consists of two counters in coinci- 
dence placed in a vertical plane, and of two side 
counters, connected in anticoincidence to the 
circuit of the two first counters. The principle of 
this selector is based on the assumption that any 
event in which the two counters in coincidence 
are discharged, but not any of the side counters, 
should be attributed to a penetrating particle, 
such as a mesotron, since a high energy soft 
particle would almost always emerge below the 
lead with a shower discharging at least one of the 
side counters. 

The purpose of the present paper is to calculate 
the mean frequency of discharges in the distant 


1 Anatole Rogozinski, Phys. Rev. 65, 291 (1944). 
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counter, which arise from collision electrons 
ejected in air by those mesotrons which traverse 
simultaneously the selector. 


A. NOTATIONS AND FORMULAS 


1. The basic formulas for collision processes 
between mesotrons and electrons will be applied 
here in the form given by B. Rossi and K. 
Greisen.? The notations used below are similar to 
those of the authors quoted. 

The velocity of light will be put c=1; the 
velocity of a particle: v=8c=8; 1 Mev will be 
taken as a unit of energy and 1 meter as a unit of 
length. In that way, the rest mass of an electron 
ue Or Of a mesotron yz will be expressed by the 
same symbol as their respective rest energies. 
These energies are assumed to have the following 
values: u,=0.5 and 1. = 100. Instead of the meso- 
tron momentum ), the ratio r= p/u=8/(1—8?)! 
will mostly be used. This ratio r will be called the 
equivalent momentum of the mesotron. 

2. Asa result of a collision between a mesotron 
of momentum ? and an electron, supposed to be 
free and initially at rest, the latter is projected at 
an angle 6 with respect to the initial trajectory of 
the mesotron and acquires an energy E’ given by 

p? cos? 6 
EB’ = 2p.——_____—_————. (1) 
[me?+(p? +)! ]— p? cos? 8 
ean) Rossi and K. Greisen, Rev. Mod. Phys. 13, 240 
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It will be seen later that the values of p to be 
taken into account for the present case will 
always be appreciably larger than yu. Therefore, 
neglecting yu,” as well as uw? in comparison with p’, 
and remembering that 2u4,=1, (1) will be ex- 
pressed by 

cos? 6 
E'= 
sin? 6+1/r?+1/p 





(2) 


3. It can be seen from (2) that, whatever the 
energy of the primary mesotron, one has always 


E’ <cot? 6. 


Thus, if, for example, @=2/4, the maximum 
energy which can be transferred to the collision 
electron by a mesotron of any energy is always 
less than 1 (Mev). 

When 1/r?>1/p, or p<y?=10', that is, for 
momenta p<10'° ev/c, (2) can be written in the 
simpler form 


cos? 6 








EF’ =— " (2’) 
sin? 6+1/r? 


On the other hand, for p>10'° ev/c, (2) be- 
comes 
cos? 6 


EF’! =— -, 
sin? 6+1/p 





(2”) 


For head-on collisions (@=0), EF’ reaches its 
maximum value E,,’, given by 


p* 
En’ =1/(1/r2?+1/p) =——, (3) 
we+p 
E,’=r°, for p<p?, (3’) 
E,'=p, for p>uw*. (3’") 


We shall use the expression (2’) since it can be 
applied to the majority of the mesotrons present 
in the cosmic radiation. It should be understood, 
however, that for p>y?, r? should eventually be 
replaced by p wherever it occurs in a formula 
which is derived originally from (2’). 


B. ASSUMPTIONS 


The following assumptions are made: (a) The 
energy loss of the mesotron is neglected over a 
distance H, which will be specified in the calcula- 


tions; (b) the mesotron trajectory is a straight 
line, regardless of any collision processes which 
occur over the same distance H,; (c) the scat- 
tering effects for both the mesotron and the 
collision electron are neglected; (d) the collision 
electrons lose energy only by ionization; the 
average ionization loss k per meter of standard 
air is taken as k=0.25 (Mev); the minimum 
energy E, necessary for the collision electron to 
penetrate the wall of a counter tube is of the 
order of a few Mev; (e) the atmosphere is 
homogeneous over the distance H, considered 
in (a). 


C. GEOMETRY; ENERGY FACTORS 


Let us consider now a mesotron of momentum 
p, directed downwards along the axis Oz. At A, in 
a plane perpendicular to Oz at O, a counter tube 
(C) is placed whose effective dimensions are 
assumed to be small in comparison with OA =D 
(see Fig. 1A). 

This counter can be discharged by a collision 
electron of energy E’ ejected at a point 7, 
located at a distance OT =H from the origin. 

The total energy loss (kD/sin 6+) suffered 
by the collision electron along the path 7A 
= D/sin @ in air, and by traversing the wall of the 
counter tube, should be such that 


a’ —(kD/sin 0+E;) >0, (4) 


or, replacing E’ by the expression (2’) and 
putting sin 6= x one has 


f(x) =x*+ax2+bx+c<0, (5) 
with 
a=kD/(E,+1), b=—(r?—E,)/r?(E, +1) 
and 


c=kD/r(E,+1). 


Collision electrons which satisfy (4) or (5) will 
be called “‘efficient”’ collision electrons. 

It can be proved that the equation f(x) =0 has 
in general 3 roots: one negative, x, and two 
positive, x2 and x3, and that x should lie in 
between the positive roots in order to satisfy the 
physical conditions of the problem. 

The general solution of (5) involves functions 
of the three parameters D, r, and E, requiring 
rather laborious calculations. It is more practical 
to proceed in the following way : 
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k= 0.25Mev/m air, 


r= p/p 
Es4Mev, 





, 1— xe 
(——) Car: 
(-----) EXE +k 




















— 





see -e J 


meee ee ww wend 

















0.5 
Fic. 1B 


Fic. 1A. A mesotron of momentum ? traveling along Oz 
ejects in T a collision electron of energy E’ which discharges 
the counter tube (C) placed at a distance OA = D from the 
mesotron trajectory. 

Fic. 1B. Families of curves E’ =cos? 6/(sin? 6+ 1/r*) and 
E*=E,+kD/sin @ as a function of x=sin @ for different 
values of the parameters r= p/u and D [Eqs. (6) and (7) ]. 
The minimum energy £: necessary for an electron to 
pone the wall of the counter (C) is taken E,;=4 Mev. 

n order to be able to discharge (C), a collision electron 
must be ejected at an angle @ such that sin @ lies between 
the respective abscissae of the two points at which E’ 
and E* intersect. 


Consider the two families of curves: 


E! =(1—x*)/(x?+1/r’) (6) 
and 
E*=E,+kD/x (7) 


for different values of the parameters r, D, and E, 
as shown in Fig. 1B. The curves show that, in 
general, EZ’ and E* intersect at two points, the 
abscissae of which are identical with the two 
positive roots, x2 and x3, of the equation f(x) =0. 
One can also see that for any values of the 
parameters D and £, there exists a critical value 
r-(D, E,) for which the corresponding curves E’ 
and E* are tangent to one another. The abscissa 
for the point of tangency gives directly the double 
root Xg=X2=X3, which f(x)=0 admits in this 
case. 

For r<r, the roots are imaginary : This means, 
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from the physical point of view, that mesotrons 
of momenta smaller than a certain critical 
momentum ~,=ur, will never project a collision 
electron which would be able to discharge the 
counter (C). On the other hand, it is easy tq see 
that, for a given £;, and for mesotrons of, say, 
r =r, collision electrons can discharge the counter 
only at distances D smaller than a certain critical 
distance D.(Ei,7%). Beyond D., no collision 
electrons arising from mesotrons of r <7 are able 
to discharge the counter (C). The critical distance 
D,. can be determined best from the two families 
of curves E’ and E*, plotted as a function of 
y=1/sin 0=1/x, 


E’ = (y?—1)/(y?/r? +1) (6’) 
and 
E*=E,+kDy (7’) 


since then the family E* is composed of straight 
lines (Fig. 2). 

Table I gives approximate values of r.(D, E;) 
for different distances D and for E,;=4 (Mev). 

According to this table, if, for example, D =5 m, 
only mesotrons of momenta p>4X 108 ev/c will 
have to be considered at all. 

It can be seen also from Figs. 1B and 2 that for 
r>r,. the corresponding values of x2 become small 
rapidly in comparison with x; and that x3 can be 
considered as independent of r for values of r only 
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Fic. 2. Same families of curves E’ and E* as in Fig. 1B, 
but plotted as a function of y=1/sin @ instead of x =sin @ 
[Eqs. (6’) and (7’)]. 
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Fic. 3. Curves of »/S expressed in m™™, as a function of r=/ for different 
values of the parameters D and £;. S is the effective area, in m*, of the counter 
tube (C). In case of vr the mesotron spin (0 or 4) is taken into account. When 
the spin term is neglected, and other less important approximations are made, 
v reduces to », which is practically independent of the mesotron momentum. 


slightly larger than r,. In first approximation, 
this value of xs can be taken as equal to its 
asymptotic value (x3)1im for r+. (x3)1im can be 
calculated easily in this case from (5), which 
becomes then : 











kD 1 
f(x) =x?+ x+ =(). (8) 
Eit1 £it+1 
For D< ~20 meters, one obtains 
[ k?D? kD 
(23) tim 1+ |- . (9) 
(E41) 8(£:4+1)) 2(£,+1) 


If D is small, (x3)jim—>(£1 +1). 

Table II gives certain values of (x3)iim, for 
different values of EZ, and D. 

From the fact that xz becomes small rapidly in 
comparison with x3, for r>r., one can show that, 
with a good approximation, 


kD 


: (10) 
r?—E, 





xe= 


The values of x2 and x3 determine the ‘‘useful”’ 
portion H,=H,—H; of the mesotron trajectory 
which, through the collision electrons ejected 


along H,, contributes to the discharges of the 
counter (C). The collision electrons ejected below 
H;=D cot 6; and above H,;,=D cot 6.=D/x. 
—(1/k)(r?—E,), have insufficient energies to 
reach. the counter and to penetrate its wall. Since 
cot 6; has a value comparable to unity and is 
practically independent of the mesotron mo- 
mentum for r>r., Hs; is practically of the same 
order of magnitude as D, for all momenta of the 
mesotron considered. On the other hand, H2 is a 
strongly increasing function of r, but it is practi- 
cally independent of D. For example, for r= 10, 
that is, for p= 10° ev/c, H,=400 meters. Hence, 
the useful portion H, of the mesotron trajectory 
is approximately given by H2. 


TABLE I. Values of r.(D, E;) for Ei: =4 (Mev). 








D meters 2 5 10 





le 3 4 6 








TABLE II. Values of (x3)jim. 











D meters 
E, (Mev) 2 5 10 
2 0.50 0.41 0.30 
4 0.40 0.34 0.27 














D. COLLISION PROBABILITIES 


The probability for a mesotron of momentum p 
and of velocity 8 to project, over a distance dH 
of its trajectory, collision electrons of energy be- 
tween E’ and E’+dE’ is given by 


Bil 
PdE'dH =— —(1—F)dE'dH. (11) 
B? E’2 


For standard air, if E’ is expressed in Mev and 
H in meters, one obtains B = 0.98 kX 10°10. 

The term F depends upon the mesotron spin 
which is not yet known. However, the results of 
experiments? on bursts occurring in an ionization 
chamber shielded by absorbers of different atomic 
number show a considerable disagreement with 
the hypothesis of mesotrons of spin 1; the agree- 
ment is much better for spin 0 or spin }. For these 
reasons, only the latter values of the mesotron 
spin will be considered here. 

One has 

F=6E'/E,,’ for spin 0 (12) 

and 


F=6E'/E,,'—E'/2(p?+u?) for spin }. 


In case of spin 3, the maximum value of the 
second term of F is equal to p*/2(p+y?)?(p? +p?) 
{see (3) ]; for p<p?, this term is always negli- 
gible ; it is equal to 3, for p=”, and approaches 3, 
for p>y?. When averaged over the mesotron 
spectrum, this term will constitute only a very 
small correction, and therefore, no distinction 
will be made in the present calculations between 
spin 0 and spin 3. 

For small values of EZ’, one obtains in all cases 
practically F=0, and P then reduces to the 
Rutherford formula 


P=B/s*E”. (13) 
E. FREQUENCY OF THE EFFICIENT 
COLLISION ELECTRONS 
(a) Approximate Formulas 


1. Let N be the mean number of mesotrons per 
unit time following the axis Oz and Ng(r)dr their 
fraction with equivalent momenta between r and 
r+dr. The function g(r) represents the differential 


3M. Schein and P. S. Gill, Rev. Mod. Phys. 11, 267 
(1939). R, F. Christy and S. Kusaka, Phys. Rev. 59, 405, 
414 (1941). S. Kusaka, Phys. Rev. 64, 256 (1943). R. E. 
Lapp, Phys. Rev. 64, 255 (1943). 
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mesotron spectrum. The integral mesotron 
spectrum will be represented by the symbol 
G(r) = So'g(r)dr. 

The differential spectrum g(r) is normalized by 


f g(r)dr=1. (14) 


2. It will be assumed first F=0. 

The mean number »y of collision electrons 
arising from a mesotron and able to discharge the 
counter tube (C) is given by 

He 
=f dH PdE’. (15) 
H3 (AE’) 

In (15), Hz and H; determine the useful portion 
H,=H,—H,; of the mesotron trajectory and are 
given, in first approximation, by (9) and (10), 
respectively. The energy difference (AE’) igs 
defined in the following way: 

Let us consider first two cones forming angles @ 
and @+d86 with respect to their common axis Oz, 
and having T as their common vertex (Fig. 1A). 
The energy interval of the collision electrons 
ejected along dH into the solid angle dw = 27 sinédé@ 
which is comprised between the two cones is 


1 dE’ 
dE’= 


Bn eg (16) 
2m sin 0 dé 


Let d be the effective diameter of the counter 
tube (C), /, its effective length, and S=d xl, its 
effective area. Since we have assumed that the 
two dimensions of the counter are small in com- 
parison with D, the values of E’, P, and dE'/dé 
can be considered as constant all over the area S$ 
which defines a solid angle Aw = 5S/(AT)? 
= § sin? 6/D*; thus, 


P dE’ 
Paz'= { — —dw 
(SB’) (ae) 27 sin 6 dé 








Ssin @ dE’ 
= —P—., (17) 
2xrD? dé 
On the other hand dH = — (D/sin? @)dé@; hence, 


He 


vpao= f dH | PdE’ 
H 


3 (AE’) 





—, (18) 


gD Pr? JSeuain cos® 8 
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where Omin=arcsin Xx~kD/(r?—E;) and @Omax 
=arcsin x3. The angle @max, as well as x3, can be 
considered as independent of r>r,. 

Developing in (18) cos*@ in a series, one 
obtains 


BS 1+r* 
dL —_— ——( 0 nex + $6.nax* _ 
mrD Fr 


Omint+-°*), (19) 


where the terms neglected in the parenthesis are 
of power 25 for @max, and of power 2 3 for Onin- 

It has been seen that except for values of r very 
close to r., the value of x2 Or Onin is small in com- 
parison with @x- Therefore, in first approxima- 
tion, vr-o is equal to 


BS 
Vo= ~—Omnex- (20) 


us 


It should be emphasized that within the limits 
of the approximations assumed, the value of v» 
which represents the contribution of a single 
mesotron to the number of discharges in the 
counter (C) is practically independent of the mo- 
mentum of that mesotron; this is true in spite of 
the fact that the useful portion of the mesotron 
trajectory increases approximately proportional 
to r’. 

3. The mean frequency n of discharges in the 
counter (C) is given thus by: 


a He 
n= f Ne(r)ar f dH PdE’ 


H3 (AE’) 
=n g(r)vpecdr. (21) 


In order to correlate m with those mesotrons 
which effectively discharge the selector, one 
should remark here that N is the frequency of the 
mesotrons which would be recorded by a hypo- 
thetical unshielded mesotron selector. In trav- 
ersing the 10-cm thickness of the lead shield 
which, in fact, covers the selector (see II), all 
mesotrons of a range not larger than 10 cm of Pb 
are stopped. Let r; be the maximum equivalent 
momentum of these stopped mesotrons. The fre- 
quency of discharges in the mesotron selector 
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Nee: Will then be equal to 
Nua=N | g(r)dr=NT1—G(rs)]. (22) 
ri 


The exact shape of the differential mesotron 
spectrum, at least on the ground, is not yet 
known at its very low end. There are, however, 
experimental reasons to believe that the fraction 
of mesotrons with energies smaller than a few 
times 10° ev is rather small. Since for 10 cm of Pb 
one has 7;~2.5, corresponding to a mesotron 
energy of 1.7 10* ev, G(r:) can be considered as 
small and probably not larger than a few percent. 

Taking into account (22), and using for » its 
value from (20), one finds that is given, in first 
approximation, by 


1—G(r.) BS 
no= Noee- 


——~ —“Faae (23) 
1—G(r;) rD 

According to Table I, the critical value r, 
represents mesotron momenta of the order of a 
few times 10° ev/c. Therefore, G(r.) is, like G(r;), 
small and the ratio [1—G(r.) ]/[1—G(r1) J~1 
—G(r-)+G(r;) is approximately equal to 1. 
Hence, 


BS 
Nno= N set —9 mez: (24) 


us 


This (approximate) formula, which gives the 
frequency of discharges in a counter of effective 
area S placed at a distance D from the mesotron 
selector whose proper frequency of discharges is 
equal to N,1, can be compared directly with the 
experimental data. 


(b) Correction for the Mesotron Spin 


In obtaining the above result the spin de- 
pendent factor F was neglected. It will now be 
taken into account (see Section D) in the case of 
spin 0, for which F=6°E’/E,,’, and it will be 
assumed that 8=1. 

Only the final result of the evaluation of 


He 
v= f dH P(i-—F)dE’ (25) 
H 


3 (AB’) 








will be given here. One finds: 


BS 1+r? 1 
t= = —— 6 nex’ =e Amin ~-h(r)} 
7 r? r 
(26) 


h(r) =arctan (7 sin @max) 
—arctan (7 sin Onin). 


For r > 30, h(r) can be considered as equal to its 
asymptotic value of 7/2. On the other hand, the 
formula (26) shows immediately that ve ap- 
proaches asymptotically vo as r increases. 

Table III gives h(r)/r for different values of 
E, and D. 


TABLE III. Values of 10?Xh(r)/r. 








Ei D 


r 
(Mev) meters 4 5 6 7 8 9 10 20 30 





2 2 24.5 22 19.7 17.6 15.9 144 133 7.2 4.9 
2 5 — 174 166 15.4 14.2 13.2 12.2 7.1 48 
4 2 21 19 18 16 15 14 13 7.1 48 
4 5 — 154 148 14 13 12.4 11.6 7.0 4.7 








In Fig. 3 the values of vo/S and vr/S are 
plotted as a function of r. 

When F is taken into account, a graphical 
integration of {g(r)vrdr will give the final result 
for m since g(r) is an empirical function. 


F. INFLUENCE OF THE SHOWER PRODUCTION 
BY THE EFFICIENT COLLISION ELECTRONS 
So far, we have considered single collision 

electrons losing energy only by ionization and 

discharging directly the counter tube (C). We 
know, however, that when the energy of such an 

electron is larger than its critical energy E.,, 

which is equal to ~ 10° ev for air, it can undergo 

a multiplication by successive radiation and pair 

production processes. As a result of this multi- 

plication, the average range of the produced 
shower is smaller than would be the range of its 
parent collision electron if it would not have 
multiplied. Therefore, in order to produce any 
effect on the counting rate in (C), im excess of v, 
the energy of the shower particles arising from an 
efficient collision electron must be high enough to 
reach and to traverse the wall of the counter, 
otherwise the above effect could become even 
negative. On the other hand, to obtain an ap- 
preciable effect, the number of these shower 
particles should be, in addition, not too low; 
therefore, the energy E’ of the collision electron 
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and the distance R= D/sin 6, which separates the 
point of its ejection on the mesotron trajectory 
from the counter tube, would have to be appreci- 
ably larger than £,,(~10*8 ev) and Apo, respec- 
tively, \o(~300 meters) representing the radia- 
tion unit of length in air. 

Let v, be the number of discharges in (C), in 
excess of v, arising indirectly, by the shower 
particles, from a collision electron ejected in air 
by a mesotron of momentum ?. Let, furthermore, 
o=v,/v be the relative shower effect. One has 
o=0 as long as E’XE,.,. On the other hand, it can 
be understood easily that provided R> ~ Xo, the 
factor o passes a negative minimum (a) when 
E’>E., remains close to E., and (b) for E’ re- 
maining in the neighborhood of the highest 
energy [see Figs. 1B, 2 and Eq. (10) ] which an 
efficient collision electron can possess and which 
is given by 1/(@min?+1/r?)~E,,’. Indeed, in both 
cases, if a multiplication of the collision electron 
takes place, the energy of any of the produced 
shower particles will be too small to reach and to 
discharge the counter. 

An upper limit of ¢ can be obtained as follows: 
According to the theory of cascade showers the 
average maximum number y of shower particles 
of energy larger than E.,, produced by an electron 
of energy E’>E.,, is given by y= E’/6E.,. If one 
considers only values of E’ for which ¥2 Wo, one 
has evidently E’ 2 Eyo’ = 6~oE.,. The above maxi- 
mum occurs for a thickness t™)o log (E’/E.,) of 
material in which the shower develops. In the 
present case only small values of @ will have to be 
considered. At any rate, 0<0,=D/Xo. For 
D<10 m, 6,~10-?. 

Thus [compare (17) ], neglecting 67, &, --~ in 
comparison with @, 


1 D cot @min 
a<- f dH P(i— F)¥dE’ 
D 


Vv cot 6s (AB’) 
1 BS Em’ 
--—___ f 6dE’. (27) 
v 2xD X6E-, ky,’ 


It can be shown that for p<1/8,?, that is, 
p<K~10'" ev/c, and consequently also for 
E,,’K1/0,?, the integral 


Em’ 
f 6dE’ <0,(Em’ — Ey,') 
E yo’ 
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and @ is smaller than a few percent. The factor ¢ 
increases with increasing p. For p>1/6,?, the 





integral 
Em’ = 
f 6dE’— —p', 
Ey,’ 2 
hence, 
1 BS ¢- 
a" pi. (28) 
vy 2nrD 12E., 


Since the mesotron momenta to be considered 
in this case are very high, » can be replaced by ». 
Now, BS/»2rD is approximately equal to 1. 


Therefore 
o<(m/12E¢:)p}. (29) 


It can be seen from (29) that ¢ will still remain 
smaller than 1 even for p= (12E.,/2)?, that is, for 
mesotron momenta as high as 1.510"! ev/c. 

It will be sufficient to remember that the fre- 
quency of the mesotron momentum in the 
differential mesotron spectrum decreases as 
p-(y>2) to make clear that the shower effect of 
the collision electrons, when averaged adequately 
over the mesotron spectrum, adds a very small 
contribution only to the frequency of discharges 
in the counter tube (C). 
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G. ALTITUDE EFFECT 


Concerning y, the only terms which in (20) or 
(26) depend upon the altitude, that is, upon the 
density of air, are B, @max, and Onin. The variation 
of @min, Which is anyway very small, can be 
neglected. The coefficient B is directly pro- 
portional to the density of air, and @nax depends 
upon this density only through the coefficient k. 
But, as can be seen from Figs. 1B or 2, @nax 
depends, at least for small distances D, very 
little upon k. For example, for E;=4 Mev and 
D=2 m, O@max would increase from ~0.40 at sea 
level (k =0.25), to ~0.45 only at the top of the 
atmosphere (k=0). Therefore, in first approxi- 
mation, one can admit that » is proportional to 
the density of air or to the atmospheric pressure. 
On the other hand, the effect of the altitude on 
the frequency n of discharges [see (21) ] in the 
counter tube (C) follows, in the general case, the 
combined effects of: (a) the altitude increase of 
the intensity of the mesotron radiation, (b) the 
altitude dependence of the mesotron spectrum 
and, (c) the decrease of v which, as has been seen 
above, is roughly proportional to the atmospheric 
pressure. 
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The solution of both the scalar and vector wave equations in regions which are bounded by 
irregular surfaces which have non-uniform physical properties has been reduced to the solution 
of a secular equation. The secular determinant is Hermitian. The solution to the secular 
equation has been expressed in a form suitable for obtaining its value to any approximation. 
Similar results are given for the corresponding eigenfunctions. Extension of these results to 
the problem of scattering and to the situation where the bounding surfaces move is indicated. 
The description of a source located in such a region is also discussed. 


I. INTRODUCTION 


HE solution to many physical problems can 

be reduced to the solution of a system of 
partial differential equations with prescribed 
boundary conditions. In the past, exact solutions 
were limited to those cases in which the boundary 
conditions were “‘simple’’ and were satisfied on 
“simple’”’ surfaces. ‘‘Simple’’ boundary conditions 


correspond to uniform physical properties of the 
surface involved. ‘‘Simple” surfaces are coordi- 
nate surfaces of coordinate systems in which the 
partial differential equations will separate. As a 
result, a great many physical problems of interest 
have not been treated theoretically. 

This paper will discuss solutions of the scalar 
and vector wave equations satisfying non-simple 














boundary conditions on non-simple surfaces. The’ 


scalar wave equation has been investigated previ- 
ously by means of a perturbation technique in 
which results were sought for conditions in which 
the variation from simple boundaries and bound- 
ary conditions were small. First-order corrections 
to the unperturbed eigenvalue have been found 
by Froelich,' Brillouin,? and Cabrera.’ First-order 
corrections to the wave functions, and second- 
order corrections to the eigenvalue in some cases, 
have been given in I.‘ Results have been given in 
the latter both for the case of discrete and con- 
tinuous eigenvalues. Applications of I have been 
made to the acoustics of irregularly-shaped rooms 
yielding good agreement with experiment.® 

The methods used here are not a continuation 
or extension of those used in I. Rather, a new 
method based on the use of Green’s functions is 
employed. Formulae are given permitting the 
calculation not only of improved first-order 
results, but also of all the higher orders. The 
problem of discrete eigenvalues is reduced to the 
solution of a secular determinant. 

The scalar and vector wave equations can be 
used to describe phenomena occurring in the 
fields of acoustics, elasticity, and electromagnetic 
theory. For example, our results can be used to 
treat the acoustics of irregularly-shaped rooms, 
scattering from irregularly-shaped objects, propa- 
gation down irregularly-shaped pipes. Moreover, 
simple boundary conditions need not be obeyed. 
For example, the walls of an irregularly-shaped 
room might be absorbing. 

The same technique can be applied to the 
determination of source functions satisfying cer- 
tain prescribed boundary conditions. These func- 
tions can be used in the solution of inhomogeneous 
wave equations. Application can be made, for 
example, to the calculation of the sound field due 
to a source placed in an irregularly-shaped room. 


Il. STATEMENT OF PROBLEM 


The equations to be discussed can be grouped 
according to type of boundary condition. The 
methods used in their solution vary in the same 


1H. Froelich, Phys. Rev. 54, 945 (1938). 

? L. Brillouin, Comptes rendus 204, 1863 (1937). 

3N. Cabrera, Comptes rendus 207, 1175 (1938). 

4H. Feshbach and A. M. Clogston, Phys. Rev. 59, 189 
(1941). This paper will be referred to as I. 

5 Bolt, Feshbach, and Clogston, J. Acous, Soc, Am. 14, 
65 (1942). 
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way although the general basic idea remains the 
same. In all cases, the equation is to be satisfied 
in a region R bounded by a surface S. The 
boundary conditions are to be satisfied on S. 
The scalat wave equation to be considered is: 


V-otk?o=0. (2.1) 


The boundary conditions satisfied by ¢ can be 
stated generally as 0¢/8n=Fg where F is a 
function which may vary on surface S. In 
acoustics, F is related to the acoustic impedance. 
The outwardly drawn normal to the surface is 
designated by n. We distinguish two cases: 


B.C.l dg/dn=Foe F small, 
B.C.II g=(1/F)(d¢/dn) 


(2.2) 
F large. 


Simple boundary conditions are limiting cases of 


(2.2). They are: 
B.C.1A 
B.C.11A g=0. 


dy/dn=0, 
(2.3) 


A general vector wave equation occurs in the 
theory of elasticity and may be written in the 
following form: 


a grad div s+ 8V’*s+x’s=0, (2.4) 


where a and £ are constants of the material. The 
vector Ss can always be split into two vectors A 
and B such that divA=0, curl B=0. In the 
latter case, a scalar potential can be defined and 
the vector equation reduced to a scalar one. In 
the first case, the vector equation can be written 


curl curl A—k*A=0 (2.5) 
where k?= x?/8. This equation is also satisfied by 
the vector potential of the electromagnetic field. 


The boundary conditions satisfied by A can be 
stated generally as 


(n XA) =(nXcurl A)-3 (2.6) 
where 8 is a dyadic. Two cases can again be 
distinguished 
B.C.1II] (mXA) 

=(nXcurl A)-3 38 small, 
B.C.I1V (nm Xcurl A) 
=(nXA)-~ small, 


(2.7) 
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where €=(3)-' and we have assumed that the 
determinant of the dyadic 3 does not vanish. In 
electromagnetic theory 3 is related to the dyadic 
impedance. The simple boundary conditions are 


B.C.1IIA 
B.C.IVA 


(nA) =0, 


(2.8) 
(nXcurl A) =0. 


Case III can be reduced to case IV and vice 
versa. Since div A=0, a vector D can be defined 
satisfying the relation A=curl D. The wave 
equation obeyed by D is the same as that for A. 
Boundary conditions II1 become (nXcurl D) 
=k?(nXD)-3. This is in the form given by 
boundary condition IV. 

All these cases will be considered separately for 
both the discrete and continuous eigenvalue 
problems. The general technique to be used 


consists of reducing equations (2.1) and (2.5) to 


their corresponding integral equations by means 
of a Green’s function. This has the advantage of 
introducing the boundary condition directly into 
the equation to be solved rather than its usual use 
as an auxiliary condition imposed on the solution. 
From these integral equations approximate 
formulae can be immediately derived. Solutions 
of the equations will then be derived by the 
method of successive approximations, and by ex- 
pansion of the solution in eigenfunctions. The 
latter method leads to a secular equation. 

Boundary condition I will be discussed in 
detail. The procedure to be used in the other 
cases will be outlined. Only specific differences 
from the treatment for boundary condition | will 
be noted. General remarks, unless otherwise 
noted, will be valid throughout. 


III. BOUNDARY CONDITION I 0¢/dn=F¢ 


The reduction of Eq. (2.1) to its equivalent 
integral equation is made possible by the use of a 
Green’s function satisfying the inhomogeneous 
equation 


VG, (x, £)+R?G;(x, £) =6(x—6&), (3.1) 


where x represents the observation point, & the 
source point. 6(x—£) is the Dirac 6 function. The 
subscript & is used to show that G; is the Green’s 
function corresponding to the wave number k. 
G, must satisfy Eq. (3.1) throughout a region 
Ro which includes the region R. Moreover, it can 


satisfy any convenient boundary condition on the 

surface So of Ro. Usually, So is a simple surface 

and the boundary conditions simple ones. 
Combining (2.1) and (3.1) we obtain: 


OG;(x, &) 
o(z)= f [on 


On; 


dy(E) 
tlie plas, (3.2) 


On: 


where the integration is performed over surface S. 
Introducing boundary condition | we have 


0G; ( ’ £) 
ols) = f oe] ——— 


On: 
— F()G,(x, ® [is (3.3) 


The result for the simple boundary condition IA 
follows immediately. 

An important special case occurs when S is a 
simple surface. Choosing a G; satisfying bound- 
ary condition IA on S we obtain: 


o(x)=— f o(t)F(8)Ge(x, dS; (3.4) 


Equation (3.3) satisfies Eq. (3.1) as can be seen 
immediately by substitution of (3.3) in (2.1). 
Thus we can conclude that this integral equation 
is equivalent to the original differential equa- 
tion and automatically includes the boundary 
condition. 

Equation (3.3) states that it is possible to pre- 
dict the required wave function in the interior of 
R by means of a proper distribution of sources on 
the surface S. This is essentially a statement of 
Huygens’ principle differing from the usual 
formulation in the use of a general Green’s func- 
tion and the direct introduction of the boundary 
conditions. Huygens’ principle has been previ- 
ously used in discussing scattering and diffraction 
problems.® 

Equations (3.3) and (3.4) can be reduced to a 
“true” integral equation involving only values of 
the unknown function on the surface S by 
noting that they can be used to predict the value 

6 See, for example: H. Lamb, Hydrodynamics (Cambridge 
University Press, 1932), p. 517; A. Sommerfeld, Die 
Differential und Integralgluchungen der Mechanik und 


Physik (Friedr. Vieweg. und Sohn, Braunschweig, 1935), 
Vol. 2, p. 853. 








of g on the surface. The resulting equation has 
a symmetric kernel and will involve one di- 
mension less than the original integral equation. 
This procedure is used in discussions of potential 
theory.’ It is particularly advantageous in those 
cases for which the boundary perturbation has a 
simple geometry, for then the equation may often 
be solved in closed form. It is also useful in two- 
dimensional problems since the resulting integral 
equation will be in one dimension and therefore 
susceptible to numerical techniques. Finally, if 
the perturbation is small compared to the wave- 
length, the kernel can be replaced by the simpler 
kernel of potential theory and the problem re- 
duced to a potential problem. Once the value of 
g on S is known, the value of ¢g in the interior can 
be immediately calculated from (3.3). 

This method is complicated in the case of (3.3) 
by the fact that ¢, being zero outside of R is, 
discontinuous at the surface S. (Cf. Appendix I. )° 
We shall, therefore, specify the value of ¢ at a 
point on the surface by approaching the point 
from the interior of R. In the case of (3.4), the 
solution ¢g has a discontinuity in slope at S. These 
discontinuities limit the rapidity of convergence 
of an expansion of ¢ in orthonormal functions. 

It is very useful to separate the integral in 
Eq. (3.3) into two parts, one of which is known. 
This is particularly true if the known part is a 
good approximation to the unknown ¢. 

Suppose the known part y, is a member of an 
orthonormal set of functions 0 satisfying (2.1) in 
a region Ro including region R and satisfying con- 
venient boundary conditions on So. G, can then 
be expressed by means of an expansion in set 0. 
The term containing y, can then be separated 
from the remaining terms of the expansion and its 
coefficient adjusted to one by means of a proper 
choice of normalization condition. Equation (3.3) 


becomes 


o(x) = a(x) + f o(E) 





OG; (x, £) 
x|~ ~ F(8) G(x, ® fs. (3.5) 


On: 
where the prime on the integral sign indicates 


70. D. Kellogg, Foundations of Potential Theory (Fred- 
erick Unger Publishing Company, New York, 1929), p. 286. 
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that the term involving y, in the surface integral 
is to be omitted. The normalization condition is 


f evnd V =1. (3.6) 
Ro 


This is stated correctly only for the case of dis- 
crete eigenvalues. It can be extended to continu- 
ous eigenvalues by integrating the integral in 
(3.6) over a small volume element in wave 
number space corresponding to the wave numbers 
associated with Wn. 

A second method employs the fact that (3.2) is 
satisfied by any solution of the wave equation 
and therefore by ¥,. We must note, however, that 
(3.2) yields yy, only inside R and yields zero 
outside R. Let us call this discontinuous function 
V,. We now find that 


Gi(x, &) 


ne 


o(x) = Va (x) + f | tt0~0d0} 


dS;. (3.7) 





Wn 
~ Gale, | Feo -—| 


ne 


Forms (3.5) and (3.7) yield exactly the same 
values for g within region R. However, they 
behave differently outside of R leading to a 
different analytical behavior at S. Expansions of 
¢ in orthonormal functions will have different 
rates of convergence depending on whether (3.5) 
or (3.7) is used. 

Our attention, so far, has been focused on the 
wave function ¢. However, the eigenvalue ? is of 
equal importance since it is needed for the 
specification of G;. Suppose the eigenvalue asso- 
ciated with y, is K,?. Using Green's theorem, it 
can be shown that 


M=K ot f[esc—vee s/f edV. (3.8) 


Introducing boundary condition | we find: 


b= K,2+ f o( "Fv. )as / f VaedV. (3.9) 
R 


A similar formula can be derived for the case of 
the continuous spectrum by performing a second 
integration in wave number space as described 
for formula (3.6). 
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Our problem has now been reduced to the 
solution of integral equation (3.3) which im- 
plicitly contains the boundary condition to be 
satisfied. Approximate formulae can be derived 
at once. If p~y,, then from (3.5) 


ox) Wal) + f val) 


dG, (x, §) 
Tm 


On: 


From (3.7), one finds 


o(x) a(x) — } Gi(x, 8) 


dW n(E) 


On: 





x| Feva(® = jes. (3.10b) 


Formulae (3.10a) and (3.10b) are equivalent as 
has been pointed out above. The eigenvalue 
becomes : 


OWn 
paKe+ [v.(—— Fv. as, (3.11) 
n 


This formula has been given previously.. The 
eigenvalue can be found to a second approxima- 
tion by using the first approximation results for 
the wave function. 

The use of these results depends upon our 
knowledge of G;. The latter is sometimes known 
in closed form, but then, the integrations are 
often unmanageable. In such case, it is useful to 
expand G, in a series in set 0 as follows: 


Wi(x)We(E) 
Gx, )=CO——. (3.12) 
t k?—-K;? 
This expansion will now be introduced into 
(3.10b) for this yields the more rapidly con- 
verging result. 


A tnWi(x) 


(3.13) 


¢(x)~V,(x) +>’ 


t 


where the prime indicates omission of the n=¢ 
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term in the summation and 


OW», 
An= t — Fy, }dS. 
fous vs) 


Expansion (3.13) can be introduced into (3.9) to 
yield a second approximation for the eigenvalue : 





(3.14) 


2 
tn 


k?=K,2+Ann/Nant >,’ ————_ (3.15) 
— 


9 


t 


where 


n= tyn . 3.1 
N [va (3.16) 


Formulae equivalent to (3.13) and (3.15) have 
been given previously in I. They have been 
applied to the acoustics of trapezoidally shaped 
rooms and agree with experiment within the 
range of their validity. 

The convergence of the series for ¢ given by 
(3.13) is fairly good, the ‘th term behaving 
asymptotically as 1/K?. This is sufficient to 
permit its introduction into (3.5) or (3.7) to yield 
a next approximation for the wave function. The 
resulting series will behave asymptotically as 
1/K,and could not be used to find either the next 
approximation to the eigenfunction ¢ or eigen- 
value k’. 

One method of improving convergence would 
be to find a function which converges asymptoti- 
cally as 1/K? and subtract it from (3.13). The 
function x satisfying 


xiao < f Go(x, £) 





avn(€) 
| Frewa(e) - is, (3.17) 


ne 


is of this type. On substitution of (3.12) into 
(3.17) one finds an expansion of the (3.13) type 
with k replaced by 0. x is the solution of Laplace's 
equation for region R, with a source distribution 
Fy,—(dv./dn). x can often be calculated in 
closed form especially for two-dimensional prob- 
lems. If this is so, the next approximation for ¢ 
can be found by introducing the amended (3.13) 
for V, in (3.7). Successive approximations could 
be found by continuing this procedure. 

These difficulties do not appear when S is a 
simple surface. It is possible to find the solution 








of the subsequent integral equation by means of 
the method of successive substitutions. This 
yields: 

Fin 


e= Vath eoKe 





FiuF en 
oe (PK (EK 
a FieFurFen | 
cen (B—KA)(B—-K2)(B-Ke) 
(Fin)? 
kK? 
FiFinFin 
(k?—K.2)(k?—K?’) 
FisFerFenFin 
(KA! —-KA(e—Ke) 








+: rie (3.18) 


k?=K,2+ Fant>’ 





+>" 





ot Paes (3.19) 


where 


F..= — [vFbds. (3.20) 


The primes indicate omission of the mth term in 
the summation over each index. It will be noted 
that the unknown k appears in these results so 
that (3.18) is really an equation which must be 
solved for k. It is, however, in a form which per- 
mits the evaluation of k? by means of successive 
approximations leading to a continued fraction 
solution for k? as compared to the usual power 
series. Its chief advantage over the latter lies in a 
considerably simplified formula.* Moreover, it is 
interesting to note that in acoustics F is a func- 
tion of k, and so it would serve no useful purpose 
to eliminate the explicit appearance of k on the 
right-hand side of (3.18). 

We shall now attempt a solution of the more 
general case of boundary condition |. It is pos- 
sible to write down general solutions of (3.5) or 
(3.7) in terms of iterated kernels by means of the 
method of successive approximations. As a rule, 
the evaluation of the multiple integrals involved 
in closed form, i.e., not by expansion techniques, 
is possible only by numerical methods. 

Another method which suggests itself would 


_ A short note giving similar results for secular perturba- 
tion theory is in preparation. 
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involve the substitution of a general expansion of 
¢g in orthonormal functions y, of set 0 on both 
sides of (3.3). The resulting linear equations for 
the unknown coefficients would yield a secular 
determinant. However, in this case this procedure 
is not valid since the g(x) occurring on the left- 
hand side of (3.3) is discontinuous at S. An ex- 
pansion of ¢g in terms of set 0 would therefore 
have poor convergence. It could not be substi- 
tuted in the right side of (3.3) since the expansion 
would not converge to the correct value at the 
surface S. This invalidates the above procedure 
for the general case. 

However, this difficulty does not occur when § 
is a simple surface. This yields the secular 
equation : 


| Fin — 5in(k? —K,,”) | =(). (3.21) 


A solution of this equation is given by formula 
(3.19) with associated wave function (3.18). 

It is also possible to derive a secular equation 
for the more general case boundary condition I. 
Two general requirements should be noted. First, 
the secular determinant should be Hermitian. 
Secondly, its convergence should be sufficient to 
yield a well-converging expansion for ¢, the éth 
term decreasing at least as strongly as 1/K?’. 

We shall expand ¢ in terms of the non- 
orthogonal set {W,}. Each member of this 
member of this set is discontinuous so that it is no 
longer necessary for the expansions to represent 
discontinuous functions. Note that this set is 
complete in region R. 

Let g= >, a:¥:. Equation (3.3) yields 


s aA in— Nin(k?—K,?) ]=0. (3.22) 


The resulting secular equation is: 
|Ain— Nin(k?—K,,”) | =(). (3.23) 
From the equation: 


An—Atn 
i.e (3.24) 
K,?-K?? 


it follows that the secular determinant is 
Hermitian. Secondly, good convergence can be 
expected since the values of the non-diagonal 
elements decrease as their distance from a 
diagonal element increases. 
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Equation (3.22) can be solved by the method 
of successive approximations for the coefficients 
a; yielding : 


g=V¥.t>d Wait , i WrisWar¥: 





tn Sin 
+ > WrirWrsWaiVit sins (3.25) 
r-AsHKrixAn 
where 
A nt— Nnilk?—K,’) 
Wiur=- : (3.26) 


Au—Nu(k?—K,’) 


The corresponding eigenvalue can be obtained 
from (3.9). 


2=K,? 
AntD WA int + WasWesA int eee 


in Sin 


+ — 
Nantd WirNint 2 WisWeiNint:*? 


in Sin 





. (3.27) 


The rapid convergence previously demanded is 
exhibited by series (3.25). All the terms involved 
in these formulae involve only surface integrals 
except the normalization integrals Nn. These 
formulae reduce in first approximation to those 
given by (3.13) and (3.15). 

The results given above apply to those situa- 
tions for which the eigenvalues are discrete. A 
great many problems, e.g., scattering, require the 
construction of a theory on which the eigenvalues 
have a continuous spectrum. The above theory 
can, with suitable modification, be applied to this 
case. Integral equation (3.3) applies to the con- 
tinuous spectrum case as well as to the discrete 
spectrum. However, in view of the usual bound- 
ary condition of incident plus outgoing scattered 
waves, forms (3.5) and (3.7) are more convenient. 

The first approximation (3.10) can be used 
directly. The usefulness of this approximation 
depends upon the proper choice of the Green’s 
function and the closeness of y, to the solution 
yg. No general rules can be made as to the latter 
since it depends mostly on physical intuition. 
However, as to the former, it can be said that G; 
should be as close as possible to the Green’s 
function for region R. This applies particularly 
over the region where the perturbation is most 
important. If the G, chosen were actually the 
Green’s function for region R, the first approxi- 
mation would yield exact results. 


In many cases, a second or higher approxima- 
tion will be needed. The formulae developed 
above can be used if all the summations are 
replaced by integrations over the continuous 
spectrum. The function y, will satisfy the 
normalization condition characteristic of continu- 
ous spectrum theory: 


J vated V=6(n—n’). 


The integrations over the continuous spectrum 
must be carried out in such a manner as to yield a 
solution satisfying the boundary conditions men- 
tioned above. This usually requires the use of 
contour integrals in the complex plane for one of 
the wave numbers. 


IV. BOUNDARY CONDITION II ¢=1/F(0¢/0n) 


As in the case of boundary condition I, this 
problem can be reduced to the integral equation : 


d¢(é) 
g(x) = f * 
é 
1 dG,(x, &) 
F(E) On; 








—Gi(x, ® es (4.1) 


The solution of this equation will be a function 
which has a discontinuous value at surface S. 
In the special case boundary condition IIA, the 
solution is discontinuous in slope at S. 

Formulae corresponding to (3.5), (3.7), (3.9) 
can be derived for this case. They are 


d¢(E) 
On: 
1 0G; (x, £) 
| 
F(E) On 





¢(x) =¥,(x) +f 





~Gilx, » is. (4.2) 
0G;(x, f) 
o(x) = ®,(x)+ f |——— 
3 


1 d¢(€) 
F(&) On; 











dS;, (4.3) 


d¢(£) -*) 


— Gi (x, | 


On: On; 











HERMAN 
ft 
x = anvn)as/ fi. eV. (4.4) 


If one assumes that gy, within the surface 
integrals of (4.2), (4.3), (4.4) one can find 
approximate formulae analogous to (3.10) and 
(3.11). 

Using expansion (3.12) we obtain the ex- 
pression : 


-_— - 


n 


Bu 
¢(x)—~vn(x) +d’ rer ys (4.5) 


where 





OWn 1 Oy: 
Bum [ ——"—y,) dS (4.6) 
On \F on 


Convergence of this series is poor. Equation 
(4.5) can no longer be substituted in (4.4) to 
obtain the next approximation for the eigen- 
value, nor can it be introduced with either (4.2) 
or (4.3) to obtain higher approximations for the 
eigenfunctions. Moreover, it is not possible to 
improve the convergence by means of an ex- 
pansion in terms of the discontinuous set {| V,}. 
For example, in the special case B.C.IIA, one 
needs an expansion in terms of functions which 
have discontinuous slope rather than value. In- 
stead of considering such a set directly, it is more 
convenient to consider the vector v= —V¢/k. 
It satisfies the vector wave equation : 


VV-v+k?*v=0. (4.7) 
Its boundary conditions are: 


= —k?/F(n-v). (4.8) 


In order to convert (4.7) into an integral 
equation, it is necessary to develop a vector 
Green’s theorem and define the proper Green’s 
function. The former is: 


fta-v(v-v)—v-0(9-0) RV 
= { [(n-u)(v-v)—(a-v)(0-uw) Ms. (4.9) 


The Green’s function to be used with a vector 
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wave equation is a dyadic satisfying 


Ve(Ve- Gu(x, &)) +h?Gi(x, &) =5(x— £)3 


(4.10) 


where 9 is the idemfactor. In acoustics, a-@ 
gives the velocity due to an elementary source 
located at £, having the direction and magnitude 
of a. &,(x, —) can be expanded in terms of the 
complete orthonormal set found by the vectors: 


u,= —Vy./K, (4.11) 
satisfying the orthogonality relation /R,Un:U,dV 
=5nm. Then: 

(4.12) 


Combining (4.7), (4.9), and (4.10) we obtain 


v(x) = f [img-v(£)Ve-Gu(x, 2) 
~V_-¥(E)mg- Gy (x, £) Sy. 


Introducing boundary condition (4.8) we obtain 
the integral equation : 


(4.13) 


2 


= { n;-v(d)] —ny- G(x, 
v(x) for w(@| ome (x, &) 
+V_-Gi(x, d fis. (4.14) 


Comparison of (4.14) and (3.3) reveals their 
strong similarity enabling us to use procedures 
similar to those of Sec. III. The first approxima- 
tion to V is 








Cas 
v~u,+ >’ —u, (4.15) 
t K,?-K?, 
where 
= {=(- k? dy, ae 
——y, 4.16 
EAP on MYO: (410 


Convergence of (4.15) is sufficient to permit its 
introduction into (4.2) and (4.4) to yield higher 
approximation for both eigenvalue and eigen- 
function. The eigenvalue result is of value: 


K,, CrtBin 
k®?~K,2+Ban/Nant >,’ — . 
K, K,?-K?, 








(4.17) 


The particular case of S a simple surface 
permits the choice of @, so that V-G,.=0 (i.e., 
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¥»=0on S). The resulting integral equation can 
be solved exactly. 


Gut 
v=u,+>! ——1, 




















t k?-K? 
GarGre 
" u;,+--- (4.18) 
rt (k?— K,*)(k? —K,’) 
where 
k? Oy, 1 dy, 
Gu=—— [ — —dS. (4.19) 
K.K; on F On 
Introducing (4.18) into (4.2) one finds 
y +5! K,.K, Gat 
9 =Yn — —y 
. oR? (k—K;) 
- KK iGarGrt y (4 20) 
v7 a + ‘ 
r,t k?(k?—K,*)(k?—K,”) , 
The value of &? is: 
p K 24G 4 x Gat Tin 
n nn} - kK? 
Gulla 
v —___- (4.21) 


fees, 
nt (k®?—K,")(k®—K,’) 


Result (4.21) can be derived from the secular 
equation satisfied by k?: 


|Gin— (k? — Kn*)bin| =0. (4.22) 


A secular determinant can be developed for the 
more general case. Vv is expanded in terms of a 
complete but non-orthogonal set {U,}. U, is a 
vector function equal to u, inside R but zero 
outside. The resulting secular determinant is: 


| Cin — (Rk? —K,.?)Min| =0 (4.23) 
where 
Car— Cin 
Mu= f wq-udV=——. (4.24) 
R K,2-K? 


This determinant is Hermitian. Solution of (4.23) 
can be given as before [cf. (3.25) ], 


v=Unt>D VarUit SY VasVerUr+--- (4.25) 
tn isn 
where 
Catr— Mar(k? —K,’) 
Vu=— ' (4.26) 





Cu — M,,(k?- K,*) 
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This permits the compufation of ¢ from (4.2): 


nt 








g=4.+> 
tn k?-—K;?? 


K, VarBre 
Lt 
swixnK, (k?—K,’) 


ve 





+ (4.27) 

Finally, the result for the eigenvalue k? has 
the same form as (3.27), with A:» replaced by 
Cin, Wate bY Vaz, and Nin by My,. Convergence 
of (4.25) and (4.26) is good. 

Remarks given in the preceding section on ex- 
tension of these results to the problem of scatter- 
ing hold here. 

V. BOUNDARY CONDITIONS III AND IV: 
(nXcurl A) =(nXA)-%. 

The reduction of wave equation (2.5) to an 
integral equation can be accomplished through 
the use of a dyadic Green’s function satisfying 
the equation 


curl; curl; G(x, £) —k®@,(x, &) =5(x—f)Y. (5.1) 


We also need a vector Green’s theorem :° 


J (A-cun curl B—B-curl curl A ]dV 


= f[A-(axcur B)—B-(nXcurl A) dS. (5.2) 
It then follows that: 


A(x) = f [A(®)-(apxcurl G,(x, £)) 


— (nz Xcurl; A(E)) - G(x, &) JdS;. (5.3) 


Introducing B.C.IV we find the integral equation 
satisfied by A(x): 


A(x) = — f (m,XA(£)) 
-Lourl, Ge(x, &)+0(€)-Gelx, MS: (5.4) 


This equation is very similar to (3.3) and per- 
mits the use of the techniques developed in 
Sec. III. 


*J. A. Stratton, Electromagnetic Theory (McGraw-Hill 
Book Company, Inc., New York, 1941), p. 250. 

















Formulae corresponding to (3.5), (3.7), and 
(3.9) can be derived. They are: 


A(x)=D,(x)— ff (mx Ao) 
*Lourl; G(x, &) +C(E) - G(x, —) JdS;, (5.5) 
A(x) =A,(x)— f {curl G(x, &) 


-(m;XA(E) —m; XD, (€) ] 
+[m;eX A(E) -¢(€) 
—(n;Xcurl; D,) ]-Gi(x, &)}dS;, (5.6) 


where D, is any solution of Eq. (2.5). A, is a 
vector equal to D, within R and equal to 0 
outside of R. If K,? is the eigenvalue correspond- 
ing to D,, then, 


W= K+ { (axa) 


-feurl D.+¢-D, WS / [ A-D.aV. (5.7) 
R 


First-approximation results can be immediately 
obtained from (5.5), (5.6), (5.7) if A(x)~D,(x) 
by replacing A(x) by D, within the surface 
integrals. 

Suppose D, is normalized so that it is a 
member of an orthonormal set. The dyadic 
Green’s function can then be expanded in terms 
of the set {D,}: 

D,.(x)D,(&) 


G(x, §) = — 2 — mas (5.8) 
t k?—K,? 





Introducing this expansion into (5.6) and re- 
placing A by D, we obtain 


Din 
ap te) (5.9) 


A(x)—=~A,, : 
(x) (x) +h x 


where 
Dim f (nXD,)-[curly Da +¢(£)-DaJdS;. (5.10) 


Asecond approximation to the eigenvalue follows: 


Bua (Din)? 
K,?-K?? 





(5.11) 


nn 
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where 


Diar—Din 
K,2-K?2 








Pu= [D,-DadV= (5.12) 
R 

Poor convergence prevents calculation of 
further approximations from (5.9). However, in 
the special case where S is a simple surface 
choose (n Xcurl G,) =0 on S (i.e., aXcurl D,=0 
on S). The resulting formulae have exactly the 
same form as (3.18). F;, should be replaced by 


m= f (axD,) -¢-DdS. (5.13) 
The corresponding secular determinant is 
[Sin — (R®?—Kn”)bin| =O. (5.14) 


A secular determinant can be developed for the 
more general case. A is expanded in the complete 
set of discontinuous vectors {A,}. The resulting 
secular determinant is 


| Din— (Rk? —K,,?)P in| =(). (5.15) 


The determinant is Hermitian. Solution of (5.15) 
can be found as in Sec. II]. Formulae (3.25) and 
(3.27) will apply to this case if ¥, is replaced by 
A:, Ant by Dut, Nin by Pen throughout. 


VI. GREEN’S FUNCTIONS 


The formulae developed in the preceding sec- 
tions apply only when there are no sources 
present in the region R being considered. This is 
usually not the case experimentally since sources 
must be present to excite the modes being in- 
vestigated, or in the scattering situation the 
source is frequently not at infinity. 

The description of the effect of a source is 
given by the Green’s function which satisfies 
the boundary conditions imposed by the physical 
situation at the surface S of R. This Green’s 
function can be constructed in those cases for 
which the boundary conditions are such that the 
solutions of the corresponding wave equation 
form an orthonormal set. The Green’s function 
can then be expanded in terms of this ortho- 
normal set according to Eqs. (3.11), (4.19), or 
(5.12). The various possible solutions of the 
wave equations yielding the desired set are given 
in the discussion above. It should be noted that 
these form an orthogonal set, but are not normal- 
ized as given. 
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However, for some boundary conditions the 
solutions are not orthogonal, or the expansion is 
unwieldy. It is then interesting to note that the 
Green’s function also satisfies an integral equa- 
tion from which approximate results for the new 
Green’s function may be deduced. We shall give 
the details for the case of boundary condition I, 
giving only the results for the others. 

Let the desired Green’s function be G;(x, &). 
Let any other satisfying the wave equation for 
other boundary conditions be L,(x, £). Then 


OL, ’ f) 
Gile, D= Late, + f | Gute —— 
Np 


dG, (x, p) 


Np 





—Li(p, f) 





es. (6.1) 


The introduction of boundary conditions | on G, 
yields 


Gu(x, )=Li(x, + f G(x, p) 


OLx(p, £) = 
x|——- F(p)Lifp, » ks. (6.2) 


On, 


For the other boundary conditions we get 








dG,(x, p) 
B.C.11 = Gy(x, &) =Li(x, +f ———— 
On, 
1 OL.(p, f) 
x. —_———L,(¢, » bs, (6.3) 
F(p) On, 
B.C.1V G(x, &) =Li(x, &) 


= f 2, xis, p) 


7 [curl, L.(p, £) +€(p) : L.(p, §) MS,. (6.4) 


Any attempt to solve these integral equations by 
means of expansions leads again to the expansion 
in eigenfunctions mentioned above. It is possible, 
however, to find an approximate value of G;, if 
G.~L, by substituting the known L, for the 
unknown G; in the right side of Eqs. (6.2), (6.3), 
(6.4). This represents the chief use of these 
integral equations in practice since usually it is 
impossible to evaluate a second approximation 
by the reintroduction of the first approximation 
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into these equations. A similar approximation 
used in a discussion of the Laplace equation has 
been given by Hadamard.'® 


VII. CONCLUSIONS 


It has been found possible to reduce the scalar 
and vector wave equations to the solutions of a 
corresponding secular determinant by means of 
the following procedure. 

1. Reduction to an integral equation through 
the use of a Green’s function. This equation 
includes the boundary conditions. 

2. For boundary conditions | and IV, expan- 
sion of the solution in terms of a discontinuous 
set of functions which vanish outside region R. 
This leads to a secular equation. 

3. For boundary conditions II and III, the 
introduction of another auxiliary quantity whose 
boundary conditions and equations are similar 
to I and IV is necessary. The technique employed 
for these boundary conditions is then valid. 

We note that the secular determinant is 
Hermitian. Its non-diagonal terms involve only 
surface integrals whereas the diagonal terms in- 
volve also volume integrals of the NV, type. 

Extension of these results to the scattering 
problem has been indicated at the end of Sec. III. 
Similar considerations as applied to source func- 
tions are given in Sec. VI. The solution of these 
equations in regions bounded by moving surfaces 
is found to the second approximation in Ap- 
pendix II. 

The procedure developed here appears to 
apply to any system of partial differential equa- 
tions whose solutions form a complete ortho- 
normal set or are orthogonal to another set. 
There is, however, a restriction for boundary 
conditions of type B.C.II. In this case, it is 
necessary for the solution of the auxiliary equa- 
tion to satisfy some orthogonality relation. This 
condition has not been found in the case of the 
Schroedinger equation. However, for B.C.1, no 
such restriction applies. In this problem, the 
solutions found for the wave equation apply 
equally well to the Schroedinger equation. It is 
this boundary condition which occurs in the 
Wigner-Seitz theory of metals. 

The author wishes to thank Professor P. M. 


10P, Levy, Legons d'Analyse Fonctionelle (Gauthier- 
Villars, Paris, 1922), p. 181. 
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Morse for his active interest and aid. He is also 
indebted to Professor J. C. Slater and Professor 
L. I. Schiff for some illuminating discussions and 
criticisms. 


APPENDIX I 


Here we shall show that the surface integral (3.2) repre- 
sents any solution of the wave equation within R but is 
zero outside. Any function f(x) which is zero outside R can 
be expanded in terms of the set 0 as follows: 


f(x) =2( J fad Vv. 


If, furthermore, f is a solution of the wave equation with 
eigenvalue k?, it follows that 


avn 
f(x) =z Sin on vin) (1.2) 


Using the expansion of the Green's function (3.12) we 
obtain Eq. (3.2). 


(1.1) 





APPENDIX II 


In this appendix, we shall consider the problem of 
moving boundary surfaces. We shall restrict the discussion 
to the scalar equation and B.C.I since the extension to 
other cases involves the same calculations given in the 
main body of this paper. The function ¢ will now satisfy 
the time dependent equation 


Vig = (1/c*) (a p/dt*). (IT.1) 


The following discussion will be restricted to cases in- 
volving nearly simple harmonic time dependence. The 
extension to a more general type of time dependence can 
be easily carried out in terms of a Fourier integral super- 
posing the various possible frequencies. We, therefore, 
consider the Green's function G;(x, £)e~**' = gx which satis- 
fies the equation 


1/2284 5(2— zeit, (11.2) 


20), = 
Vege at 


Combining (I1.1) and (11.2) we have 
0 0 
o(xjeitt= f[¢ (a : — 8k oe Jase 
~ @ 
voi [« Fett) ae 


(e) Ja Ve. 


HERMAN FESHBACH 


Another result which can be derived in the same manner is 
P(X) —iwt y(t) Age I ge wie . 
ys -f| a? ane of dSy. (IL-4) 


of ane 

It will be noted that the time enters explicitly in these sur- 
face and volume integrals as a parameter. Equation (11.3) 
will be used in the present discussion since it yields a 
value of ¢ directly, whereas (11.4) would involve further 
integrations before value of ¢ is found. However, the 
latter equation has definite advantages because it involves 
only surface integration. 

Introducing boundary conditions and time dependence 
of gx we have 


(2) =f [ (EO rercnte, & Jas; 


te) 
+1/2 fGr(x, 0[ =P + ue) ]ave. (Ls) 
The surface integral term is identical to that which occurs 
in (3.3). As a first approximation let g<y,e~**'. The 
volume integral drops out. The surface integral can be 
evaluated in two forms: 


fs 
—iwt >’ ar 


an‘ »,Anar— 1 nr(k?— K,?) 
—_ twt o_O 
sis [v.  Arr— Ner(k?—K;*) 


Formula (11.7) is taken from (3.25) and converges more 
rapidly than (II.6). Note that the matrix elements A,,, N,, 
are functions of time. It is most convenient to introduce 
(11.6) into the volume term since the y, are independent of 
time whereas ¥, in (II.7) are not. The second approxima- 
tion for the surface terms will just be the second term in 
(3.25). We finally find 


(11.6) 





v,|: (11.7) 


ene“ + z War¥-+ ¥ 3 WrsWar¥r 


; —_ 2iwA ne 
K,*)(k*—K,*) 


1 ys 
Cavin kh? — 





Neate (11.8) 


The dot notation has its usual significance of time 
derivative. These results are useful when the time de- 
pendence of ¢ is approximately simple harmonic. More 
general types of time dependence can be handled by means 
of a superposition of simple harmonic terms as suggested 
above. It would also be possible to make use of a more 
general formulation of Huygens’ principle for the time 
dependent equation. 
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To determine the amount of sputtering in a glow discharge three functions must be known: 


the number of ions of a given energy striking the cathode, the amount of cathode material 
released from the cathode by each ion, and the fraction of material released from the cathode 
which diffuses away. Expressions derived for these allow determination of the dependence of 
total rate of sputtering on the geometry of the discharge, pressure, cathode fall, current, and 
constants of the gas and cathode surface. The result is most accurate for very low voltage, 
high pressure discharges. Comparison with experimental data shows quantitative agreement 


under these conditions. 








ATHODE sputtering is the process which 
removes material from the surface of the 
cathode of a gaseous discharge and deposits this 
material as a film on the walls of the discharge 
tube or on other nearby surfaces. There is experi- 
mental evidence that sputtering is due to ionic 
bombardment of the cathode and that the ma- 
terial sputtered leaves the cathode in all direc- 
tions and in an uncharged state.' By identifying 
the spectral lines of atomic material from the 
cathode in gaseous discharges and measuring 
their intensities, Von Hippel? showed that a 
large part, if not all, of the sputtered material 
present in the region of the discharge was in an 
atomic state. There have also been numerous 
attempts to deflect the sputtered particles by 
electric and magnetic fields with the resultant 
general conclusion that they are uncharged. 
Baum® succeeded in attaching charges to sput- 
tered particles of silver and deflecting them after 
they had been extracted from the discharge 
region. He concluded that they were atoms which 
had velocities corresponding approximately to a 
temperature of 1000 degrees C. Seeliger and 
Sommermeyer* allowed a beam of ions to strike 
a silver surface in very low pressure gas and 
measured the amount of material sputtered in 
various directions, obtaining the result that the 


1In some cases chemical reactions take place between 


the gas and the cathode material, usually giving a large 
amount of cathode erosion. This type of sputtering will 
not be considered. 

* A. Von Hippel, Ann. d. Physik 80, 672 (1926). 

*T. Baum, Zeits. f. Physik 40, 686 (1926). 

*R. Seeliger and K. Sommermeyer, Zeits. f. Physik 93, 
692 (1935). 


319 






sputtered particles leave the cathode according 
to the cosine law. 

If the mechanism for sputtering indicated by 
these results is accepted, i.e., that positive ions 
bombarding the cathode dislodge atomic ma- 
terial, the theoretical problem of calculating the 
amount of sputtering may be resolved into the 
calculation of three functions: 

1. The fraction of atoms released from the 
cathode surface which find their way to other 
surfaces in the discharge. 

2. The number of ions of given energy striking 
the cathode. 

3. The number of atoms dislodged by an ion 
of given energy. 

Most previous work on sputtering has dealt 
principally with the effect of function (3) and 
with some justification because sputtering under 
high voltage low pressure conditions has been 
principally discussed. It will be shown, however, 
that under certain conditions in low voltage dis- 
charges, function (2) giving the distribution of 
ionic energies plays a dominating role. 


DIFFUSION OF SPUTTERED MATERIAL 


If the mean free path of the sputtered particles 
is much larger than the distances between the 
cathode and the surrounding surfaces, then ma- 
terial which leaves the cathode may never return ; 
but if the mean free path is much smaller than 
these distances, a large fraction of the atoms 
which leave the cathode diffuses back to the 
cathode surface. To obtain ‘quantitatively the 
fraction of material diffusing away from the 
cathode, let 
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ds,=vector element of cathode surface 5S), 
d82= vector element of surface S, surrounding 
cathode on which sputtered material is 
deposited, 
Jo=rate at which atoms leave cathode sur- 
face, 
J=rate of deposition on So, 
N=concentration of cathodic atoms between 
S; and So, 

D, v, \=diffusion coefficient, average velocity, 
and mean free path, respectively, of 
sputtered atoms in gas of discharge, and 

Ci, C2=the fraction of the atoms which adhere 
on collision with S; and S2, respectively. 

Now 


J=D [ ww-as, 
is the rate of deposition and 
v 
dienes” fas, 


gives the rate at which sputtered atoms are 
being returned to the cathode. If J» is constant 
so that a steady state is reached, 


V*°N=0 


with the boundary conditions 


v ° 
ce f Nas.=D | VN -di. 
4vs, So 


and 
Cw 


=f Nds+D [ VN -dBz. 
4 S1 S2 


If Jo>J and if atoms are emitted uniformly over 
the cathode, the boundary conditions become 


and 


Jo=iCwN,S; 


where N, and Nz are the values extant at S; 
and S2, respectively. Then N may be written 


4J, 
N= f(x, y; 2)Ni=——f (x, 3; 2) 


1010 


4DJ, 
J= ff v4-as,, 


and 





7 Ci Sy 


Ss 2 
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where f is a solution of Laplace’s equation which 
is zero on ‘Sz and unity on S;. From kinetic 
theory, D=2/8vX. Therefore, 


Pa Vf-d 
=—— | Vf-das. (1) 
2C:Si4 g, 


For the’ case where S; and S2 are two parallel 
planes separated by a distance d, the number of 
atoms deposited per second per sq. cm becomes 


1 jor 
~ 2C.d 





j 


from which it is seen that the approximation 
Jo>>J is justified when \<d. For a spherical 
cathode of radius a surrounded by a concentric 
sphere of radius ), the total rate of deposition 


aJ odd 
~ a(a—b) 


The inverse dependence of j on d for the plane 
case was found experimentally by Giinther- 
schulze’ and by Blechschmidt,* and derived theo- 
retically by von Hippel.’ It may be seen from 
these formulae that for a fixed rate of evapora- 
tion of atoms at the cathode, the actual rate of 
disappearance of material from the cathode is 
inversely proportional to the gas pressure and 
depends as well on the distance to the nearest 
surfaces where deposition may take place. A large 
part of the evaporated atoms returns to the 
cathode surface and perhaps forms the powder 
of small crystals which is often observed there. 
The evaporated atoms may, of course, condense 
in the discharge, form large aggregates, and thus 
materially change the rate of deposition, but in 
most cases the density of sputtered material in 
the discharge is too small to allow appreciable 
condensation. 


ENERGIES OF IONS STRIKING THE CATHODE 


To find the energy of a typical ion in an electric 
field E, let the ion start from rest and have a 
mean free path in the gas \. It will move in the 
direction of the electric field an average distance 


5 A. Giintherschulze, Zeits. f. Physik 38, 575 (1926). 
6 E. Blechschmidt, Ann. d. Physik 81, 999 (1926). 
7A. v. Hippel, Ann. d. Physik 81, 1043 (1926). 
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i, when it will have an energy Ed electron volts, 
and collide with a molecule having a mass 
approximately equal to that of the ion. If the 
collision takes place as between hard spheres, the 
average ionic energy immediately after collision 
is EX. Since most of the velocity before the 
second collision is in the direction of E, the 
average energy just before the second collision 
will be approximately EA+3EX. Thus, if the 
field varies, we may write for the average total 
energy at a given collision the expression 


@~Ea+3EA+---(1/2)"E,A+---, 


where E,, indicates the average value of the field 
between collisions and n+1 before the given 
collision represented by the subscript 0. If E is 
constant over several mean free paths, then 


é~ 2d. 


This value is an upper limit for é since all of the 
ionic velocity is not in the direction of E. 

To find the energy distribution function, con- 
sider variations in the energy contributions of 
only the last two free paths. This is justifiable 
because the average contribution of the next 
free path is only } the total energy, and varia- 
tions would be correspondingly small. The aver- 
age energy contribution of all but the last two 
free paths is }£\. During the last free path, 


there is a probability 


——— Aeéo 
exp (=) 
EX/EX 


that an energy between €9 and ¢€)+Aeo will be 
contributed to the total energy. Variations in 
the energy contribution ¢,; of the next to the last 
free path occur both as a result of variations in 
the free path and of variations in the energy 
retained after collision. It may be shown that 
if f is the fraction of the precollision energy which 
an ion retains after collision, then all values 
of f from 0 to 1 are equally probable if the ion 
is considered a hard sphere. The number of ions 
having energy less than ¢ becomes, then 


No e—( FE X/2) : ee 
mad | J 0 


éot+e, 
Xexp ( a ) ded fe 
Ex 








n(€) = 
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for 
e>Ed/2 
where 
No = total number of ions, 


dn(e) mo 
ene exp ( —€/EX+ 1/2) 


de EX 
i-—fs-—e 
1—exp i—( +1/2)] 

: f \EX 

<I f “|. 

0 1 —f / 
A tractable approximation for the above integral 
of form 








f e(f)df is 4¢(0) +4e(1/2)-+4e(1) 


which gives 

dn(e) no 

— ~=— exp (—e/E\+1/2) 
de EX 


9 € 
x (<+ -exp (-/-+1/2)). (2) 
8 AEX 
In order to determine EX, the gaseous dis- 
charge in which the sputtering takes place must 
be examined. If the discharge is at very low 
pressure so that the cathode dark space is large, 
E gnay actually be measured, or at least the 
cathode potential fall and length of cathode dark 
space may be measured. At higher pressures the 
cathode dark space may be too small to permit 
accurate measurement or may not even be seen. 
Considerations here will be confined to the 
common case when \ is a small fraction of D, 
the length of the dark space, and D is much less 
than the radius of curvature of the cathode 
surface. A number of experimental and theo- 
retical investigations® have shown that the elec- 
tric field in the cathode dark space of a glow 
discharge may be represented by a linear function 
which becomes zero at the edge of the negative 
glow, 
E.(D—x) 
Pp emanations 
D 


where Ec is the field strength of the cathode. 





®M. J. Druyvesteyn and F. M. Penning, Rev. Mod. 
Phys. 12, 87 (1940). L. B. Loeb, Fundamental Processes of 
Electrical Discharge in Gases (John Wiley & Sons, New 
York, 1939), p. 578. 
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On this assumption, it may be shown that 


miV 
E=2 , ; 
(1+y7)k 
where =current density at the cathode, V = cath- 
ode potential drop, y=number of electrons re- 


leased per ion impinging on the cathode, and 
k=ionic mobility. If 





A=A/p and k=K/p, 


Er=24| IC =) 
K(1+v7) 


The behavior of K with varying average energy é 
is in general complicated, but for the values of 
E/p usually encountered in the cathode dark 
space, K is thought to be inversely proportional 
to the mean ionic velocity, or 


Thus, 


A3 


Qa! 2/5 
mare [ Vi/p? }?/5 (3) 
K 


Ea=2 
(1+y7 


where the first bracket may be considered con- 
stant since y<1 and y does not change rapidly. 
For many glow discharges the cathode fall V is 
very nearly equal to the total voltage across the 
discharge, so that the above formula gives the 
dependence of the energy with which ions strike 
the cathode on the voltage, current, and pressure. 
For known discharges ¢~2E. obtained from 
this expression varies from less than 1 ev to 
several hundred ev; for the higher voltage part 
of this range the approximations made in the 
derivation are poor. 

When E) is only a few volts or less, it might be 
supposed that the work function of the cathode, 
or the energy acquired by the ion from its 
attraction to the cathode surface, is a very im- 
portant part of the ionic energy.’ Thus if the 
force of attraction of an ion to the cathode sur- 
face is simply an image force and the ion ap- 
proaches to 0.5A from the surface, the resulting 
energy is 7.2 ev. Massey® has shown, however, 


*H. S. W. Massey, Proc. Camb. Phil. Soc. 26, 386 
(1930) ; 27, 460 (1931). 
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that an ion is usually neutralized before it strikes 
a metallic surface by capturing an electron from 
the surface. A rough calculation of the type 
made by Massey shows the probability that an 
electron has been captured in an nS level when 
a 1-volt proton approaches to within a distance } 
of the surface is 





n 
mre " exp (—2b/0.54n X 10-8). 
n 2 


Considering capture in a 2S state only, the proton 
would be neutralized when it reaches a distance 
of about 4A from the surface, having acquired 
an energy due to image forces of 0.9 ev. While 
the kinetic energy acquired by the ion due to 
image forces is not known exactly, it is probably 
small and not strongly dependent on ion veloci- 
ties, so that it will be neglected. 

The neutralized ion, to be sure, has consider- 
able potential energy, but this energy will not be 
effective in releasing atoms from the surface. 
Massey’? has computed the probability per 
second of the neutralized ion becoming de- 
excited and delivering its potential energy to one 
of the conduction electrons of the cathode. He 
obtained about 10'*/sec., so that this process is 
much more likely than the emission of radiation. 
Although the rate of transfer of excitation energy 
to a lattice point is unknown, it would be less 
than the rate of transfer of energy to the con- 
duction electrons because of the larger momenta 
or shorter wave-lengths associated with the 
lattice points. It is assumed therefore that the 
excitation energy is delivered to one or more 
conduction electrons. Since the conduction elec- 
trons have a mean free path of the order of 
10-* cm, they are able to escape from the region 
of a few atomic radii about the center of dis- 
turbance without delivering energy to the lattice 
points. 

Before attempting a calculation of the third 
function, the number of atoms released per ion 
of given energy (d/o/dn), it is helpful to consider 
the expected general behavior of d/o/dn and its 
influence on the desired total rate of sputtering. 
For ionic energy ¢ less than ¢€, the energy re- 


1H. S. W. Massey, Proc. Camb. Phil. Soc. 26, 386 
(1930). 

















quired to release one atom from the cathode 
surface, dJ,/dn, will be zero. From this point 
dJ,/dn may be expected to rise rapidly at first 
and then continue to increase more slowly with 
increasing e. If Keo, then the exact shape of 
dJ,/dn has little influence on Jo, and this func- 
tion may be idealized as in Fig. 1(b). This is the 
case for very low voltage discharges and is the 
condition for which the present theory is most 
accurate. If €~«€, as in Fig. 2, both dn/de and 
dJo/dn must be known with some accuracy, and 
a quantitative treatment is difficult. For > eo, 
the idealization of Fig. 3(b) may be made. This is 
the condition under which most observations on 
sputtering have been made, i.e., in high voltage, 
low pressure discharges. 

An expression for dJo/dn may be calculated 
with the assumptions of von Hippel’s" evapora- 
tion mechanism. Suppose that an ion delivers its 
energy ¢ almost instantaneously to a small 
hemispherical region of radius ro and that this 
energy is propagated just as is heat. An approxi- 
mate solution for the temperature at any distance 
r from the center of disturbance and at time ¢ is” 


1.7X10-*'e ct 
T =—————_- exp (—r’ps/4ct)+T7T,) if —>ro?, 


ct\! ps 
(!) 
ps 


where 7)>=temperature of surface before col- 
lision, p, s,c=dénsity, specific heat, and con- 
ductivity of cathodic material. This formula may 
be used when ¢>f) where t) corresponds to the 

















Fic. 1. Ke. 


time of collision and its magnitude in seconds is 
to be chosen about equal to 7? in cm’, to give 


1A. v. Hippel, Ann. d. Physik 81, 1043 (1926). E. 
—- and A. v. Hippel, Ann. d. Physik 86, 1006 
1928). 

2H. S. Carslaw, Mathematical Theory of the Conduction 
of Heat in Solids (Macmillan Company, London, 1921), 
p. 149. 
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Fic. 2 (Above). E~eo. 
Fic. 3 (Below). €>¢€. 


the correct initial distribution of energy. If 


to=4X10-', 
e=4 ev, 

ps=0.5, 
c=0.5, 


then 
T =1700 exp [—r2/(1.6 10-5) ]+T7o. 


There is over a region of a few atomic radii and 
for a time of the order of 10-'* a very high tem- 
perature. It may be easily shown that heat 
radiation during this time is negligible and that 
no other ion is likely to strike within this area 
and time. 

The rate of evaporation of the cathode ma- 
terial may be found if the vapor pressure is 
assumed to be 


p=exp (b—a/kT), 
where a and b are constants characteristic of the 
cathode material, a being the heat of vaporiza- 


tion of one molecule. The number of molecules 
leaving unit surface area per second is 


v=C2p/mv», 


where m is the molecular mass, v the average 
velocity of molecules in the vapor state, and Cy 
the fraction of molecules which stick to the 
surface on striking it. Hence 


C: exp (b—a/kT) 


2 4 
2(=mir ) 
TT 





v 
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To find the total number of molecules evaporated 
by collision of a single ion with the surface, the 
maximum rate of evaporation per unit area will 
be multiplied by the length of time At required 
for the temperature to fall to such a value that 
the rate of evaporation is one-half its maximum, 
and by the area AS over which the rate of 
evaporation is not less than one-half maximum. 
This yields, on neglecting 7», 


dJy 1.5X10-*Cye! 


dn a°c*!4(psto)'m! 





a (clo) ; 


- | 4) 
2.3 10-e(ps)! 





Xexp |>- 


The usefulness of this expression is limited by an 
inexact knowledge of several quantities involved 
in addition to the assumptions under which it is 
derived. The thermal constants s and c may have 
values considerably different from those meas- 
ured under ordinary conditions. The proper value 
of to determines the radius over which the kinetic 
energy of the impinging ion is assumed to be 
distributed by the collision. Substituting reason- 
able values 


ps=1.0 cal./degree; c=0.5 cal./sec. cm? de- 
gree ; f9=3X10~"* sec. ; a=6X10-" erg (€9 =3.8 
ev); b=43; m=10-" g; C.=1. 


The above formula gives 


dJ» dn 





4 ev | 2x 10-* atom /ion 
6 | 0.2 
8 | 2 


| 


This shows that the above calculation gives a 
correct order of magnitude and demonstrates 
the rapid variation of d/)/dn near the energy 
required to remove one atom from the surface. 
The total rate at which atoms leave the 
cathode surface may be found by integrating 
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- dJ» dn 
o= —— ——ae 
€0 dn de 
IG(EX)! c® 79 x 
-—— | (<+-—exp (—x+1/2) xt 
1+y¥ o/EX\8 4 
H 
Xexp (-*-—+1/2)as 
Edyx 
where 
1.5 10-*Ce 
G= . 
a*c5/4( psto) tm! 
a(cto)! 
H= —, 
2.3 10-*7(ps)! 
I 
———=jonic current in ions/sec. 
1+y¥ 


This integral can be expressed in terms of error 
functions, but the result is complicated, and 
sufficient accuracy may be obtained without the 
complete expression. For the low pressure high 
voltage case where EA> €o, 


21(re)'G I 


For the opposite case, where EAKe, 


Gleo! €0 H €0 7 
Jo~ exp ( - -—+1/2)(—+-). (6) 
i+y EX € 4EX 4 


As was expected, the expression for Ed>e 
depends entirely on the function giving the 
number of atoms knocked off the cathode per 
impinging ion, and not on the distribution of 
ionic energies, whereas the rapidly varying ex- 
ponential part of the expression for EAKe, 
depends entirely on the energy distribution func- 
tion. Since dn(e)/de is much better known than 
dJo/dn, the second expression should be corre- 
spondingly more accurate than the first. 
Using (1), (3), (5), and (6), we have 





if Ed> € 
Als Vivy* 
Ja—{—) (7) 
Se pb \p? 
if EXKeEo 


BI p? 2/5 C€o p? 2/5 7 
J=— exp (—ceo| —] (4—] +-), (8) 
p Vi 4LVi 4 
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Fic. 4. Neon argon discharge. Ne, 99 percent; A, 1 
percent. Rate of sputtering per unit current curve: 
J/I=const. [7+ (20.5/J*/*)] exp [—(20.5/J?/*)], I=total 
current in milliamp., J=rate of sputtering measured as 
1000 /tube life in hours. 


where A and B may be roughly determined from 
geometry of the discharge and from constants of 
the gas and cathode material, 


(+)KPs 
gee 
2rA' 


and J is the total discharge current. Although 
formula (7) agrees qualitatively with experi- 
mental results, it is formula (8) which is of par- 
ticular interest here and will be compared in 
detail with experiment. 

Sputtering measurements in the region EAKeo 
have been very meager because the rate of sput- 
tering in this case is extremely slow. Recently, 
however, Rockwood" has measured these slow 
rates of sputtering by measuring the rate of 
disappearance of the thin layer of active material 
on the cathodes of commercial cold cathode gas 
tubes. He showed that the aging of cold cathode 
tubes is due to sputtering of active material from 
the cathode, and by measuring the life of a 


8G. H. Rockwood, A.I.E.E. Technical Paper (May, 
1941), pp. 4-117. 
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number of tubes obtained the points in Fig. 4 
for sputtering rate vs. current. A continuous 
fixed current was passed through each tube until 
it showed a characteristic increase in voltage 
drop. This increase in voltage marked the end 
of useful tube life and the disappearance of a 
large part of the cathode coating. Each point on 
Fig. 4 represents the mean of three tubes. These 
tubes contain neon gas at 60 mm Hg with a small 
admixture of argon and operate at 60 volts with 
constant cathode area (2.7 cm?) over a range of 
current from about 5 to 200 milliamperes, making 
good the approximations involved in deriving (8). 
Agreement between theory and experiment is 
obtained as shown in Fig. 4 by properly choosing 
two constants to fit an expression of the form 


J=al~'?"(7+b/ I?!) 


to the experimental points. 
Using the experimentally determined value of 
b and comparing its theoretical expression 


2/5 
[I 
“1 


we may evaluate ¢. Since it is thought that 
most of the ions striking the cathode are argon," 
moving in a gas which is essentially neon, we 
estimate from Table I, the values 

TABLE I. Value of constants A and K where A=Ap 
(A=mean free path in cm), p=pressure in mm Hg, and 


K=kp,/e (k=mobility of potassium ion in electrostatic 
units, €=ionic energy in electron stat-volts). 

















* E. H. Kennard, Kinetic Theory of Gases (McGraw-Hill Book Com- 
pany, Inc., New York. 1938). 
** A.V. Hershey, Phys. Rev. 56, 908 (1939). 


A=7.0X10™, 
K=7.0X10', 
c=8.3X 10° 


and 

so that 

if p is measured in mm of Hg, V and ¢ in static 
units. This gives ¢)=0.013 electron stat-volts or 


4M. J. Druyvesteyn and F. M. Penning, Rev. Mod. 
Phys. 12, 87 (1940). 








TABLE II. Sputtering data. 











Rate of 
Discharge Discharge sputtering 
pressure current P 
p arbitrary 
mm Hg milliamperes units I/p* J/p 
5 2.6 <4.0 0.104 <0.8 
10 11.5 5.8 0.115 0.58 
15 27 14.2 0.120 0.95 
20 47 20 0.118 1.0 








3.9 electron volts which is in agreement with the 
observed values of ¢)9 for metals varying from 
about 1 to 10 ev, and in particular with 2.1 ev 
given for pure Ba or Sr."® The cathode material 
is a mixture of Ba, Sr, and their oxides. 

The second constant which is determined by 
the experimental data may be expressed in terms 
of the average number of atoms released by 
each ion of energy greater than é9. A tube operat- 
ing at 25 milliamperes has a lifetime of 400 hours, 
during which time 1.8 X 10* ions will have struck 
the cathode, and of these ions 3.610” will 
have energy greater than 3.9 ev [Eq. (2) ]. The 
cathode coating consists of about 1.0 mg of 
material of average atomic weight 100, and if 
we suppose that } of this must be sputtered off 
to end the useful life of the tube, the number of 
sputtered atoms is found to be 1.5X10'%. The 
shape of the cathode and tube walls allows no 
simple evaluation of the fraction of atoms re- 
leased from the cathode which are deposited on 
the tube walls, but a rough calculation from 
formula (1) gives 10-*, so that about 1.510” 
atoms are evaporated by the cathode or 0.4 atom 
per ion of energy greater than 3.9 ev. The fraction 
0.4 seems physically reasonable and also agrees 
with the values found above from expression (4) 
for dJo/dn. 

Rockwood also made measurements on a series 
of tubes containing pure argon at different 
pressures and operating at currents which were 
nearly proportional to the square of the pressure, 
so that the expression (p?/ Vi) would remain con- 
stant.!® In this case, it is easily seen that accord- 
ing to Eq. (8), the rate of sputtering will be 
simply proportional to pressure. Table II shows 
that this is true within experimental error. This 


%F, R. Bichowsky and F. D. Rossini, The Thermo- 
chemistry of the Chemical Substances (Reinhold Publishing 
Corporation, New York, 1936). 

16 Private communication. 
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fact is not a test of the specific form of Eq. (8), 
for by making ?/Vi constant, equivalent dis- 
charge conditions are obtained with the same 
ionic energies, so that it is only a test of funda- 
mental assumptions and of Eq. (1) for the 
fraction of evaporated material which escapes 
from the cathode. 

The rate of sputtering in a pure gas is compli- 
cated by the fact that ions exchange their charges 
rapidly in the neutral atoms of their own kind. 
This phenomenon probably does not occur in 
the mixture of gases cited above because the 
argon ions carry a large part of the discharge 
current, and argon is present in such small 
quantities that exchange of charge between an 
argon ion and an argon atom may be shown to 
be quite improbable. The effect of such an ex- 
change would be to decrease the ionic mean 
free path and to make it strongly dependent on 
ionic velocity. Kallman and Rosen" show that 
the mean free path is decreased by a factor of 
two by the exchange of charge in He. The results 
quoted above for pure argon with varying pres- 
sures would not be affected because (p?/ Vi) and 
hence the ionic energies and free paths were 
kept constant. However, since the free path 
probably changes rapidly with change in ionic 





I 


Fic. 5. Pure argon discharge. 100 percent A. Rate of sput- 


tering per unit current curve: J/I =const. [7+ (30.6/J*)] 
xXexp [— (30.6/J?/5)], I=total current in _ milliamp., 
J=rate of sputtering measured as 1000/tube life in hours. 


17H. Kallman and B. Rosen, Zeits. f. Physik 64, 806 


(1930). 
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energy, formula (8) will not be accurate under 
more general conditions. 

With this reservation, it is interesting to com- 
pare measurements made by Rockwood'* on 
variation of sputtering with current for argon- 
filled tubes. His points are fitted to Eq. (8) in 
Fig. 5. Treating the experimentally determined 


18 Private communication. 
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constants as before, and assuming that the ionic 
mean free path is just half the atomic mean free 
path because of electron exchange and that 
consequently the mobility of an argon ion in 
argon gas is half that of a potassium ion in 
argon we obtain ¢9=3.1 ev and dJo/dn=0.27. 
These quantities agree as well as can be expected 
with their more accurate determination from the 
neon-argon discharge. 





PHYSICAL REVIEW 


VOLUME 65, NUMBERS 11 AND 12 


JUNE 1 AND 15, 1944 


On the Extraction of Electrons from a Metal Surface by Ions and Metastable Atoms* 
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An elementary theory is developed for the process by which a helium ion or a 24S metastable 
helium atom may extract an electron from a metal surface. The helium ion after collision may 
become a neutral atom although sometimes an excited or even a metastable atom. In the case 
of the metastable atom the collision involves an exchange of electrons between the metal and 
the atom, the excess energy being carried away by the ejected L shell electron. The average 
distance of transition § for He* and for the 2*S metastable helium atom is calculated for 
different velocities. For v=105 cm/sec., &met=2@o and Sion = 11.5¢9 where ap =0.528 X 10-* cm. 
For v=3X 107 cm/sec., &met comes out less than 0.1a@9 and Sion = 6.0ao. These values are different 
from the values obtained by Massey in a similar computation. The theory accounts for some 
but not all of the experimental evidence. The probable shape of the potential that the metal 


surface offers to a metastable atom is found. 


I. INTRODUCTION 


POSITIVE ion or a metastable atom col- 
liding with a metal surface may undergo a 
collision of the second kind with the aid of one 
of the metallic electrons. In the case of a meta- 
stable atom whose radiative transition to the 
ground state is forbidden by the spin conserva- 
tion selection rule AS =0, the collision can involve 
an exchange of electrons between the metal and 
the atom which leaves the atom in a state of 
different multiplicity—let us say, the ground 
state. The positive ion simply captures the 
* Submitted by Amador Cobas in partial fulfillment of 
the requirements for the degree of Doctor of Philosophy, 
in the Faculty of Pure Science, Columbia University, New 
York, New York. Publication assisted by the Ernest 
Kempton Adams Fund for Physical Research of Columbia 
Universiy. 


** On leave from the University of Puerto Rico. Guggen- 
heim Fellow 1943-44. 


metallic electron, thereby becoming a neutral 
atom, although sometimes an excited or even a 
metastable atom. In the latter case, the meta- 
stable atom may decay with secondary electron 
emission. 

The electron emission due to impact of mercury 
metastable atoms on a nickel surface was first 
definitely proved by Webb.' This was followed up 
by the work of Messenger? and Coulliette.’ 
Sonkin‘ studied this electron emission using a 
tungsten surface under various conditions. 
Oliphant® also studied the emission of electrons 
from metal surfaces as a result of impact by 
helium metastable atoms. Besides the above- 


1H. W. Webb, Phys. Rev. 24, 113 (1924). 

?H. A. Messenger, Phys. Rev. 28, 962 (1926). 

*H. J. Coulliette, Phys. Rev. 32, 636 (1928). 

*S. Sonkin, Phys. Rev. 43, 788 (1933). 

*M. L. E. Oliphant, Proc. Roy. Soc. Al24, 228 (1929). 
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mentioned direct studies, the extraction of elec- 
trons from metal surfaces by metastable atoms 
has been suggested as a mechanism to ac- 
count for various phenomena. Uyterhoeven and 
Harrington® made use of the action of meta- 
stable atoms on metal surfaces to explain the 
emission of electrons from cold electrodes in dis- 
charges of neon, argon, and helium. In a similar 
manner Found and Langmuir’ tried to explain 
the secondary emission from a negative probe in 
a neon discharge. Other authors*-” considered 
the action of metastable atoms on metal surfaces 
in their discussion of electrical and resonance 
radiation phenomena in gases, making use of 
some of the before-mentioned investigations. 
Penning" investigated the electron emission pro- 
duced by positive neon ions which pass down a 
fine metal canal. Oliphant’ pointed out that 
Penning’s results probably include a contribution 
from metastable atoms produced as in his work. 

The extraction of electrons by positive ion 
bombardment has been studied by many more 
people than the case of the metastable atoms. We 
shall mention only those where the energy of the 
bombarding ions is of the order of a few hundred 
volts. The reason for this limitation is that at 
higher energies the positive ion is able to eject 
secondaries from the metal by virtue of its 
motion besides merely capturing electrons due to 
its electric field. Oliphant" studied the electron 
emission from a metal surface due to He*+ bom- 
bardment. Chaudhri" investigated the secondary 
emission from a nickel surface hy Hg*+ bombard- 
ment. He admitted the possibility that metastable 
mercury atoms could be responsible for part of 
the emission. Ramsauer and Kollath" studied the 


*W. Uyterhoeven and M. C. Harrington, Phys. Rev. 
35, 438 (1930) ; 36, 709 (1930). 

7C. Found, Phys. Rev. 34, 1625 (1929); I. Langmuir 
and C. Found, Phys. Rev. 36, 604 (1930); C. Found and 
I. Langmuir, Phys. Rev. 39, 237 (1930). 

8 K. K. Darrow, Electrical Phenomena in Gases (Williams 
and Wilkins, 1932). 

*L. B. Loeb, Fundamental Processes of Electrical Dis- 
—_ in Gases (John Wiley & Sons, Inc., New York, 
1939). 

0A, C. G. Mitchell and M. W. Zemansky, Resonance 
Radiation and Excited Atoms (The Macmillan Company, 
New York, 1934). 

11F, M. Penning, Vers. K. Ak. Amst. 36 (1927); Physics 
8, 13 (1928). 

2M. L. E. Oliphant, Proc. Roy. Soc. A127, 373 (1930). 

18 R. M. Chaudhri, Proc. Camb. Phil. Soc. 28, 349 (1932). 
( oa Ramsauer and Kollath, Ann. d. Physik 16, 560 
1933). 
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reflection of protons from solid bodies. Brasefield!5 
measured the efficiency of argon ions and argon 
metastables in extracting electrons from a metal 
surface. Healea and Chaffee" studied the electron 
emission from a hot nickel target bombarded 
with molecular hydrogen ions. Veith’? measured 
the ion reflection and the electron emission while 
bombarding metals with K+. Guntherschulze and 
Bar'® studied the emission of electrons from a 
MgO cathode in the glow discharge using He, Ne, 
A, He, Ne, and Os. Healea and Houtermans!® 
measured the electron emission and the ion 
reflection from a degassed nickel target bom- 
barded by He, Ne, and A ions. Korff and Presen t® 
made use of the transition of ions into metastable 
atoms near a metal surface in their discussion of 
the discharge mechanism of Geiger counters. The 
energy involved in this process is thermal energy, 

Oliphant and Moon*! explained the neutraliza- 
tion of positive ions during a collision with a 
metal surface by the auto-electronic emission 
into a state of equal energy in the atom under the 
influence of the electrostatic field of the ap- 
proaching ion. This theory was admittedly not 
made very precise. Massey” made a quantum 
mechanical calculation of the neutralization of a 
positive hydrogen ion when colliding with a metal 
surface and found the average distance at which 
this transition takes place. He assumed that for 
the case of a He* this distance is larger by a 
factor of two or more. He also made a similar 
calculation for the emission of electrons during 
the impact of metastable helium atoms and a 
metal surface. When the helium ion is neutralized, 
it may assume the triplet 2°S state. He used 
perturbation theory to find the probability of a 
transition as a function of the distance from the 
metal surface and obtained the result that the 
probability of transition for the positive ions 
approached unity at a smaller distance from the 

146 C, J. Brasefield, Phys. Rev. 44, 1002 (1933). 

16M. Healea and E. L. Chaffee, Phys. Rev. 49, 925 
(1936). 

17 W. Veith, Ann. d. Physik 29, 189 (1937). 

18 A, Guntherschulze and W. Bar, Zeits. f. Physik 107, 
11-12, 730 (1937); Zeits. f. Physik 108, 11-12, 780 (1938). 

aan Healea and C. Houtermans, Phys. Rev. 58, 608 
eS A. Korff and R. D. Present, Phys. Rev. 65, 253 


(1944). 
1M. L. E. Oliphant and P. B. Moon, Proc. Roy. Soc. 
A127, 388 (1930). 
2H. S. W. Massey, Proc. Camb. Phil. Soc. 26, 386 
(1930) ; 27, 460 (1931). 
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metal surface than for the case of the metastable 
atoms. This means that if helium ions and meta- 
stable 2°S helium atoms are approaching a metal 
surface, the average distance at which the meta- 
stable atoms will decay is larger than the average 
distance at which the ions will be neutralized. 


II. METASTABLE CASE 


In this paper we shall show that in Massey’s 
analysis of the metastable case the wrong matrix 
element was taken for the calculation of the decay 
constant. Taking the correct matrix element, one 
obtains the more reasonable result that the 


1*19.77v 


1s 


Fic. 1. Energy balance diagram for the extraction of 
electrons by metastable helium atoms from molybdenum. 
The region to the left of AB represents the potential well 
equivalent to the metal surface; the lines represent the 
continuous energy levels of the metallic electrons, the 
spacing of the lines roughly representing the density of 
the levels. The region to the right of AB represents the 
levels of the atomic electrons. The thick arrows indicate 
the transition which involves an exchange of electrons. An 
empty circle represents a vacancy in one of the atomic 
energy levels. A black dot represents an occupied level. 


average distance for transition is larger in the 
case of the ion than in the case of the metastable 
atom. We proceed to develop our theory for the 
case of the triplet 2*S metastable helium atom 
colliding with a metal surface. 


a. Energy Balance 


Figure 1 shows the energy levels of the helium 
electrons, the metallic electrons, and the final 
free electron for the case of a Mo surface. The 
thick arrows indicate the transition which in- 
volves an exchange of electrons. The lines to the 
left of AB represent the continuous energy levels 
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of the metallic electrons, the spacing of the lines 
roughly representing the density of levels. An 
empty circle represents a vacancy in one of the 
energy levels of the helium atom. A black dot 
represents an occupied level. Equating initial and 
final energies, we get for the transition shown 


—M+Erm—e= —(1+6)+ Fi, 
E.=I-—-M+Ekxo, 
where 


I=energy difference between the 2°S and the 
1S states, 
e=energy of the 2°S state, 
Exo=kinetic energy of the metallic electron, 
E,=energy of the free electron after 
transition, 
M =energy of the potential well which represents 
the metal. 


the 


From this equation and the values shown in 
Fig. 1 it is clear that any metallic electron is 
capable of making the transition in the case of a 
Mo surface. 


b. Decay Rate 


From the quantum mechanical theory of col- 
lision, one finds that the rate at which metastable 
atoms decay as a function of the distance s from 
the metal surface is 


d(s) =[nmkQ?/xh? || H\?, 
where 


m = mass of the electron, 

h=h/2x=Bohr unit of angular momentum, 

k =(2mE)*/h, where E is the kinetic energy of 
the emitted electron, 

n=density of free electrons, 

Q=volume used for normalizing 
functions, 


the wave 


and H is the matrix element of the Hamiltonian 
for the transition where the initial state is a 
triplet 2°S metastable helium atom and a free 
electron of momentum /Ako inside the metal and 
the final state is a helium atom in the ground 
state and a free electron of positive total energy 
and momentum Ak. 




















A. 








c. Matrix Element 





The matrix element is the following: 


H= ff fusryur(rdlexp (—ik-r,) ]/V2 





2e7_ e? 
x]. 1s(7%2)M28(P1) 


r3 T13 





xX [exp (tko- 3) ]//2 dr,drod73, (1) 


where the u’s represent hydrogenic wave func- 
tions and the subscripts 1, 2, and 3 attached to 
the r’s refer to the three electrons involved in the 
transition. e is the electronic charge. 
U1s(13)Ure(%2) = (1/4/m)(Z'/ao)! exp (—Z'rs/ao) 
X (1/4/23) (Z'/ao)! exp (—Z’r2/ao) 
=([a* exp a(—re—rs) |/a, 
U1s(P2)U2e(71) = (1/4/4)(Z1/a0)! exp (—Zir2/ao) 
> 4 (1/4/22) (Z2/ao)*(2 —Zor;/ao) 
Xexp (—Zer1/2a0) =[(bc)!/x ] 
X(1—r,) exp (—bri—cro), 













where 
ay=one Bohr radius=0.528 X 10-* cm, 
a=Z'/ao=1.69/do, 
b=Z2/ao=1.19/ao, 
c=Z,;/a9=2.04/ao. 


The numerical values assigned to Z’, Z;, and Z2 
are those used by Bethe.” 

This matrix element has the physical interpre- 
tation that in the initial state electron 3 is a metal 
electron, and electrons 1 and 2 are the atomic 
electrons; and that in the final state, electron 1 is 
a free electron, and electrons 2 and 3 are the 
atomic electrons. 

Massey” starts with an expression which is 
identical to Eq. (1), but in his subsequent calcu- 
lations he actually uses the following matrix 
element: 


H= ff furcrdustrdCexp (—tk-r,) ]//2 


2e?— e? 
x [== furs) 


Ts 118 
x [exp tky +13 |/+/Qdrdredr3. (2) 
%H. Bethe, Handbuch der Physik, Vol. 24/1, p. 365. 
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This matrix element implies that electron 3 has 
initially both the wave function exp iko-r3//Q 
and the wave function wu2,(r3) and that simi- 
larly electron 1 has finally the wave function 
[exp (—7k-r) ]/\/Qand the wave function 1;,(r;). 

The,contribution from 2e?/r; in Eq. (1) is not 
very small because u2,(7;) is not orthogonal to 
[exp (—tk-r,) ]/(Q)4. If we were to use a more 
exact wave function for the free electron, the 
orthogonality would be better, and for this 
reason we neglect this contribution. 

If we substitute for the u’s the hydrogenic 
wave functions they represent, Eq. (1) becomes: 


H=Af f [Cexpa(—r—re) thers) 


all space 
X[1/ris }(1— bri) 
Xexp (—bri—cre+iko-r3)dridredr3. (3) 
A =[e*a*(bc)#]/r°Q. 


Now we undertake the evaluation of H. The 
integration with respect to drz can be carried out 
immediately. This gives a factor 8x/(a+c)*. The 
remaining integral 


J [texe (—ar3—ik-r;) | 


all space 
X(1/|r1—r1s| ](1—4r1) 
Xexp (—britiko-rs)dridrs 3’) 


is rather complicated, but one can get a fairly 
good approximation, valid for large distances s 
of the atom from the metal surface; s being large 
compared to an atomic radius. Then exp (—ars) 
falls off much quicker than 1/|r,—r3|, and one 
can substitute 1/|r,—s| for the expression 
1/|ri1—rs| ° 

The integral with respect to dr; is the following: 


2r r/2 x 
{ f [ [exp tkors(cos 83 cos Oko 
0 0 o/s sec 03 


+sin 03 sin 09 Cos $3) — ars |rs*dr3 sin 03403ds, 


where we have taken s as the polar axis and 
whére the integration over 63 is taken from 0 to 
m/2 to take care of the fact that the wave 
function of the metallic electron is zero outside 
the metal. The integration with respect to $3 
gives 2nJo(kors sin 03 sin 0x9). Because of the 
overlapping of the wave functions, the integrand 
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will have non-vanishing values only for 6; small. 
Therefore we can replace the Bessel function by 
unity and do the rest of the integration im- 
mediately obtaining 


2x[exp (iko cos Ok9—a)s |/as 
X[—s?/(iko cos 0k9—a) 
+2s/(iko cos 09 —a)?—2/(iko cos 6ko—a)* }. 


The square of the absolute value of this expres- 
sion is a factor in | H|*, and we are interested in 
knowing how it depends on 6x9. For s =3a» and 
§k9 = 2/2 this absolute value squared is bigger by 
a factor of 1.44 than for the same s, but 6%. =0. 
For s = 4d» the ratio is 1.40. If we assume kp to be 
zero, i.e., if we take the initial metallic electron to 
be at rest inside the metal, the integral becomes 
just 


[ex (—ar3)dr3 
=[22/as |(as+2) exp (—as). (4) 


The square of this integral for s = 3a» gives a value 
1.32 times the value obtained when we assume ko 


" to be perpendicular to the metal surface. Thus it 


is reasonable to believe that it does not matter 
very much whether the initial electron is at rest 
or moving in any direction whatsoever. The case 
of the electron moving perpendicularly to the 
metal surface but away from the surface gives a 
result identical to the case where the electron is 
moving towards the metal surface. 

We now consider the integration with respect 
to dr;. The integral we have to do is the following : 


[Sf forinst10irs 


x [exp (—br; —ik-r;) |ri2dry sin A; de, dd). 


We can expand 1/!r,;—s| in terms of spherical 
harmonics and exp (—7k-r,) in terms of spherical 
harmonics and Bessel functions. The integration 
over 6, and ¢; can be done immediately. 

The integration over r; reduces to 


f DX (—1)'Pi(cos &)[r1'+?/s'**] 
0 3 
Xexp (—bra)[1 — bry \(m/2kri)§T i44(kriddrs 
+f X (-d'Pi(cos o)Es!/ns1 
8 l 
Xexp (—bri)[1—bri (a /2kri)*J 4(Rriddri. 


The integration with respect to dr; can be 
evaluated numerically for arbitrary values of s. 
For large s, one can obtain a good approximation 
by considering only the first integral and taking 
the limits of integration over r; from 0 to © 
instead of from 0 to s. 

We assume k to make an angle 6 with the 
perpendicular to the metal surface and estimate 
its numerical value by letting E, = (19.77 —4.3) 
ev, where 19.77 ev™ is the difference in energy 
between the ground state and the 2°S metastable 
state of helium, and 4.3 ev®® is the approximate 
value of the work function of Mo, which is one of 
the metals Oliphant used. Neglecting all terms 
for which />2, we get the following result: 


(1.077 K 10-**) /s —7(8.519 XK 10-8 cos 6,.) /s? 
+ (6.520 x 10-*# — 1.955 K 10-*! cos? 0) /s°. 
Multiplying this expression by its conjugate, 
integrating over all space, and putting together 


our previous calculations, we find that the rate at 
which metastable helium atoms will decay is: 


\(s) = 8nmkQ*A 2x?/h®at[8x/(a+c)* } 

Xexp (—2as)[(3.956 X 10-*+-0.796s)? 

+ (3.268 X 10-*+-0.404s)?] 

X [ (2.320 X 10-8) /s? 

+ (5.033 K 10-*) /s* 

+ (7.649 X 10-88) /s*]. (5) 

This expression is a good approximation when s 
is of the order of 3ao or greater. 

One can find a reasonable value for A(0) by 
solving the problem anew for this special value of 
s. In this case one can no longer substitute 
1/\r,—s| for 1/|r,;—rs| in the integral over dr; 
and drs. If we assume kp to be zero, the integral is 
as follows: 


ff exp (—ar3—ik-r,)[1/|r1—rs| J 
xX (1 — br;) exp (— bri) f(rs)dridrs, 


where f(r3) is equal to unity inside the metal and 
is equal to zero outside the metal. We assume 
that half of the atom is inside the metal and the 
other half is outside, and take the polar axis per- 


% Handbuch der Physik, Vol. 24, p. 348. 
% Handbuch der Physik, Vol. 14, p. 69. 
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pendicular to the surface. We expand 1/|r;—rs| 
in terms of spherical harmonics and exp (—7k-r;) 
in terms of spherical harmonics and Bessel 


functions. The limits of integration for 63 are 
from 0 to 2/2; this takes care of the discontinuity 


of f(r;). The integration over 6, 03, ¢1, and ¢3 can 
be done immediately. Part of the integration over 
r, and rs; was done analytically, and part was 
done numerically. When the values of the con- 
stants as already defined are used, the integral 
reduces to the following expression : 


1.589 K 10-*!+-2[ 5.707 X 10-*2P; (cos 6,) 
+9.077 X 10-**P3(cos 6;) | 
+higher order terms. (6) 


Squaring this expression, integrating over all 
space, and substituting in the expression for \(0), 
one finds \(0) =1.97 X10" sec.—'. We thus have 
found a value for \(0) and an analytic expression 
for \(s) which is a good approximation when s is 
of the order of 3ao or greater. Figure 2 shows 
logio A(s) plotted against s. The heavy line was 
obtained from the expression for \(s) and there- 
fore was extended into s=3d 9. Point A corre- 
sponds to logio A(0), and point C corresponds to 
logio A(3a0); these two points are known fairly 
well. We have drawn a dotted line to show the 
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Fic. 2. Graph of logio A(s) as a function of s for He* 
and for 2*S metastable He. 


probable value of logio A(s) for s<3ao. For s = 2a9 
we took ky perpendicular to the metal surface and 
performed the integration over drs exactly and 
the integration over dr; numerically for this 
particular distance. In this integration we again 
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substitute s for rs in the denominator of Eq. (3’), 
which is not too bad an approximation for s = 2ap, 
Using our results we can get a fairly good value 
for logio A(2a9). Point B corresponds to this value 
and falls on the extrapolated curve as shown. 


d. Probability of Transition 


Actually the motion of the atom would make 
it necessary to derive a more complicated ex- 
pression involving the time-variation of the 
perturbing potential. For velocities of collision 
such that the time-Fourier analysis of the pertur- 
bation does not contain appreciable amplitudes 
for frequencies of the order of an atomic fre- 
quency, one may regard the collision as adiabatic 
and use the time-independent perturbation at 
each distance. However, it is necessary to allow 
for the possible decay of the atom at some larger 
distance. Let P(s)ds be the probability that the 
atom decays at position s in a distance ds. The 
probability that the atom reaches position s 
without decaying is 


1- [Peas 


P(s) =[N6s)/oo] 1 ~ [ Pesras as, 


where % is the velocity of the helium atom, which 
we shall assume constant. The solution of this 
integral equation is 


P(s)=B exp frr-a ‘ds[v0/d(s) ]}A(s)/vods, 


where B is a constant of integration. 

Approximating A(s) by the simple exponential 
form \(s) =A exp (—as), where A = 9.626 X 10'* 
and a=7.303 X 108 we find the expression for P(s) 
becomes 


P(s) =F exp [—(A/voa) exp (— as) —as J, 


where F is a constant. 
From the integral equation it is easy to see that 
for large s 


P(s) =X(s)/vo=[A/v0 ] exp (—as) 


and from the solution for P(s) it is obvious that 
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for large s 
P(s) =F exp (—as), 
F= A /vo 


and 
P(s) =[A/vo0 ] exp [—(A/v0) exp (—as) —as ]. 


Figure 3 shows a graph of P(s) against s for 
vo =105 cm/sec. It should be noted that the 
exponential form taken for X(s) keeps on in- 
creasing for values of s<2a, while the graph of 
logio A(s) indicates that A(s) probably levels off. 
For this reason the graph of P(s) may probably 
be less steep on the side of small s than is shown 
on Fig. 3. To get & the average distance for 
transition, one has to do the integral 


f sP(s)ds, 


which can be done graphically or analytically. 
For v9 = 105 cm/sec., § = 2a which is very near Sm, 
the value of s for which P(s) is a maximum. For 
vo = 3X10’ cm/sec. s comes out less than 0.14. 
This last value of § does not have very much 
meaning because of the breakdown of the pertur- 
bation theory and the crudeness of the model 
used for the metal. For this velocity Massey” 
obtained a value of & of the order of 3a, which 
according to our calculations is certainly too 


TABLE I. Relation between the velocity component 
perpendicular to the metal surface and the percentage of 
metastable atoms reflected. 





Vel. perpendicular 
_ to surface 
in 10% cm/sec. 


Average 
percentage 


60 to 90 75 
40 to 60 50 
40 to 70 55 
>50 

40 40 
15 to 40 28 
40 to 50 45 
10 to 30 20 
10 to 20 15 

5 


Percentage of ref. met. 
(Oliphant’s results) 





NSN ENON NEN 
coooooenoun 
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high. A very good approximation to § can be ob- 
tained by finding the s for which 


f A(s)/vods = 1. 


For vp = 105 cm/sec. one gets 1.9ap for 8. 
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From Oliphant’s” results (shown in Table !) of 
the percentage of reflected metastable atoms at 
different velocities and different angles of inci- 
dence, it can be seen that the percentage of re- 
flected metastable atoms varies fairly linearly 




















B) 
s in Bohr rodi 


Fic. 3. Graph of P(s) as a function of s for metastable 
helium atoms of velocity 10° cm/sec. 


with the component of velocity perpendicular to 
the surface. This might indicate that the per- 
centage of reflected metastable atoms depends 
mainly on how close the atoms get to the metal 
surface. 

From these data one could calculate the dis- 
tance from the metal surface at which most 
atoms turn back. Let us call this distance sg. An 
approximate value of se can be obtained by 
solving for sz in the equation 


R=1-2f A(s)/vods, 
SR 


where R is the fraction of reflected metastable 
atoms. Knowing the initial energy of the meta- 
stable atoms and knowing where they turn back, 
one can get the shape of the repulsive potential 
that the surface offers to the metastable atoms. 
Figure 4 shows this potential as a function of s. In 
this analysis we have neglected the fact that we 
should use P(s) instead of \(s) and v instead of 
vo, where v is the instantaneous velocity. When 
we found P(s) previously, we assumed that the 
atom decayed on its ingoing trip, but if we now 
consider the reflection coefficient R, we must 
solve the problem more rigorously. Let us call 
Ps) the P(s) we already calculated, meaning 
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that it is the probability that the atom decays on 
its ingoing trip. To get an expression for Po(s), 
the probability that the atom decays on its 
outgoing trip, one must solve the integral 


equation 
Pa(s)=|1 —- “Pi(sids— fPa(sdds Ws) ‘Vo. 
—_ 5R 
The solution of this equation is 
Po(s) =[A/v0] exp [—as+(A/av0) 
Xexp (—as) —(2A/avo) exp (— asp) ]. 
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Fic. 4. Potential offered by a metal surface 
to a 28S helium atom. 


! 
$s in Bohr radi: 


Now Sp can be obtained by solving the equation 


R=1 ~j Pisyas— { P(s)ds 
=exp [ —(2A/avo) exp (— asp) ]. 


Calculating sz from this formula, one finds that 
the curve representing the potential that the 
surface offers to the metastable atom is closer to 
the surface than what was previously obtained as 
shown by points A, B, C, and D on Fig. 4. We 
now correct for the fact that the velocity of the 
metastable atoms is not a constant vp but a 
function v(s) which can be obtained from v and 
the potential curve already found. By a process 
of iteration, the potential curve can be improved. 
Plotting v(s) against s, one can find an analytic 
expression that fits the curve obtained. One can 
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now find P,(s) and Po(s) by again solving the 
integral equations 


bei f "P,(s)ds )\(s)/v(s) = Pils), 





j1- “Pa(s)ds— [“Po(s)ds W\s)/0(s) = Pals). 
sR sR : 


Again sr was obtained from the equation 


R=1 - [Pass — f/Palsyas. 


This second correction brought the potential 
curve almost to coincide with the curve obtained 
by using 

R=1- 2f A(s)/vods. 
®R 
This means that apparently the two errors in- 
volved in using this approximation cancel each 
other. It is unfortunate that both Oliphant’s 
results and our theory are rather crude and that 
therefore the curve we obtained is only a rough 
approximation. It does, however, indicate a 
method to obtain information about the repulsion 
of atoms by metals. The theory can probably be 
improved considerably when there is need for it, 
so that more precise data on the reflection of 
metastables from metal surfaces are to be desired 
highly. 
e. Angular Dependence 


Equations (5) and (6) show that if the transi- 
tion takes place at the metal surface or at large 
distances from the surface, there is practically no 
angular dependence. No detailed experiments 
have been carried out to test this angular depend- 
ence, but Oliphant informed Massey™ that he 
found no electrons ejected at an angle less than 
15° with the surface. Our theory does not predict 
such angular distribution, and for this feature a 
more refined theory or experimental procedure is 
needed. For the special value s = 2d, the integral 
over dr gives 


4mn X 10-"8{5.246P (cos 6;,) +2.916P2(cos 6,) 
—0.407P,(cos 6.) +0.017P.¢(cos 4) 
+i[2.370P1(cos 6) —1.350P3(cos 6) 
+0.093P;(cos 0.) —0.003P;(cos 6) }} 

+higher order terms. 








1 





An expression whose absolute value squared 
has an angular dependence is shown in Fig. 5.” 
For thermal velocities this is the average distance 
of transition, and hence from our theory we 
should expect some angular dependence for these 
velocities. Oliphant’s velocities were of the order 
of 10’ cm/sec., and thus the atoms penetrated 
closer than 2a) before decaying, making the 
angular distribution even less than that shown 
in Fig. 5. 

Massey” explains the angular distribution 
found by Oliphant by assuming that “‘the portion 
of the wave function which gives rise to the 
ejection of electrons is that which falls within the 
boundary of the metal.’’ Hence when a transition 
occurs, he regards the atomic electron as over- 
lapping the metal, and when this electron is 
replaced by the metallic electron, it must there- 
fore pass through the metal surface. Massey 
further states that the electron will be refracted 
similarly to a ray of light going from one medium 
to an optically denser one. There seems to be no 
basis for the assumption that the atomic electron 
must always be inside the metal for the transition 
to take place, as an inspection of Eq. (1) clearly 
shows. Again, we would expect that when the 
electron is refracted as it passes out through 
the metal surface, it is refracted away from the 
normal because its initial velocity parallel to 
the surface is unaltered after going through the 
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Fic. 5. Graph of P(8, 6) as a function of 6 for 2°S 
helium atoms of velocity 10° cm/sec. 


26 For this computation and for most of the other com- 
gers the WPA tables of exponential function and the 


VPA table of spherical Bessel functions were used. 
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surface but its normal velocity is diminished 
because of the work function. Hence if we were 
to include this effect in our theory, it would 
decrease the preference for electron emission 
perpendicular to the metal surface. 


Ill. HELIUM ION CASE 
a. Energy Balance 


The case of the helium ion colliding with a 
metal surface can be solved in a way analogous to 
the case of the helium metastable atom. Figure 6 
shows the energy balance diagram for a Mo 
surface. This diagram shows that inside the 
metal there are electrons whose energy level is 
identical or very nearly identical to the 2°S orbit. 
When the ion gets near the metal surface the 
potential field to which the metallic electrons are 


Zero Energy 


2's 
23s 








eo 1s 


Fic. 6. Energy balance diagram for the extraction of 
electrons by helium ions from molybdenum. This diagram 
should be interpreted as explained under Fig. 1. 


exposed is altered, and cne of the above-men- 
tioned electrons may pass through the metal 
surface and occupy the 2°S orbit. In the case of 
Mo there are no electrons whose energy level is 
identical to the 24S or 1S, and for this reason the 
only probable transition is the one to the 2°S 
state. In the case of Pt, the work function is 5.01, 
5.1, or 6.0 v?? as obtained by different investi- 
gators, and our diagram would show that then 
the transition to the 2*S would be impossible. 
This is not exactly correct because as the temper- 
ature of the surface increases, some of the higher 
electronic cells will be occupied. As shown by 
Sonkin‘ the behavior of a metallic surface is very 
easily affected by the conditions of the surface. 


27 Handbuch der Physik, Vol. 14, p. 69. 
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b. Matrix Element 


In this case the matrix element is 


H= f f f [u3.(r1)130(r2) —use(r2)uta(r1)] 


X [—Ze?/re+e?/ris |uis(r1) 
X Lexp io: r2]/+/Qdridr2 


where the variable r. must remain inside the 
metal. A good approximation for H can be ob- 
tained for large s by replacing r;2 by re. If we 
make this assumption, H reduces to 


H= { {1-0 exp (—bre—cr;)[1/re } 


Xexp [—ari+iko(cos 6k cos 62 
+sin 0% sin 02 cos $2) \dridr. 


The integration over dr; can be done immedi- 
ately. The integration over dr, can be done 
without much difficulty if we expand the Bessel 
function obtained from the integration over ¢2 
and take only the first term of the expansion. 
The final result is 


H=16r?/(a+c)*[exp (iko cos Oos — bs) ]/bs 
X [(bs?—s)/iko cos Oeo—b 
+(1—2bs)/(iko cos Oko —b)? 
+2b/(iko cos Oko—b)? }. 


The value of ||? can be found for any large dis- 
tance s and any direction of ko. For s of the order 
of 10a) we found that the value when kp was 
parallel to the metal surface was about four 
times the value when ko was perpendicular to the 
metal surface. We took an average of these two 
values. 


c. Decay Rate 


The rate at which ions decay as a function of 
the distance from the metal is now given by 
A(s) = [mk Q/xh* || H|?. This expression unlike the 
corresponding one for the case of the metastable 
atom does not have a factor nQ, the number of 
free electrons. This extra factor arose because in 
the averaging over kp all free metallic electrons 
were available. Any excess energy could be 
carried away as kinetic energy by the ejected L 
shell electron. In the case of the ion there is no 
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final outgoing electron, so the energy of the 
metallic electron capable of performing the trans- 
ition is uniquely defined by energy conservation 
within the limits of the uncertainty principle. 
When the numerical values of the constants are 
substituted, 


A(s) = 2.390 K 10” exp (— 2.254 x 108s) 
X [s/b+1/b?+1/b's }2, 


Figure 2 shows a graph of logy» A(s) as a function 
of s. 
One can find § the average distance of transi- 


tion from the integral 


§ = [ sP(s)as 


As previously shown, a good approximation is 
obtained by finding the s at which 


{ A(s)/vods = 1. 


For vp = 105 cm/sec., = 11.5a9 and, for v9 = 3 X 107 
cm/sec., §=6.0dp. 

Massey solved the case of a hydrogen ion which 
collides with a metal surface and performs the 
transition to the 28S state. He found that for 
vo = 3X10? cm/sec. 8+ is of the order of lao, and 
then assumed that sq.+ might be bigger by a 
factor of two or more. We calculated 84+ for 
vo=3X10’ cm/sec. and obtained an § of the 
order of 5a9. For the case of the hydrogen ion 
Massey used the correct matrix element but did 
not evaluate the integrals correctly. 


IV. DISCUSSION 


The extraction of electrons by He? colliding 
with a metal surface may be explained as the 
result of two separate processes. First the ion is 
neutralized, becoming a 2°S metastable atom, 
and then if the metastable atom can get closer to 
the metal surface than Smet, it will have a very 
high probability of assuming the ground state, 
the resulting energy being used in extracting an 
electron from the metal. Figure 7 shows that our 
values Sye+ and Smet for thermal velocities permit 
such a double process. For 1» =3X10' cm/sec. 
Smet Might seem too small to permit ions to be 
reflected as neutrals, but one must remember that 
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Fic. 7. for 28S helium atoms and for He* for 
two different velocities. 


this is an average distance, which means that 
some ions after being changed into metastables 
will decay at a distance greater than 0.1dp. 
Besides, at these very small distances our theory 
can give only very approximate results. 
Oliphant and Moon”! in their discussion of the 
helium ion proposed the theory that the inverse 
transition metastable—ion was highly improbable 
because once the transition ion—metastable 
takes place the remaining metallic electrons 
would fill the cell previously occupied by the 
electron which is extracted from the metal. This 
would be true if the relaxation time were of the 
order of 10-5 sec. or less. A rigorous calculation 
of the relaxation time®® is too involved, and we 
shall not attempt it. Nevertheless, our value of 
Sue+ (6.0a9 for v9 = 3X10’ cm/sec.) indicates that 


28 From dimensional considerations one may guess the 
following expression for r. 
7r=A/(vn'(e*n!/mv*)*], 
where v is the velocity of the metallic electrons, n their 
density, m the electronic mass, and e the electronic charge. 


A is an undetermined constant. The numerical computa- 
tion shows that r~A10-" sec. 


a beam of Het ions colliding with a metal surface 
will be almost completely transformed into 
metastable and neutral atoms unless the meta- 
stable—ion transition is also highly probable. 
When Oliphant” produced a beam of metastable 
atoms by sending a beam of helium ions at a 
small angle with the metal surface, he also ob- 
tained many reflected ions which he removed 
from the beam by means of an electric field. This 
experimental evidence indicates that the relaxa- 
tion time is too long to prevent the metastable 
—ion transition to take place. 

We will therefore assume that the metallic 
electron will oscillate between the metal and the 
ion (Schiittelwirkung) once the ion gets closer to 
the metal surface than 8y.+. Thus if a beam of 
ions is sent against a metal surface, once they are 
closer to the metal surface than the 8.+, they 
can no longer be identified as ions or metastable 
atoms.”® The probability that an ion finally may 
achieve the ground state depends on how close 
the particle can get to the metal surface and on 
the number of ion@metastable jumps that can 
occur. These two factors depend on the energy of 
the ion, but the problem of finding this functional 
dependence is too involved. 
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29 Most metals have a work function less than 5.0 volts 
and for this reason if 23S helium atoms are sent against 
these metals at absolute zero, all the electronic cells are 
filled up and the metastable atoms would remain as 
metastables until they perform the transition to the 
ground state and extract an electron. 
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In the presence of foreign gases, a narrow and structureless absorption maximum was 
observed on the short wave-length side of certain members of the principal series of alkali 
metals. A résumé of the results obtained so far is given. It is suggested that the band is due 
to the increase in the energy of the optical electron through the process of collision with a 
foreign gas atom when the electron is at that energy level whose orbital radius is close to the 


distance of optical collision. 





INTRODUCTION 


HE influence of surrounding atoms on the 
frequency of the absorbing or emitting 
atoms as exhibited by broadening, asymmetry, 
and shift is well known and gives a new way to 
study the perturbations of energy levels; but 
during the study of the intensity distribution of 
the line due to pressure broadening, the narrow 
diffuse band':? appearing very close to certain 
members of the principal series is very easily 
overlooked, yet it appears to be of considerable 
importance. Ny and the author first observed the 
band on the short wave-length side of the second 
member of the rubidium principal series in ab- 
sorption.! Later on the corresponding band was 
also observed for caesium.’ It was found that the 
width of the band and its position with respect 
to the corresponding series member depend on 
the nature of both the colliding and the absorbing 
atoms. The lighter the colliding or the absorbing 
atoms the greater the width and the separation 
between the series line and the band. In the in- 
vestigation of the pressure effects of foreign 
gases on spectral lines the concentration of the 
foreign gas is always sufficiently high for the 
band to appear, even when the temperature of 
the absorption tube is rather low. The corre- 
sponding phenomenon was also studied for so- 
dium.* It is very interesting to note that, for 
sodium, this sort of absorption band appears 
only in the neighborhood of the third member of 
the principal series. Very recently the phenom- 
enon was also observed by the author and his 
* In memory of the author's mother. 
1Ny Tsi-Zé and Ch’en Shang-Yi, Nature 138, 1055 
OD Ny Tsi-Zé and Ch’en Shang-Yi, J. de phys. et rad. 9, 
169 (1938). 


* Ke T’ing-Sui and Ch’en Shang-Yi, Chinese J. Phys. 5, 
publication delayed by the war. 


collaborators for potassium. The narrow and 
structureless absorption maximum was observed 
on the short wave-length side of the second 
member of potassium principal series, and no 
bands of this kind were observed in the vicinity 
of the other members of the series. It is the 
purpose of this paper to present the results 
newly obtained for sodium and potassium and 
to give a résumé of the results obtained so far 
with a theory to account for the phenomenon. 


EXPERIMENT 


The experimental arrangement was essentially 
the same as that previously described.‘ Nitrogen 
and hydrogen were used as foreign gases in the 
presence of sodium and potassium. Pure nitrogen 
was prepared by decomposing ammonium nitrite. 
Hydrogen was prepared by the electrolysis of 
water whose contained air was first removed by 
repeated boiling and cooling in “vacuum.” 
Several stages were also used to remove various 
kinds of possible impurities. Before the admission 
of hydrogen or nitrogen to the absorption tube, 
the gas was allowed to pass through a vapor trap 
of the same alkali as that whose absorption 
spectrum was to be studied, so that the residue 
of impurities (especially water vapor and oxygen) 
present in the gas were removed. A high intensity 
tungsten filament lamp and a hydrogen discharge 
tube were employed as light sources for the 
visible and the ultraviolet parts of the spectrum. 
A Hilger E-316 and a 10-foot grating spectro- 
graph were used to register the entire absorption 
series of the alkali. The appearance of the band 
was studied by means of an automatic recording 
microphotometer. 

Spectra were taken when the absorption tube 


‘Ny Tsi-Zé and Ch’en Shang-Yi, Phys. Rev. 51, 567 
(1937). 
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was heated to various temperatures. Great care 
was taken, so that the possibility that weak lines 
were overlooked owing merely to improper ex- 
posures is very small. 


RESULTS AND DISCUSSION 


Table I gives a résumé of the data obtained so 
far on the absorption bands of alkalis in the 
presence of foreign gases. It has been found that 
the pressure of the foreign gas does not appre- 
ciably affect the appearance of the band as long 
as the pressure is much higher than that of the 
alkali vapor. The pressure of the foreign gas is 
not so important and, therefore, is not listed in 
the table. In fact, for potassium the pressure in 
the absorption tube was 34 and 45.5 cm Hg for 
nitrogen and hydrogen, respectively. The corre- 

TABLE I. New absorption bands of alkali metals in 


the presence of foreign gases. 








ture of 
absorp- 
tion tube 


A °C 
41.7 (3rd member) 544 
4007 .2 38.5 (2nd member) 450 
Rb 4178.7 7 23.1 (2nd member) 300 
Cs 4541 : 14 (2nd member) 280 


Rb 4182.5 19.3 (2nd member) 220 
Cs 4543 12 (2nd member) 200 


Wave- Width 
length of 
of band band 


Separation between 
the band and the 
series member 


Foreign 
gas Alkali 





A 
Na 2811.2 
K 


34.8 (3rd member) §22 
27.0 (2nd member) 484 
18.8 (2nd member) 240 
(2nd member) 210 


Na 2818 
K 4018.7 
Rb 4183 
Cs 4544 11 


7.3 (2nd member) 220 
(2nd member) 


Rb 4194.5 
Cs 4549 6 


no band 
Cs no band 











sponding values for sodium were 48 and 30 cm 
Hg, and for all the rest of the data the pressures 
were about one atmosphere. 

The last column of Table I gives the most 
favorable temperatures of the absorption tube 
for observation of these bands. Certainly the 
lowest temperatures to render the band per- 
ceptible are much lower than those listed. For 


potassium, the lowest temperatures at which the 


band became visible were 410°C for hydrogen 
and 380°C for nitrogen. In all cases the tempera- 
tures were rather low, so that no molecular bands 
of alkalis or of alkali-foreign-gas compounds 
appeared. 

The readings in the fifth column of Table I 
show in each case the distance between the 
maximum of the band and the *P3;2 component 
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of the member of the series. The readings for 
the width (i.e., the half-width) of the band are 
only approximate. 

Oldenberg, Kuhn, and others® have studied 
some bands which may conceivably be a similar 
phenomenon in emission spectra. They used 
various other metals and observed one or more 
bands appearing near all spectral lines. Obviously 
owing to the many complicating factors, such as 
the non-uniformity of the temperature of the 
discharge and of the vapor pressure of the metal, 
the pressure of ions, dipoles, etc., the condition 
for the production of these bands in emission is 
not so simple as that in absorption. The bands 
were interpreted by Oldenberg, Kuhn, and others 
as due to loosely bound molecules resulting from 
close collisions. According to their theory, bands 
of this type would occur with all metal lines, in 
agreement with their results. 

Their theory is evidently not confirmed in the 
present work, in that this type of band appears 
only in the neighborhood of certain spectral lines. 
Nor is there any obvious rule for selecting these 
lines, as illustrated by the fact that in the series 
of sodium the absorption maximum appears by 
neither of the first and second members nor by its 
weaker members, except the third, 3*°S1j2—5*°P 3/2. 
Again for K, Rb, and Cs, the band appears only 
beside the second member of the principal series. 
Another fact also unexpected is that argon pro- 
duces no band for Rb and Cs. Whether the bands 
observed by Oldenberg and Kuhn are of the 
same nature as those observed here is quite 
doubtful. 

In the light of a systematic observation of this 
phenomenon as shown in Table I, it is easy to 
single out five points: 

(1) The lighter the absorbing or perturbing 
atoms, the greater the separation between the 
band and the series line. Hence for a certain 
foreign gas the separation corresponding to Cs 
has the lowest value while that corresponding to 
Na possesses the highest value. For a certain 
absorbing alkali the separation corresponding to 
the lightest gas has the highest value and that 
corresponding to the heaviest gas has the lowest 
value. Consequently the separation correspond- 

5 Oldenberg, Zeits. f. Physik 47, 184 (1928); 55, 1 (1929). 
Kuhn and Oldenberg, Phys. Rev. 41, 72 (1932). Krefft and 


Rompe, Zeits. f. Physik 73, 681 (1932). Preston, Phys. 
Rev. 51, 298 (1937). 
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ing to H2—Na has the highest value of all and 
that corresponding to Ne—Cs has the lowest 
value of all; as to the extreme case of argon and 
Rb or Cs, no band can be resolved. 

(2) The lighter the absorbing or perturbing 
atoms, the greater the width of the band. Ex- 
amples given in (1) can be similarly applied here. 

(3) The wider the band, the higher the tem- 
perature of the absorption tube or the concentra- 
tion of absorbing atoms that is needed to render 
the band perceptible. 

(4) For heavier alkalis, viz., K, Rb, and Cs, 
the band appears on the short wave-length side 
of the second member of the principal series, 
while for a lighter alkali, Na, the band appears 
on the short wave-length side of the third 
member of the series. It can be predicted without 
hesitation that the corresponding band for 
lithium will appear on the short wave-length 
side of the third or still higher member of the 
principal series. 

(5) The band observed is narrow, structure- 
less, and much fainter than the intensity of the 
series lines. 

One more point should also be mentioned here.*® 

(6) When the pressure of the perturbing gas is 
very large in comparison with that of the ab- 
sorbing vapor, the intensity of the band depends 
only upon the concentration of the absorbing 
vapor. If the pressure of the perturbing gas is 
reduced so that the concentration of the absorb- 
ing vapor becomes comparable with that of the 
perturbing gas, the intensity of the band will 
decrease with the decrease of concentration of 
either the absorbing or the perturbing atoms. 

The production of the band may be suggested 
as due to the increase in the energy of the optical 
electron through the process of collision with the 
foreign gas atoms when the electron reaches that 
energy level whose orbital radius is close to the 
distance of optical collision. Foreign gas atoms 
can be regarded as evenly distributed in the 
surroundings of the absorbing atom. It is the 
effect, say the polarization effect, of the introduc- 
tion of the foreign atoms within the space of the 
possible orbits of the optical electron of the 
alkali that contributes certain features of the 
well-known “pressure effects.”” As the spacing 
between foreign atoms is quite large, and the 


6 Cf. references 1 and 2. 


foreign atoms are moving with an enormous 
speed, the chance of collision between the excited 
electron and the foreign gas atom is rather small. 
But the chance of collision becomes considerable 
when the foreign gas atom reaches the location 
where collision between the foreign gas atom and 
the absorbing atom takes place because the 
colliding atom will spend the major part of its 
time around the district of closest collision. If 
that distance of closest collision happens to be 
near a certain energy level of the alkali atom, 
only the electron at that state will have appreci- 
able chance to change energy. With this picture 
in mind the above six points can be accounted for 
without difficulty. 

The magnitude of the separation between the 
band and the series line is a measure of the energy 
change of the optical electron at the instant of 
collision with the foreign gas atom. The lighter 
the colliding particle, the higher the relative 
speed of collision; consequently more energy 
will: be given to the electron when it is excited 
to the energy level whose orbital radius is close 
to the distance of collision. For heavier alkalis, 
Rb and Cs, and a heavier gas, say argon, the 
separation might be so small that no band is 
perceptible. This is illustrated in point (1). The 
distribution of velocities has a wider range for 
lighter absorbing or perturbing particles, so the 
structureless band is wider for lighter atoms. 
This explains point (2). Point (3) is simply a 
matter of photometry because the amount of 
absorption depends only upon the concentration 
of the absorbing atom. For a broader band the 
energy is distributed over a wider range, thus a 
higher concentration of the alkali vapor is re- 
quired to obtain the perceptible intensity. 

Point (4) tells the approximate distance of 
optical collision. By applying the quantum 
mechanics expression for the average value of 
the electron nuclear distances at various energy 


levels, it is possible to tell the approximate value 


of the distance of collision between the nuclear 
and the colliding atom. The value for the third 
member of Na principal series, 3s — 5p, is 19A. For 
the second members of K(4s—5p), Rb(5s— 6), 
and Cs(6s—7p) series the values are 19A, 28A, 
and 38A, respectively. Points (5) and (6) are 
simply a reflection of the probability of this kind 
of transition. 
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Cosmic-Ray Meson Spectra 


D. M. Bose, BinHA CHOUDHURI, AND M. SINHA* 
Bose Institute, Calcutta, India 


(Received April 22, 1944) 


N a recent paper' (denoted hereafter as 1), 
Hamilton, Heitler, and Peng have proposed a 
comprehensive theory of the origin of cosmic-ray 
mesons, which appears to provide at least a 
qualitative interpretation of a large number of 
observations on cosmic-ray penetrating particles. 
These authors aim at giving a connected account 
of cosmic-ray and nuclear phenomena, and have 
selected that form of the meson theory which 
gives the best account of nuclear forces, viz., that 
proposed by Mller and Rosenfeld.? According to 
Heitler et al.,* the meson field of a fast proton of 
energy E>>uc? when scattered in the field of a 
stationary nucleon, radiates an amount of meson 
energy e~£E. We take the meson mass energy uc’ 
to be 10* ev. The released energy ¢ may (1) result 
in a single meson of energy e which in the process 
of being scattered in the field of the stationary 
nucleon gives rise to a single or a multiple of n 
mesons. The cross section y, for each multiplicity 
n has a maximum at some energy ¢, which in- 
creases with nm. After the maximum, y, decreases 
rapidly, and the asymptotic expressions are 


yi~we?; ya~we' for n#¥1. (1) 


Alternately (2) it may also happen that only a 
part e’ of the radiated energy is used up in meson 
production, and the balance is taken up by the 
recoiling nucleon; the latter can then produce 
fresh mesons by nuclear collisions. This process 
is analogous to cascade production of electrons. 

On the assumption of a suitable primary 
proton energy spectrum incident on the top of 
the atmosphere, an expression is deduced for the 
energy distribution in the single-meson spectra 
f(e, x), where x is the depth from the top of the 
atmosphere. It is shown that the function f(e, x) 
for a given value of x has a maximum value for 
e~3yuc* below which it falls rapidly to zero. These 


critical values of e~1/f where f=f’~} are the 


* At present at the Indian Institute of Science, Banga- 
lore, India. 

1 Hamilton, Heitler, and Peng, Phys. Rev. 64, 78 (1943). 

2 Mdller and Rosenfeld, Kgl. Danske Vid. Sels. Math.- 
Fys. Medd. 17 (1940). 

* Heitler and Peng, Proc. Camb. Phil. Soc, 38, 296 (1942). 
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coupling constants between a nucleon and the 
vector and pseudoscalar components of its meson 
field. The decay constants of the corresponding 
mesons are taken to be 10-* and 10-® sec., 
respectively. 

Given below are some of the results we obtained 
which are in agreement with the predictions of 
the theory. 


MESON SPECTRA 


Several investigators have used Wilson cham- 
ber photographs of penetrating particles in mag- 
netic fields to measure the meson energy spectra. 
The most reliable measurements for the low 
energy spectra appear to be those of D. J. 
Hughes.‘ His energy distribution curve shows a 
maximum near e~8 X 10° ev, and the number of 
penetrating particles of lower energy diminish 
rapidly. This is not in agreement with the 
distribution f(e, x) in which the maximum occurs 
for e~3X10* ev. 

Photographic plates exposed to cosmic rays at 
high altitudes are found to contain single and 
multiple star-like tracks. A detailed account of 
the results will be published by one of us (B.C.) 
elsewhere. In Table I are given our findings on 
the total frequency of multiple track distribution 
observed in photographic plates, kept at altitudes 
varying from 7000 ft. to 14,500 ft., under (1) air, 
(2) absorbers like water, paraffin, and wood, and 


TABLE I. Frequency of multiple tracks. 








Absorber 


n= 1 2 3 4 5 
Air 192 30 87 63 38 
Hydrogeneous 
substances 86 14 
Lead 
0.5cmto5cm 69 


Frequency of tracks 





22 20 20 0 O 


19.6 42.3 34.2 12.9 2.2 0 








(3) under lead. In every case the frequency distri- 
bution has a principal maximum at m=1, a 
second maximum at »=3; and in no plate was 
any multiple track with »>7 observed. The first 


‘D. J. Hughes, Phys. Rev. 57, 592 (1940). 
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maximum is not significant, as it contains mainly 
tracks of ionizing particles originating outside the 
photographic emulsion. In a paper® on meson 
spectra, reasons were given for the conclusions 
that (1) the star-like tracks on photographic 
plates are mainly due to multiple meson produc- 
tion, the energy associated with each multiple n 
is €,=n-yuce?+k, where k~10’ ev represents the 
kinetic energy of the particles; (2) protons 
(neutrons) are primarily responsible for such 
multiple production, (3) a transition effect is ob- 
served in lead, with a maximum under 1.5 cm 
Pb; in this case photons directly or indirectly are 
responsible for the additional multiples produced. 
The occurrence of a maximum at n=3 was 
noticed, but no satisfactory explanation for it 
could be suggested. 

According to I, the frequency distribution 
curve for the occurrence of multiples of 2 mesons 
is approximately given by 


fle, x) + (¥n/¥1)- (2) 


The asymptotic expression for y,/y1 is «~*~, 
which shows that within the energy interval for 
which the values of y have been calculated, the 
ratio diminishes exponentially, i.e., the expecta- 
tion of occurrence of high multiplicity becomes 
vanishingly small with increasing m. According to 
our findings (2) has a maximum for n=3 corre- 
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Fic. 1. Track in photographic emulsion traced by use of a 
camera lucida. Magnification 760. 


5D. M. Bose, Trans. Bose Research Inst. 15, 55 (1942). 
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Fic. 2. Wilson cloud-chamber photograph showing 
production of mesons by a cascade process. 


sponding to e~3yuc?, and this may be taken as a 
qualitative verification of the occurrence of the 
maximum predicted according to |. Further it 
appears to us to be more than a mere coincidence 
that the region of maximum occurrence of single 
mesons, e~8 X 10° ev, should also coincide with 
the region where the production of multiple 
mesons cease (i.e., > 7yuc*). 


CASCADE PRODUCTION OF MESONS 





According to I a cascade production of 
mesons can also be expected. Hazen® finds no 
evidence of their occurrence in his high altitude 
Wilson chamber photographs. On the photo- 
graphic plates referred to in the previous para- 
graph, records of at least six such events have 
been found. A camera lucida tracing of one of 
them is given in Fig. 1. Since photographic plates 
are exposed to cosmic rays for at least six months, 
it is difficult to adduce conclusive evidence that 
all the tracks on a plate which are considered to 
be due to cascade meson production have been 
traversed by the mesons at the same instant of 
time. A Wilson chamber photograph obtained by 
one of us (M.S.) is given in Fig. 2. It shows the 






®°W. E. Hazen, Phys. Rev. 64, 257 (1943). 
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LETTER 


production of mesons by a cascade process. All 
the tracks appear to be of the same age. The 
tracks in Group 1 have originated outside the 
chamber; the other two groups have originated 
within the chamber. The left-hand track of group 
2 is absorbed in the lead plate (2.2 cm), and its 
energy is <4X10’ ev. The four tracks in group 3 
appear to start as a pair which is doubled again. 
These tracks are not appreciably scattered during 
their traversal through the lead plate, but they 
have produced knock-on electrons. 


BURST PRODUCTION 


One of us (M.S.)’ has recently published a 
photograph of a burst originating in the gaseous 


7M. Sinha, Phys. Rev. 64, 248 (1943). 
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volume of a Wilson chamber. It appears to start 
as a number of ionization strands, diverging from 
a small region, which subsequently broaden and 
merge into a dense volume of ionization. A 
possible interpretation of the burst according to I 
is as follows: A high energy proton has by 
process of nuclear collision given rise to a cascade 
production of fast vector mesons. The latter have 
within a short length disintegrated into high 
energy electrons, from which several cascade 
showers have started. 

The above theory enables us to correlate and 
interpret many findings on penetrating particles 
obtained from Wilson chamber photographs and 
from tracks on photographic plates. A detailed 
comparison between the predictions of the theory 
and our findings will be given elsewhere. 
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gyre publication of brief reports of important dis- 
coveries in physics may be secured by addressing them 
to this department. The closing date for this department is the 
third of the month. Because of the late closing date for the sec- 
tion no proof can be shown to authors. The Board of Editors 
does not hold itself responsible for the opinions expressed by 
the correspondents. Communications should not in general ex- 
ceed 600 words in length. 


On the Maximal Energy Attainable’ 
in a Betatron 


D. IWANENKO AND I, POMERANCHUK 


Physical Institute of the Moscow State University, Moscow, and 
Physico- Technical Institute of the a of Sciences of the 
. S. S. R., Leningrad, U. 


May 18, 1944 
Y means of a recently constructed induction accelera- 
tor-betatron, Kerst succeeded in obtaining electrons 
up to 20 Mev.' The principle of operation of the betatron 
is the acceleration of electrons by a tangential electric field 
produced by a changing magnetic flux, which is connected 
with the magnetic field keeping electrons on the orbit by a 
simple relation. In contrast to a cyclotron, whose applicabil- 
ity is essentially limited to the non-relativistic region on the 
ground of defocusing of orbits due to the change of mass at 
high energies, there is no such limitation for the betatron. 
We may point out, however, that quite another circum- 
stance would lead as well to the existence of a limitation 
for maximal energy attainable in a betatron. This is the 
radiation of electrons in the magnetic field. Indeed, elec- 
trons moving in a magnetic field will be accelerated and 
must radiate in accordance with ghe classical electro- 
dynamics. One can easily see that quantum effects do not 
play here any important role as the dimension of the orbit 
is very great. As was shown by one of us? an electron moving 
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in a magnetic field H radiates per unit of path the energy 

— (dE/dX) =2/3(e2/mc?)?(E/mc?)*((V/oH ? (1) 

where e is the charge, m the mass, V the velocity, and E the 

energy of the electron; E is assumed much greater than mc’. 

In the betatron V is normal to H and practically for the 
whole path equal to c. Then we have 

— (dE/dX) =2/3(e*/mc*)*(EH/ mc*)?*. (2) 

The limiting value of energy Eo is to be determined from 

the condition that the radiated energy (2) will be equal to 

energy gained by the electron in the electric field produced 

by magnetic flux per unit of path:* 








2, of EoH dq/ dt 
off) = Sido eR, \ | (3) 

mc 2rRoc 

H =dH/dt ro=e?/mc*. 
Here Ro is the radius of the orbit, ¢ is the induction flux.! 
Hence: 
Eo ~ (5k fy ; (4) 
me? \2retc H? 


Taking for H and E the values now being in use we get 

«5X 10* ev, which is only five times as great as the 
energy which one expects to obtain in the betatron now 
under construction. From (4) one sees that Ep is inversely 
proportional to the magnetic field applied and propor- 
tional to the square root of energy gained in the rotation 
electric field per unit of path. All this requires the using of 
smaller H or of higher frequencies with the purpose of 
getting higher limiting values of Eo. 

The radiative dissipation of energy of electrons moving 
in a magnetic field must be also of importance for the dis- 
cussion of the focusing of the electronic beam, as the energy 
of particles being accelerated will grow more slowly with 
the growth of H if the radiation is taken into account. 
This latter question may deserve a separate discussion. 


'D. W. Kerst, Phys. Rev. 61, 93 (1942). 
21. Pomeranchuk, J. Phys. 2, 65 (1940). 
3D. W. Kerst and R. Serber, Phys. Rev. 60, 53 (1941). 
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HE 260th meeting of the American Physical 

Society was held in the building of the 
Mellon Institute at Pittsburgh, Pennsylvania, on 
Friday and Saturday, the 28th and 29th of April, 
1944. This proved to be one of our strongest 
meetings since Pearl Harbor; for instance, there 
were forty-four contributed papers, though the 
authors of some of these were not able to be 
present. A good share of the credit belongs to the 
officers of our new Division of Electron and lon 
Optics, who assembled a programme of seven 
invited papers and induced the presentation of 
several of the contributed papers. This was the 
debut of the Division, and augurs well for the 
prosperity of its future. The titles of the invited 
papers are as follows. 















Address of Welcome by the Chairman or the Vice Chair- 
man of the Division. 

Historical Background of Electron Optics. C. J. CALBick, 
Bell Telephone Laboratories. 

The Electron Microscope in Metallurgical Research. 
C. S. BARRETT, Carnegie Institute of Technology. 

Extending Microscopic Examination of Metals. F. 
KELLER AND A. H. GEISLER, Aluminum Company of 
America. 

The Magnetically Focused Radial Beam Vacuum Tube. 
A. M. SKELLETT, Bell Telephone Laboratories. 

Kinematics of Reflection Oscillators. A. E. HARRISON, 
Sperry Gyroscope Company. 

The Thickness of Electron Microscopic Objects. RoBLEy 
C. WILLIAMS AND RALPH W. G. Wyckorr, University of 
Michigan. 












Another good share of the credit belongs to the 
Physical Society of Pittsburgh, whose officers 
arranged for four invited addresses as follows. 






Address of Welcome. SIGMUND HAMMER, Gulf Research 
and Development Company. 

The Mechanism of Crystal Growth and Its Conse- 
quences. WHEELER P. Davey, The Pennsylvania State 
College. 

Modern Theories of Plasticity in Solids. FREDERICK 
Seitz, Carnegie Institute of Technology. 

Hipersil, the New Transformer Iron. T. D. YENSEN, 
Westinghouse Research Laboratories. 
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Another feature of the meeting was a sym- 
posium on cosmic rays with five participants as 
follows. 


The Nature of Cosmic-Ray Phenomena. Marce. 
SCHEIN, University of Chicago. 

The Role of the Primary Cosmic Radiation. W. F. G. 
Swann, Bartol Research Foundation. 

On the Lateral Spread of Cosmic-Ray Showers. L. W. 
NORDHEIM, Duke University. 

Nuclear Disintegrations Produced by Cosmic Radiation. 
M. GOLDHABER, University of Illinois. 


Production of Neutrons by Cosmic Radiation. S. A. 
KorrFF, New York University. 


Five additional papers on a diversity of sub- 
jects, listed below, were also presented. 

Potential Barriers to Internal Rotation in Organic Com- 
pounds. J. G. Aston, The Pennsylvania State College. 


Photoelectric Measurement of Weak Emission Spectra. 
D. H. Rank, The Pennsylvania State College. 

The Philosophical Foundations of Quantum Mechanics. 
HANS REICHENBACH, University of California. 

Physics and Phenomenology. HENRY MARGENAU, Yale 
University. 

Physics Bulletins for High School Bulletin Boards (The 
Pennsylvania Conference of College Physics Teachers). 

No local committee has ever handled one of 
our meetings better than the one headed by Mr. 
A. G. Worthing and comprising several others of 
the physicists of Pittsburgh. About four hundred 
people registered for the meetings, and over two 
hundred were present at the dinner, which was 
held at the Hotel Schenley on Friday evening 
with President A. J. Dempster presiding over the 
after-dinner programme. 

At the meetings of the Council, one hundred 
and twenty-six candidates were elected to mem- 
bership and four Members were advanced to 
Fellowship; the list is appended. 

The Council approved a text of by-laws for the 
Division of Electron and lon Optics, excerpts 
from which are being distributed to the members 
of the Division. The Council also approved a text 
of by-laws for the Division of High Polymer 
Physics, contingent on certain modifications 
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which are now being made and are being sub- 
mitted to the Organizing Committee for this 
Division. A petition for a Division of Applied 
Spectroscopy having been received, the Council 
acted upon it by passing a motion of which a part 
is here quoted: ‘The Council approves in princi- 
pal the formation of a Division of Spectroscopy: 
it authorizes the President and Secretary to 
confer and correspond with signers of the petition 
of Applied Spectroscopy concerning the scope 
and name of the proposed division.” 

The Society has lost through death the follow- 
ing Members and Fellow: F. L. Cooper (Yale 
University), W. S. Day (formerly of Columbia 
University), and John L. Rose (F., New York 
University). 


Elected to Fellowship: Allen, J. S.; Burt, R. C.; Salisbury, 
W. W.; Staub, Hans. 

Elected to Membership: Alzofon, Frederick E.; Anway, H. 
Wilbur; Bader, Frank; Beer, Albert Carl; Berenda, Carlton 
W.; Berry, Chester Ridlon; Beth, Erich W.; Birzis, Lucy; 
Blum, Victor J.,S.J.; Bobrowsky, Alfred R.; Boor, Ladislav; 
Bowman, John Charles; Braendle, H. A.; Byerly, Paul R., 
Jr.; Burbank, Burr G.; Carlson, F. D.; Cocke, Albert 
Kirven; Cohen, Robert Sonne; Conner, Willard P., Jr.; 
Cooper, Marcella; Cosler, Arthur Stackdale, Jr.; Cox, 
David C.; Davis, Francis J.; Deutsch, Robert A.; Dietz, 
Frank Tobias; Di Toro, Michael J.; Dwight, Barbara; 
Eastmond, E, John; Eggen, Donald T.; Eisner, Elmer; 
Elder, Fred Kingsley, Jr.; Feitelberg, Sergei; Ferrell, O. 
Perry; Frazier, Thomas Vernon; Gans, David M.; Gettner, 
Mark Earl; Gibson, Robert John, Jr.; Gross, Milton H.; 
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Hall, James W.; Halpern, Isaac; Hart, Stephen D.; Havas, 
Peter; Hemenway, Curtis L.; Hickman, R. E.; Holton, 
Gerald James; Jeffries, Jasper Brown; Jones, V. L.; Jordan; 
A. Raymond; Ke, T’ing-sui; Kleckner, U. Frederick, 
Knight, Walter David, Jr.; Kuhlthau, Alden Robert; 
Laing, Kenneth M.; Landers, Harry M., Jr.; Lapsley, 
Alwyn C.; Lee, Donald E.; Lefler, Ralph W.; Lepore, 
Joseph V.; Lock, Charles; Loeb, Arthur L.; Lofgren E. Bel 
Lovins Gerald H.; Loughman, James C.; MacLean, 
William R.; MacRae, Duncan, Jr.; Mann, Alfred K.; 
Marcus, Jules Alexander; McDonald, Robert S.; Merritt, 
Thomas P.; Mills, Peter J.; Moeckel, Wolfgang E.; Moore, 
Donald Clark; Moore, John Wilson; Needels, Theodore S.; 
Nicoll, Frederick H.; Oetjen, Robert A.; Ohlheiser, Carlton 
E.; Osborne, M. F. M.; Owens, Ralph G.; Peacock, Wendell; 
Pederzani, Gregory; Peterson, Rolf E.; Pomerantz, Jacob; 
Pratt, Orison S.; Proskauer, Richard; Ream, Donald F.; 
Reichenbach, Hans; Richards, Albert Gustav; Rose, 
Albert; Rossi, Harold H.; Roth, Raymond E.; Rothstein, 
Jerome; Rovner, Leopold; Sable, Morris; Sachs, Donald 
Charles; Saunders, Bernard G.; Schuman, Robert Paul; 
Sedwitz, Ira M.; Smith, Paul Kenneth; Sondheimer, Ernst 
Helmut; Stafford, Jerome W.; Stanton, Henry E.; Stock, 
Charles R.; Stowens, Bernard H.; Sutro, Peter John; 
Swearingen, Thomas; Taxwood, Vincent C.; Taylor, James 
E.; Taylor, William J.; Tuerck, William, Jr.; Turkevich, 
John; Vargus, Joseph Anthony, Jr.; Vine, Benjamin H.; 
Wack, Paul E.; White, A. H.; Whitehead, Walter Dexter; 
Wightman, Arthur Strong; Wikstrand, Wm. Carl; Williams, 
Neal T.; Williams, S. E.; Wilson, William Solomon; 
Withrow, Robert B.; Witte, Robert S.; Wolfe, Ralph A.; 
Wylde, Ronald James; Zihlman, Charlotte. 


KarRL K. Darrow 
Secretary 


Contributed Papers 


Al. The Small Gulf Gravimeter. ANTHONY CARNVALE, 
Gulf Research & Development Company.—For greater 
portability under difficult field conditions, a small and 
light modification of the Gulf Torsion Spring Gravimeter' 
has been constructed. The instrument consists essentially 
of a helical ribbon with a mass hanging freely at the lower 
end. Around the moving system are built the elaborate 
details needed to insure the necessary precision. These 
include the proper choice of spring material, temperature 
control, optical system, damping, magnetic shielding, and 
the elimination of air buoyancy errors. The salient features 
of the instrument and the extreme stability requirements 
of the moving system are discussed. In geophysical pros- 
pecting with a gravimeter, it is desirable to measure gravity 
differences from point to point to a precision of 107’ g. 
The probable error of a single observation with the Gulf 
_ Gravimeter in routine field use is about 4X10~* g. For 
stationary set-up (e.g., for measuring variations in tidal 
forces) the probable error is 1X 10~® g or better. 


1R. D. Wyckoff, Geophys. 6, 13-33 (1941). 


A2. A Thermopile for Micro-Calorimetry. L. S. Mason, 
University of Pittsburgh. (Introduced by A. G. Worthing.) 
—In the study of the thermodynamics of solutions it is 
desirable to measure heat effects (of the order of 0.001 cal.) 
which occur during the dilution of dilute solutions. In the 
Lange-type calorimeter' a thermopile of approximately 
1000 junctions is used differentially. For this purpose a 
thermopile has been constructed using constantan- 
chromel-P couples? mounted in a sheet of Micarta, 5” x 10” 
x #”’, and electrically insulated with a Vinylite resin which 
has an insulation resistance greater than 25 megohms. 
This instrument appears to offer somewhat better be- 
havior than previously used units with respect to stability, 
sensitivity, and rapidity of response to temperature 
changes. The thermopile consists of two sections of equal 
numbers of junctions which are connected in parallel to a 
Leeds and Northrup HS galvanometer in series. At a 
scale distance of nine meters the sensitivity of the system 
is 0.25 microdegree per mm. Since one liter of diluent 
(water of dilute solution) is contained in each half of the 























calorimeter, this corresponds to a mean sensitivity of 
2.5 10-4 cal./mm, with a probable error of the mean of 
+0.03. The heat conductivity constant of the calorimeter 
is 0.004 min.~! and that of the thermopile is 0.03 min.“ 
Measurements with the present thermopile of the heats of 
dilution of sodium chloride solutions agree satisfactorily 
with well-established values. 


1 E, Lange and A. L. Robinson, Chem. Rev. 9, 89 (1931). 
2H. F. Launer, Rev. Sci. Inst. 11, 98 (1940). 


A3. Amplifier with Logarithmic Response. J. A. Hi1pPLe 
AND D. J. GROVE, Westinghouse Research Laboratories.— 
An amplifier with d.c. output voltage logarithmically 
proportional to the input signal over a range of three 
decades has been developed. This response characteristic 
is based on the charging and discharging of a condenser 
through a resistance and is thus independent of the charac- 
teristics of the individual tubes in the circuit. At any fixed 
setting the d.c. output is constant to better than 0.2 per- 
cent. The condenser discharge is initiated periodically by 
means of a trigger circuit which also synchronizes a saw- 
tooth oscillator with the exponential circuit. The d.c. level 
of this saw-tooth is varied by means of the input signal 
so that a thyratron can be made,to fire at any time from 
the beginning to the end of the interval. The impulse 
developed by the firing thyratron drives the grid of a 
hard vacuum tube from very negative to positive for an 
instant. The plate supply for this tube is the exponential 
voltage from the condenser-charging series resistance so 
that for the instant the grid is driven positive a signal 
appears across the load resistance equal to the voltage of 
the exponential at that instant. The output thus consists 
of a series of pulses of magnitude logarithmically propor- 
tional to the input voltage. These pulses are then averaged 
and filtered to d.c. 


A4. Deflection of.a Beam of Cesium Atoms by Gravity. 
I. ESTERMANN, O. C. Srmmpson, AND O. STERN, Carnegie 
Institute of Technology.—In the course of experiments for 
the exact measurement of the Bohr magneton! it was found 
necessary to determine the actual velocity distribution in a 
beam of cesium atoms 2 meters long with a cross section of 
2 mm by 0.02 mm. In the apparatus used, this velocity 
distribution was measured by gravity deflection. The results 
showed that deviations from the Maxwell distribution in 
the range of low velocities occur at much lower oven pres- 
sures than expected. The slow atoms remain longer 
in the neighborhood of the oven slit and have therefore a 
higher probability to be deflected from the beam by colli- 
sions. The results agree with an approximate calculation of 
this effect. 


10. Stern, Phys. Rev. 51, 852 (1937). 
2 Compare V. W. Cohen and A. Ellett, Phys. Rev. 52, 502 (1937). 


AS. Performance Tests on the Penn State Type of 
G-M Circuits. HAroLD E. WALCHLI, The Pennsylvania 
State College—Descriptions of the Penn State type of 
Geiger-Mueller circuits for the measurements of x-ray 
intensities were published in February, 1944 in The Review 
of Scientific Instruments. It is the purpose of the present 
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paper to give the results of performance tests on such cir- 
cuits. These tests are: 1. G-M high voltage supply, tem- 
perature effects, effects of changes in input voltage, 
effects of rheostat settings, improvements in circuits 
resulting from above tests. 2. Power and stabilized-bias 
supply, temperature effects, effect of load, voltage regula- 
tion, time constants. 3. Main amplifier, changes from orig- 
inal circuit, temperature effect, time constants. 4. Pre- 
amplifier, time constant, impedance, change in fixed bias 
vs. plate current. 5. Auxilliary amplifier, temperature effects 
and improvements in the circuit intended to reduce heat- 
ing effects. 6. Complete circuit, voltage recovery of G-M 
supply, oscillograph records of pulses through G-M 
counter. 


A6. The Influence of Statistical Fluctuations on the 
Observed Size-Frequency Distribution of Bursts in a 
Single Ionization Chamber. L. G. Lewis ANp R. Haypen, 
The University of Chicago.—A particular burst of ioniza- 
tion occurring in an ionization chamber will be recorded 
as a burst larger or smaller than its true size by an amount 
equal to the statistical fluctuation in the non-burst particles 
occurring within the collecting time of the ions generated 
by that burst. This phenomenon causes the shape and 
absolute values of the observed size-frequency distribution 
curve to differ from the true curve by an amount given by: 


fla)= fF) exp [—b(x—a)* dx. 


f(a) is the observed differential size-frequency distribution 
function and F(x) is the true function. The constant } is 
obtained from the conditions of observation in the ioniza- 
tion chamber. The true distribution F(x) corresponding to 
the observed size-frequency distribution curve for a single 
chamber given by Lewis and Lewis! has been determined. 


1L. G. Lewis and E. W. Lewis, Phys. Rev. 65, 63A (1944). 


A7. The Ratio of the Average Current Generated by 
Cosmic-Ray Air Showers to the Total Current Observed 
in the Ionization Chamber. E. W. Lewis AND L. G. Lewis, 
The University of Chicago.—The average ionization current 
generated in a single unshielded ionization chamber by 
air showers producing bursts of more than 50 particles has 
been computed from the size-frequency distribution curve 
given by Lewis and Lewis! for an elevation of 10,650 ft. 
The average shower current was found to be 1.810 
ampere. The current produced by all ionizing radiation in 
the chamber was 2 X 10-" ampere. This value was obtained 
by measuring the voltage drop produced by the ionization 
current flowing through a 10"-ohm resistor. The ratio of 
these currents is then 0.9 percent. 


1L. G. Lewis and E. W. Lewis, Phys. Rev. 65, 63A (1944). 


A8. Calculations on Extensive Cosmic-Ray Showers in 
Air. LINCOLN WOLFENSTEIN, University of Chicago. 
(Introduced by M. Schein.) —Large cosmic-ray showers in 
air, investigated with ionization chambers and coincidence 
counters, have been explained hitherto as originating from 
primary electrons of very high energy. The recent measure- 
ments by Lewis! of the frequency of coincident bursts in 
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two unshielded ionization chambers provide a more com- 
plete picture of these showers. The theoretical coincidence 
frequency as a function of the particle density and the 
separation of the chambers has been calculated using the 
cascade theory of showers and the theory of the lateral 
spread of showers.? The calculations give much smaller 
values than the experiments for (1) the absolute value of 
all the frequencies observed, (2) the drop in frequency as 
the chambers are separated, and (3) the increase in fre- 
quency with elevation of bursts in a single chamber. These 
discrepancies strongly indicate that the bursts observed 
in the experiments are not due to showers from primary 
electrons. It is considered possible that the bursts may be 
explained by showers originating not at the top of the 
atmosphere, but relatively close to the chamber. 


1L. Lewis, Phys. Rev. 65, 63A (1944). 
2 J. Roberg, Phys. Rev. 61, 735A (1942); 62, 304A (1942). 


A9. Theory of Counter Experiments on Collision Elec- 
trons Ejected in Air by Mesotrons. ANATOLE ROGOZINSKI, 
University of Chicago.—It was found in recent counter-tube 
experiments! that mesotrons are associated with particles 
present at several meters distance from the mesotron 
trajectories. These particles together with the mesotron 
possibly belong to an Auger-shower, but at least part of 
them are collision electrons (CE) ejected by the mesotron 
along its path in the air. The results of a theoretical calcu- 
lation show that: (a) Only mesotrons of momenta p> po 
~5X 108 ev/c can eject in air CE sufficiently energetic 
to reach and to traverse a counter (C) placed at a distance 
D (a few meters) from the mesotron trajectory. (b) The 
number v of CE, which arise from a single mesotron and 
which are able to discharge (C), is practically independent 
of the mesotron momentum 9, if the spin dependent term 
is neglected, and is given, in first approximation, by 
vo=(BS/xD) X@max(D, E;), where: B=10-? for standard 
air, S=effective area of (C) in m?*, D is expressed in meters, 
and E,=minimum energy necessary for an electron to 
traverse the counter wall; for E;=2 Mev, e.g., the function 
Oinax decreases from 0.50 for D=2 m, to 0.30 for D=10 m. 
(c) When the mesotron spin (0 or 4) is taken into account, 
v<4yvo for the lowest mesotron momenta, but approaches 
asymptotically vo as p increases. (d) Shower production in 
air of the CE adds a very small contribution to the count- 
ing rate in (C). 


1 P, Auger and An. Rogozinski, Phys. Rev. 65, 64 (1944); An. Rogo- 
zinski, Phys. Rev., in print. 


Al0. Measurement of High Energy Electrons in A 
Showers. R. E. Lapp, Ryerson Physical Laboratory, 
University of Chicago.—A discussion is given of the energies 
associated with high energy particles (electrons and pho- 
tons) present in the cores of A showers at sea level. The 
approximate values of these energie are obtained from 
an analysis of the number of cosmic-ray bursts found to be 
coincident with A showers as a function of the iron shield- 
ing the ionization chamber in which the bursts occur. 
The A showers are detected by lateral arrangements of 
G-M counter coincidence sets. 


Bl. The Mechanics of the Metal Cutting Process. 
I. Simple Straight-Edged Tool and Continuous Chip. 
M. EvuGENE MERCHANT, Cincinnati Milling Machine 
Company.—The analysis of the mechanics of a simple 
straight-edged cutting tool removing a continuous chip 
in the absence of a built-up edge, and traveling in a direc- 
tion perpendicular to the cutting edge, as introduced by 
Ernst and Merchant,! has been expanded and correlated 
with experiment. In this special case, the chip is formed by 
a process of shear in a direction perpendicular to the cutting 
edge, in a plane extending from the cutting edge to the 
work surface, and the resultant force on the tool or on the 
work is coplanar with all its significant components and 
with the direction of relative motion of tool and work. 
Relationships pertinent to the geometry of chip formation 
are derived. Significant relationships are derived between 
the various force components and such quantities as the 
shear strength of the metal, the coefficient of friction be- 
tween chip and cutting tool, and the work done in shearing 
the metal and in overcoming friction on the tool face. 
Relationships are presented showing how the direction 
of shear of the metal ahead of the cutting edge is deter- 
mined by the geometry of the force system, for certain 
special cases. 


1H. Ernst and M. E. Merchant, Surface Treatment of Metals (Ameri- 
can Society for Metals, 1941), pp. 299-378. 


B2. Anomalous Behavior of the Electrical Resistivity 
of Some Iron-Cobalt Alloys. Sipney SIEGEL AND ROBERT 
McGreary, Westinghouse Research Laboratories—In the 
course of the development of iron-cobalt alloys for mag- 
netic purposes, erratic variations in the resistivity of 
alloys of presumably similar composition and thermal 
history had been observed. Systematic attempts to de- 
termine the cause of these variations were made. The 
possibility of these being due to an order-disorder transi- 
tion was investigated. The equilibrium electrical resistivity 
of alloys containing 42 percent Co, and 30 or 42 percent Co 
plus 2 percent V, was measured over the range 20°C to 
1100°C. The resistance rises monotonically up to about 
950°C, the phase change temperature, where a discontinu- 
ous increase occurs. Above this temperature, the resis- 
tivity again increases, with a smaller constant, temperature 
coefficient. No evidence for an ordered structure was found 
from these measurements. A series of experiments attempt- 
ing to relate the room temperature resistivity of these 
alloys to thermal history showed that the cold-rolled 
alloy has a lower resistivity than the , sumbably equi- 
librium value of the alloy slowly cooled from high tempera- 
tures; the alloy quenched from above the point has higher 
resistivity than the cold-worked alloy, but still is lower 
than the annealed material. Cold work following quenching 
reduces the resistivity again; repeated quenching and 
annealing result in cyclic decrease and increase in resis- 
tivity. These effects are only present in alloys containing V, 
and are absent in the 30 percent Co alloys. They are diffi- 
cult to account for with any of the usual hypotheses con- 
cerning the resistivity of alloys. 
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B3. Magnetic Permeability at Very Rapid Rates of 
Change of Induction. H. L. GLicK AND SipNEy SIEGEL, 
Westinghouse Research Laboratories ——The permeability of 
several magnetic materials in the form of strip 0.003 inch 
thick or less was measured at rates of change of induction 
as high as 10'° gauss/sec. The measurements were made 
from cathode-ray oscilloscope observations of B and H, a 
square pulse of voltage about one microsecond in duration 
being applied to coils surrounding the sample. During this 
pulse, the induction rises linearly. The eddy current effects 
present in such a situation have been calculated by T. D. 
Holstein and the B-H loops observed are in qualitative 
agreement with those predicted with these calculations 
from the d.c. magnetic properties, the electrical resistivity, 
and the thickness of the strip material. Quantitatively, 
however, the observed loops indicate a permeability, dur- 
ing the rise of induction, about 30 to 50 percent lower than 
the d.c. permeability. This indicates that the mechanism of 
the magnetization process, at rates of change of induction 
so great that the average velocity of domain boundary 
movement must be of the order of 10° cm/sec., is appre- 
ciably different from that obtaining in slow field changes. 


B4. Curie Transition in Ferromagnetic Alloys and Glass 
Containing Nickel. W. R. Ham anp Howarp S. COLEMAN, 
The Pennsylvania State College—The experimental results 
obtained by Bennett and Ham in their studies of hydrogen 
diffusion through ferromagnetic binary alloys showed 
breaks not only associated with a Curie transition for a 
particular alloy, but also Curie transitions corresponding 
to the metals from which the alloy was made. It has been 
suggested that such Curie transitions corresponding to the 
components of the alloy were caused by segregation. 
Magnetic experiments have been performed on the same 
alloys used by Bennett and Ham and the positions of the 
Curie transitions reported by them have been confirmed; 
however, it was found that the intensity of the Curie 
transitions corresponding to the components of the alloy 
were several orders of magnitude less than the Curie of 
the alloy itself. Further experimental work involving 
magnetic torsional measurements of glasses containing 
various amounts of nickel-oxide show a persistence of 
Curie transition at 360°C, even though the nickel-oxide 
content was less than one percent. This latter fact again 
suggests that the Curie transition in some cases is not 
altered by a chemical combination of ferromagnetic ele- 
ments with certain other elements, since it is rather un- 
likely that there would be any appreciable segregation of 
nickel in a glass having a low nickel content. Although the 
experimental work to date does not justify saying there 
could be no segregation in the cases cited, it is felt that the 
evidence is against segregation as an explanation of the 
results. 


BS. Crystal Structure of Tellurium. Joun D. ORNDoFF, 
Purdue University—Results of Hall effect, thermoelec- 
tric power, and resistivity measurements on tellurium have 
been interpreted as evidence! that tellurium exists in two 
crystalline modifications, present in various proportions 
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depending on the temperature and the heat treatment of 
the sample. To decide whether these anomalies in electrical] 
behavior as a function of temperature are really due to the 
development of various modifications, the crystal struc- 
ture of tellurium as a function of temperature has been 
investigated. X-ray diffraction pictures taken with Seeman- 
Bohlin‘ and Debye-Scherrer cameras at temperatures from 
— 192°C to 160°C fail to show any evidence of two modifi- 
cations, nor do they show any peculiar change in lattice 
parameters with temperature. Tellurium evaporated in 
vacuum upon a glass plate produces two types of deposits: 
in contact with the glass a thin homogeneous metallic 
layer of tellurium, on the outside a second film finely dis- 
persed, black and sooty. X-ray investigations of both 
layers give the ordinary tellurium pattern. Electron diffrac- 
tion patterns’ show complete but diffuse electron diffraction 
rings for the metallic film, while the tellurium black‘ yields 
a very sharp orientation pattern. This orientation and the 
non-symmetry of the pattern can be explained as due to 
the fact that the crystals did not have random distribution 
about the orientation axis but tended toward an alignment 
in a preferred direction. 


' Haken, Ann. d. Physik 32, 291 (1910). 

2Smith, Phys. Rev. 1, 352 (1913). 

3 Wold, Phys. Rev. 7, 169 (1916). 

4 Lark-Horovitz, Science 64, 303 (1926). 

5 We are indebted to Mr. T. S. Renzema for taking the electron 
diffraction pictures. 


B6. Spontaneous Time-Dependent Changes in Crystal- 
line Material. D. D’Eustacnio AND SELMA BLAZER 
Bropy, Brooklyn Polytechnic Institute—Freshly broken 
crystalline material undergoes a change during the first 
few hours after being prepared. The phenomenon may be 
described tentatively as a further breaking up of freshly 
broken material. In one experfiment particles of adularia 
(orthoclase) between 200 and 300 mesh in size were photo- 
graphed a few minutes after preparation with a focusing 
camera using unfiltered copper radiation. After an interval 
of several hours, during which the sample was not dis- 
turbed, another photograph was taken. A marked increase 
in the number of spots was observed. The original spots 
were still present. The experiment has been tried with other 
materials. There is good reason to believe that the change 
is associated with the surface rather than the interior of the 
material. 


B7. Operating Characteristics of Ceramic Dielectrics 
with Dielectric Constants in Excess of 1000. R. B. Gray, 
Erie Resistor Corporation.—Within the past several years 
there have been introduced ceramic condensers which have 
extremely high dielectric constants. These materials are 
somewhat similar in their properties to Rochelle salt. 
However, their high dielectric constant obtains over much 
wider temperature and frequency ranges. There are quite 
a great many possible combinations, and curves showing 
some of the temperature and frequency variations are 
presented. Materials are subject to certain aging and volt- 
age effects which produce deviations from the nominal 
capacity. The temperature variation of capacity usually 
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exhibits a peak which may be explained by referring to the 
thermal expansion curve of the material. More exact 
explanation in terms of composition and crystal structure 
must be deferred until a later date when conditions will 
permit such discussion. Materials are somewhat affected 
by humidity and at least a certain fraction of them is water 
soluble so that the production of commercial condensers 
requires fairly stringent waterproofing methods. Probably 
the greatest difficulty with these materials is the fact that 
they are subject to ionic conduction if certain trace im- 
purities are present. This results in failure on life test, 
particularly at slightly elevated temperatures. Data show- 
ing the effects of several impurities will be presented. 


B8. A Resonant Cavity Method for Measuring Dielec- 
tric Properties at Ultra-High Frequencies (50-1000 
Megacycles). C. N. Worxs, T. W. DaAKIN, AND F. W. 
BocGs, Westinghouse Research Laboratories.—A re-entrant 
cylindrical cavity has been adapted to measure the dielec- 
tric constant and power factor of small disk samples of 
insulating materials. The disk of dielectric is placed be- 
tween the flat ends of the cylindrical posts which form the 
re-entrant part of the cavity. The circuit is analogous to a 
less than }-wave-length section of air-filled coaxial trans- 
mission line shorted at one end and ending in a capacitor 
with dielectric at the other end. This cavity with the dielec- 
tric inserted is brought to resonance by tuning the fre- 
quency of a loosely coupled oscillator. The power factor of 
the dielectric is determined from the breadth of the reso- 
nance curve and the dielectric constant is obtained by 
removing the disk of dielectric and decreasing by a meas- 
ured amount the separation of the flat ends of the re- 
entrant posts until resonance is again obtained at the same 
frequency as previously. The precision differential screw 
for varying the extension of the re-entrant post and other 
mechanical details are described. Electrical coupling and 
detecting methods are also discussed as well as the theory 
of the measurements. Each individual cavity can be 
utilized only over a frequency range of about 1.5 ratio from 
lowest to highest frequency. Therefore, several different 
sized cavities would be required to cover a range of fre- 
quency from 50 to 1000 megacycles. An accuracy of 
+0.00005 in tan 6 and +1 percent in dielectric constant 
may be obtained in routine measurements. Because the 
cavity has a very high Q (>2000), it is much more sensitive 
to low power factor dielectric samples than any conven- 
tional coil and capacitor resonant circuit. Operation of the 
apparatus is simple and rapid and very similar to the 
susceptance variation technique now used at lower radio 
frequencies. 


B9. Decomposition of Water with the Permanent 
Magnet. F. EnRENHAFT, New York City—Permanent 
Alnico magnets, ‘Blue Streak,”’ are used. They are semi- 
circular rings with outer radius of 5.1 cm and inner radius 
of 1.58 cm, 1.75 cm thick. These were fitted with pole pieces 
of pure soft Swedish iron having circular cylindrical ex- 
tensions 1 cm in diameter, 12.5 mm long, terminating in a 
truncated cone, the end 6 mm in diameter. The extensions 
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passed through rubber gaskets into a glass cell divided by a 
partition so that the gas evolved from either pole could be 
collected separately. This partition had four holes 1 mm in 
diameter plugged with asbestos. The pole faces were about 3 
mm apart. The field strength was much less than that in 
the earlier described experiments with the electromagnet.* 
Using 4 percent H.SO, in the cell, oxygen was found to be 
coming from both poles together with the hydrogen formed 
by chemical action. More oxygen was found to be coming 
from the north pole than from the south. For instance, 1.43 
percent and 1.12 percent in one experiment, 0.87 percent 
and 0.57 percent in another. Using the pole pieces alone 
without applying the magnet, 0.00 percent oxygen was 
found. Measuring the pole strength of one magnet by 
means of a search coil and ballistic galvanometer before 
and after one magnetolysis, the decrease in strength was 
from 16,980 to 15,380 maxwells. The flux value of this 
magnet was at a steady state before and after this test. The 
evidence here given shows that water can be decomposed 
by a permanent magnet as well as by the electromagnet, 
that magnetism is polar, and that the energy for this process 
is supplied by the magnet as shown by its loss in pole 
strength. 


* Felix Ehrenhaft, Phys. Rev. 60, 169 (1941); 61, 733 (1942); 63, 
216, 461 (1943); 64, 43 (1943), etc. 


Cl. Magnetic Potential Energy and the Zeeman Effect. 
F.W. WARBURTON, University of Kentucky.—The quantiza- 
tion of orbital angular momentum in units of 4/2x remains 
unchanged in obtaining the Schroedinger equation from the 
Hamiltonian function which includes magnetic potential 
energy. When intrinsic electron spin is omitted and the 
magnetic potential energy of an orbit used instead, the 
vector model description of the Zeeman effect needs re- 
vision. Bohr’s original assumption that the magnetic 
field affects the mechanism of transitions between energy 
levels, together with the mutual electrostatic oscillating 
dipole force of neighboring atoms during a transition, 
provides a qualitative explanation of certain anomalous 
Zeeman patterns, and it avoids the incongruity of includ- 
ing classical magnetic force in setting up the energy states 
then suspending this force during a transition between 
states. 


C2. Sensitized Fluorescence of Copper and Test of 
Selection Rule for Collisions of the Second Kind. F. K. 
McGowan AND J. G. Winans, University of Wisconsin.— 
Sensitized fluorescence and Tesla discharge spectra of 
copper vapor were obtained with the simplified procedure 
used for other elements.' Particles of copper were heated 
to 1377° without collapsing the fused quartz tube. The 
Tesla discharge spectrum showed 28 copper lines and the 
sensitized fluorescence showed four. These were the reso- 
nance doublet 3273 and 3247 and lines at 2492 and 2441. 
Microphotometer comparison of intensities of lines in 
sensitized fluorescence and arc showed 2492 greatly en- 
hanced in fluorescence. The selection rule? found for colli- 
sions of the second kind for Sn and Pb was found to apply to 
copper excited by collision with Hg *Po. The rule states 
that for atoms A and B colliding J=Ja+Jp remains 












constant during collisions of the second kind. This rule 
requires Cu 3273 to be stronger than 3247 for excitation 
by Hg *Po. 3247 is stronger than 3273 in the arc spectrum. 
Fluorescence excitation by Hg 6 *P» alone was produced by 
introducing a time lag of about 10~ sec. between excitation 
and emission. 


1J. G. Winans and R. M. Williams, Phys. Rev. 52, 930 (1937). 
2J. G. Winans, Phys. Rev. 60, 169A (1941). 


C3. Raman Spectra of Normal-heptane and Normal- 
hexane. Forrest F. CLEVELAND AND DONALD E. LEE, 
Illinois Institute of Technology.—Continuing work on 
hydrocarbons, Raman frequencies, relative intensities, 
and depolarization factors have been obtained for n-hep- 
tane and n-hexane. The relative intensities and depolariza- 
tion factors were obtained by use of a Gaertner microden- 
sitometer. The Raman frequencies will be compared with 
the previous results of Kohlrausch and Képpl and of 
Rosenbaum, Grosse, and Jacobson. No previous polariza- 
tion data seem to exist and the previous relative intensities 
are only visual estimates. The results will be discussed in 
the light of the selection rules for various possible equi- 
librium configurations. For example, if m-hexane has a 
symmetry corresponding to the point group C;, 27 of the 
54 fundamentals are forbidden in the Raman spectrum 
and all of the allowed Raman fundamentals should be 
polarized. Furthermore, for this symmetry, combination 
frequencies of two A, or two A, frequencies, overtones of 
A, frequencies, and first, third, fifth, etc. overtones of A, 
frequencies, are all allowed in the Raman spectrum and 
could bring about Fermi resonance interactions if they 
should fall sufficiently close to one of the type A, Raman 
lines. 


C4. Some Double-Bond Absorptions in the Vacuum 
Ultraviolet. JoHN R. Piatt, Northwestern University.— 
Recent measurements in heptane solution! are compared 
with the theoretical predictions of Mulliken? and McMurry® 
on the absorption below 2000A of conjugated and uncon- 
jugated C=C and C=O groups, with very satisfactory 
agreement. The data shed light on the theory of —-COOR: 
the first true maximum above N—Q (54,000 cm) must be 
above 66,000 cm= (by extrapolation) with f>0.1, and so 
is probably the N—V transition. The splitting and red 
shift of the 1850A ethylenic peak with increasing number of 
double bonds (unconjugated) may be attributable to 
hyperconjugation (f values increasing almost linearly, as 
predicted), but the reason for the red shift (about 1000 
cm of the trans-ethylenes from the cis-ethylenes is not 
clear. Secondary maxima near 55,000 cm™ in the conju- 
gated acids and dienes may be the N-—> V2,; transition. 

1 Platt, Rusoff, and Klevens, J. Chem. Phys. 12, 535 (1943); Rusoff, 
Platt, Klevens and Burr, in preparation. 
2? Mulliken, Reports on Progress in Physics (Physical Society of 


London, 1941), Vol. VIII, p. 231. 
+H, L. McMurry, J. Chem. Phys. 10, 655 (1942). 


CS. Apparatus for the Determination of Dispersion at 
Supersonic Frequencies. L. N. LIEBERMANN, University of 
Kansas.—Variation of the velocity of sound with frequency 
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yields useful information on the rate of transfer of trans- 
lational energy to molecular vibrational energy. Variations 
have been measured by a method which allows the direct 
observation of small velocity changes. The apparatus 
permits comparison of the phases at various frequencies 
in two gases, one being a non-dispersive gas. If the phases 
are measured in the two gases at positions which correspond 
to equal numbers of wave-lengths from the source, then 
A¢/¢=AV/V where A¢ is the phase difference in the two 
gases and V is the velocity in the dispersive gas. A 90- 
kilocycle quartz bar is excited at its fundamental and upper 
harmonics to produce the various frequencies. The sound 
is received by an L-cut Rochelle salt microphone and the 
phase determined by a 360° electrostatic phase shifter 
accurate to +3°. @ was 3600° (10 wave-lengths) at the 
fundamental. V is taken from the literature. Five fre. 
quencies between 90 and 450 kilocycles were used; these 
differed by approximately 90 kilocycles. Insufficient 
amplifier gain at higher frequencies prevented extension of 
this range. Measurements on the dispersion of COs, at 
atmospheric pressure are in agreement with previous 
values. 


C6. Precision Measurement of the Velocity of Sound at 
Supersonic Frequencies Using a Microphone. E. J. Grorn 
AND L. N. LIEBERMANN, University of Kansas.—The 
most precise measurements on the velocity of sound at 
supersonic frequencies have been made with the Pierce 
interferometer. The present method which is adapted from 
audible frequency technique offers greater precision to- 
gether with certain advantages. The horizontal plates of a 
cathode-ray tube are connected to a radiofrequency oscil- 
lator which excites a quartz bar; the vertical plates are 
connected through a video amplifier to a microphone. 
The ellipse formed indicates wave-length as the micro- 
phone position is changed. The microphone is an L-cut 
Rochelle salt crystal. The microphone screw, 8 cm long, 
is calibrated to 0.001 mm. Frequencies are measured to 1 
part in 10° using standard equipment. Measurements made 
on the velocity of sound in air (uncorrected for CO, and 
H.0 content) at 0°C indicate a value of 332.99+0.01 
m/sec. at a frequency of 477.72 kilocycles/sec. Values at 
other frequencies are obtained by exciting the quartz at 
various harmonics. .The apparatus can readily be adapted 
to measure sound absorption. 


C7. The Kinetics of the Thin Oxide Film Formation of 
Iron and Copper. EarL A. GULBRANSEN, Westinghouse 
Research Laboratories. (Introduced by S. Siegel.)—The 
oxidation measurements are obtained by means of a sensi- 
tive vacuum microbalance. The temperature regions of 100 
to 200°C for copper and 200 to 450°C for iron are covered. 
The pressure is varied from 0.2 to 0.0002 atmosphere of 
oxygen. Several types of empirical equations are shown to 
fit the data. The parabolic equation is found to fit after a 
certain initial period of the reaction. The form of the 
equation can be correlated with diffusion phenomenon. 
The influence of temperature pressure and surface effects 
are discussed. 
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C8. Mechanism of Electrolytic Oxidation of Aluminum. 
Scott ANDERSON, Carleton College—When an aluminum 
sheet is made the anode in certain electrolytes, a layer of 
oxide forms on its surface and constitutes a barrier between 
the metal and the electrolyte. Further oxidation can occur 
only as a result of (1) the aluminum or (2) oxygen ions 
penetrating this film. A reaction mechanism is proposed in 
which one of three aluminum ions diffuses from the metal 
to the electrolyte and the other two occupy aluminum sites 
in an extension of the oxide lattice. The oxygen sites so 
created are filled by ions from the oxide and the oxygen 
ion vacancies diffuse to the oxide-electrolyte interface. 
The mechanism is successful in explaining the strong ad- 
herence of the films, and in predicting the experimental 
ratio of the weight of oxide formed to the weight of metal 
removed, as well as other properties. 


C9. A New Research Tool—Irreversible Solution 
Potential Measurements. R. H. Brown AND R. B. MEArs, 
Aluminum Research Laboratories—Reversible solution po- 
tential measurements are contrasted with irreversible solu- 
tion potential measurements. The advantages and limita- 
tions of each are discussed. The importance of selecting 
suitable conditions, especially methods of surface prepara- 
tion and composition of the solution used for making the 
measurements, is emphasized. It is pointed out that the 
limitations of the method are largely dependent on the type 
of solution used. The interpretation of results obtained 
using solutions producing relatively high anodic polariza- 
tion is presented and compared to results from solutions 
causing relatively high cathodic polarization. Methods are 
given for applying irreversible solution potential measure- 
ments to studies of: (a) solid phase diffusion of metals, 
(b) solid phase metallurgical transformations, (c) distribu- 
tion and mechanism of corrosion, (d) cathodic protection 
and galvanic action, (e) metallic and nonmetallic coatings. 
Each of the above is illustrated by specific examples. In 
most cases the results of solution potential measurements 
are correlated with results obtained by other methods. 

1 Chemical Metallurgy Division, Aluminum Research Laboratories. 


* Chief, Chemical Metallurgy Division, Aluminum Research Labora- 
tories. 


D1. Deviations in the Classification of Stable Nuclei. 
J. D. Kursatov, The Ohio State University.—One case of 
deviation from the rule imposed in the classification is in 
even nuclei of period four, J=20. The mass number of the 
last nucleus of the group is higher than that of the first 
of the following group, J=28. This is due to 6.5m". The 
isobaric nucleus ¢oNd'* and the next in the sequence 
5m" of P, are said to be radioactive. Other gaps in 
groups of even nuclei are found mainly in periods 6 and 8. 
The deviations in odd nuclei are in groups of J=1, 11, 23, 
and 25 and are due to the following shifts: 


Actually existing 


Non-existent stable nuclei expected 
stable nuclei 


by the classification 


Species I Species I 
isA” 3 wk 1 
«Ma” il «Mo” 13 
ssCel? 23 s7Lal™ 25 
s1La!*? 23 ssBal? 25 
all 23 cooNd'# 25 
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The position of s;La™® remains unique in the classification. 
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This odd-even nucleus is located in the midst of a group of 
the even-odd series. Three of four odd isobars could not be 
placed in the classification. 


D2. Classification of Stable Nuclei by Means of 
Particles of Mass Numbers Three and Four. M. H. Kursa- 
TOV AND J. D. Kursatov, The Ohio State University.—The 
requirements imposed in the classification are: (1) that the 
difference between isotopic numbers of groups of nuclei 
having common characteristics remain constant and (2) 
that the mass number of the last nucleus in a group be no 
higher than that of the first of the following group. The 
odd-even and even-odd series of nuclei form 22 groups 
which satisfy these requirements. With the even-even 
unequal series the conditions are fulfilled by the division 
of its 22 groups into four periods: P:, J =8a+2, A =4b+18: 
P., [=8a'+6, A=4b’'+46: Py, I[=8a+4, A=4b"+36: 
Ps, [=8a’'+8, A=4b'’+48: where a=0,1---5; a’= 
0,1---4;b=0,1 --- 47; b'=0,1 --- 38; 6”=0,1 --- 43; 
b’’ =0, 1 --- 39. AJ =8 is constant and the periods 2 and 6 
and 4 and 8, respectively are isobaric. The nuclei of a given 
sequence are built up by the addition of the group (22m). 
The use of the particles (p2”) and (m2p) results in the 
continuity of groups of even-odd and odd-even nuclei. 
A group of even-odd nuclei always precedes one of odd- 
even nuclei having the same number of particles of mass- 
number three. The isotopic number of even-odd nuclei, 
therefore, is two units lower than that of odd-even nuclei 
in the same sequence. A periodicity of nuclei with common 
characteristics is observed in the appearance of odd se- 
quences alternating with consecutive periods of even nuclei. 


D3. Radioactive Silver. L. K. Hurst anp M. L. Poot, 
The Ohio State University.—In the continuance of the study 
of the radioactive isotopes of silver, variations from former 
results* were observed. Bombardments of silver with high 
energy deuterons have given an extremely long silver period 
of the order of 225 days. Results have been obtained to 
support the postulation of K-electron capture. X-rays of 
0.6A which are close to the K-a x-rays of Pd and a gamma- 
ray of 1.0 Mev which could be expected from the excited 
Pd have been observed. Low intensity electrons of 0.5 
Mev were also obtained. During this study neutrons which 
are produced from deuteron bombardment of silver were 
observed to be energetic enough to produce the n, 2m reac- 
tion leading to the 24.5-minute Ag’ isotope. In conse- 
quence of this observation the production of the Ag’ 
isotope by the d, H* or d, p2n reaction* must be seriously 
questioned. 


* R. S. Krishnan, Proc. Camb. Phil. Soc. 36, 500 (1940). 


D4. Transmutation of Titanium. Cart T. HIBpon, 
M. L. Poot, anp J. D. Kurspatov, The Ohio State Uni- 
versity.—A study of the 16-day period of V“ produced by 
activation of titanium with deuterons has led to a re- 
evaluation of the energies of the positrons and gamma-rays. 
The positron upper limit has been found to be 0.58 Mev. 
Lead absorption measurements indicate a gamma-ray of 
1.50 Mev. The probability of V“ decaying by K-electron 
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capture to that of positron emission is approximately 18 
to 1. The scandium fraction from titanium bombarded 
with deuterons indicates that the (d, a) reaction is un- 
usually strong. The isomers of Sc“ and the isotopes Sc“, 
Sc, and Sc*® were produced. It has been found that the 
yield of the 85-day period of Sc“ from the Ti* (d, a) 
reaction is greater than from the Sc* (d, p) reaction. Two 
new reactions, Ti*® (d, a) and Ti®™ (d, a) are established. 
The previously reported 72-day period of Ti* was looked 
for but has not been produced in observable quantities. 
Alpha-particle bombardments of titanium free of calcium 
have failed to produce the 3.7-hour period of V*. 


DS. The Regularities in Energy Levels of Light Nuclei. 
W. Y. CHANG, Princeton University.—It is found, from the 
energy levels of about fifty nuclei, that in each of the four 
series, 4n, 4n+2, 4n+1, the nuclei with the same value 
of the isotopic number N-Z have similar level systems, and 
that the general distribution of the levels reveals a strong 
dependence of the level density on the values of N-Z. 
These two points have also been discussed from the stand- 
point of the ‘“‘alpha-particle model.’’ In each series the 
levels of the nuclei having the same value of N-Z can be 
well expressed by the following formula: E=as—b/s, 
where s is an integer, and a and b are two constants different 
for different series and N-Z. It is found possible to interpret 
the first term as the contribution of the nuclear particles 
as if they were not interacting with one another, and the 
second term as the correction for the interactions. The 
results stated above are for the levels below dissociation 
energy. The data for those above dissociation energy is 
insufficient to make any correlation. However, one point 
is noted: In a given group of nuclei (i.e., of the same N-Z) 
the energies of the level systems seem to decrease with 
increasing mass number as would be expected generally 
from theoretical consideration. 


D6. Excitation Functions Obtained by Bombardment 
of Sodium Bromide with 14-Mev Deuterons. E. T. CLARKE 
AND J. W. IRVINE, JR., Massachusetts Institute of Tech- 
nology.—Excitation functions for production of the radio- 
active isotopes Na™ (14.8 h), Br® (34), and Kr (34h) 
were obtained by bombardment of NaBr films obtained by 
evaporation in high vacuum onto aluminum foil. Deuteron 
energy determinations were made by measurement of 
range in aluminum, using the relationship given by Living- 
ston and Bethe.' The three activities were separated by 
solution of each foil in water to liberate the krypton, and 
precipitation of the bromine as AgBr, leaving the radio- 
sodium in solution. The activities of each fraction were 
measured on a calibrated gamma-ray counter, which 
yielded three excitation curves and the absolute disintegra- 
tion rates for the Na* and the Br® whose disintegration 
schemes are known. The reaction Br*®!(d, ») Br* exhibits an 
excitation function which begins at about 3 Mev, rises toa 
maximum absolute cross section of 3.8 10-** cm? at 8.5 
Mev, and decreases with further increase in energy. The 
excitation of the reaction Na**(d, p)Na™ shows the same 
shape but begins at about 1 Mev and reaches its maximum 
value of 4.7 X 10-** cm? at 6 Mev. The reaction Br(d, 2n)Kr 
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shows a definite threshold at 5.3 Mev and rises almost 
linearly with increase in deuteron energy. 
1M. S. Livingston and H. A. Bethe, Rev. Mod. Phys. 9, 261 (1937), 


D7. The Variation of the ~-Ray Activity of Radioactive 
Minerals with Particle Size. F. MorGan, Gulf Research 
and Development Company.—In oil well logging, carnotite 
or other radioactive material is frequently added to cement 
or mud in order to determine subsequently the location of 
the mixture by y-ray counters or equivalent devices. To 
permit a quantitative interpretation of the resulting y-ray 
logs, it is necessary to know the strength of the y-ray 
source. This information also permits a choice of the 
optimum quantity of the radioactive tracer to be added to 
the cement or mud. The strengths of a variety of radio- 
active minerals, including a number of carnotites, were 
accordingly measured. As was expected from previous 
field observations, the results show large differences in 
activity. For example, one carnotite sample was found to 
have at least six times the strength of another of approxi- 
mately the same cost. For the carnotites tested, the activity 
was found to be a function of particle size of the mineral, 
In most cases the activity was a minimum at an intermedi- 
ate particle size, but for all the samples the finest particles 
showed by far the greatest activity. Samples of gummite 
and uraninite gave a more uniform distribution of activity 
with particle size. These differences can be correlated with 
the geological origin of the minerals involved. 





D8. Charge Distribution around Protons. J. M. Biarrt, 
Cornell University —The calculations of Frohlich, Heitler, 
and Kahn! for the deviation from the Coulomb law for a 
proton due to mesons are reexamined and extended to the 
scalar and pseudoscalar meson theory. It is sufficient to 
calculate the charge density around the nucleon. The 
potential energy U(r) for a test particle is then obtained 
as the solution of Laplace’s equation, AV =4ep. A per- 
turbation calculation is used and only the terms propor- 
tional to the square of the coupling constant are considered 
It is found that in all standard theories (scalar, pseudo- 
scalar, vector) the total charge diverges with some nega- 
tive power of r. Some physical meaning may be attributed 
to the charge density in a distance larger than the Compton 
wave-length of the proton. The same problem was in- 
vestigated using the recent theory of Dirac* which involves 
negative energy states of the mesons. It turns out that the 
total charge is exactly zero for all types of meson theories 
if the recoil of the nucleon is neglected. It has a finite 
value when the recoil is taken into account. 

' Fréhlich, Heitler, and Kahn, Proc. Roy. Soc. A171, 269 (1939). 


?P. A. M. Dirac, Proc. Roy. Soc. A180, 1 (1942); W. Pauli, Rev. 
Mod. Phys. 15, 175 (1943). 


El. Paths of Electrons and Ions in Non-Uniform 
Magnetic Fields. N. D. COGGESHALL AND M. Muskat, 
Gulf Research and Development Company.—The integration 
of the Lorentz force equations to give electron or ion 
paths has been reduced to simple quadratures for systems 
in which the electric field is zero and the magnetic field 
is a function of one Cartesian or cylindrical coordinate. 
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For several interesting types of magnetic field variation the 
quadrature can be carried through analytically, and even 
for complicated magnetic fields, or such as are known only 
empirically, the numerical integration can be effected 
without difficulty. From general considerations of the 
functions involved, it is possible to determine the exten- 
sion and periodicity of the orbits for any set of initial 
conditions. The representation used is also convenient for 
obtaining information regarding the dispersion and focusing 
characteristics of the trajectories, such as may be used for 
instrument design. Illustrations of such results will be 
presented, including the manner in which the nature of 
the orbits varies with changing initial conditions. 


E2. Energy Distribution of Electrons Within Dense 
Electron Beams. C. J. Catsick, Bell Telephone Labora- 
tories.—An electron beam injected into an enclosure having 
equipotential walls sets up a non-uniform potential dis- 
tribution within the enclosure because of space charge.! 
The non-uniformity increases with beam current until an 
unstable condition is reached when current falls rapidly 
and a virtual cathode forms. Calculations are presented for 
beams and enclosures of two- and three-dimensional sym- 
metry. In the first case, the maximum possible beam cur- 
rent is i=9.4X 10-*y@*/? when the beam, of half-width y, 
completely fills the enclosed space whose walls are at 
potential @ relative to the electron source. When the virtual 
cathode forms, the maximum current drops to 1=4.65 
X 10~*y@*/2, In the second case, the corresponding values 
are 1=3.8X 10-°@?/? and i= 3X 10-*$* and do not depend 
on beam radius. If the beam does not fill the enclosure, 
the maximum currents are reduced, especially in the 
second case. With dense beams, deflection of an electron by 
electric or magnetic fields depends on its trajectory within 
the beam, and thus deflection distortion of an electron 
“spot’’ may become serious. 


1A. V. Haeff, Proc. I.R.E. 27, 5 (1939). 


E3. Limiting Stable Current in Electron Beams in 
the Presence of Ions. J. R. Pierce, Bell Telephone 
Laboratories.—If an electron current at a given voltage is 
injected into the space between two grids, there is a 
limiting stable value of electron current beyond which a 
virtual cathode will form. There is a similar limiting current 
for electrons restrained by an infinite magnetic field to 
axial motion through a conducting tube. There are similar 
limiting stable currents even if the electronic space charge 
is neutralized by positive ions of infinite mass. Such a 
limiting current is a constant times the limiting current in 
the absence of ions. For parallel grids the limiting current 
density is 

im = 27 2?(e/m)%¢ V 9/8 /x?, 


For a conducting tube of infinite length capped by grids 
the limiting current is 
Im = 23/2(2.4)*x(e/m)* Vo. 


E4. Mass Spectrometer with a Small Magnet. R. E. 
Fox, J. A. Hippte, anp T. W. WILLIAMS, Westinghouse 
Research Laboratories.—A five-inch radius, 60° type mass 
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spectrometer with a unique magnet design has been con- 
structed. The magnet pole faces are sealed into the analyzer 
tube by means of a Kovar cylinder to which the glass is 
sealed, thus establishing the vacuum seals. The magnetic 
flux is transmitted to the pole faces through the Kovar 
cylinder from an external source which is either a perma- 
nent or an electromagnet. Uniform magnetic fields of 
about 1500 gauss in a quarter-inch gap have been obtained 
with an instrument in which the complete analyzer tube 
plus the magnet weighed less than twenty-five pounds. The 
design and construction of the magnet allows the pole faces 
to be accurately aligned to less than 0.001 inch. With a 
larger external magnet it is possible to obtain stronger 
fields, since the size of the magnet is not fundamentally 
limited. Two metal cylinders bolted to the 60° magnet 
faces act as shields for the ion beam and also serve to 
mount the ion source and collector. These metal cylinders 
are so constructed that when bolted to the magnet the ion 
source and collector slits are automatically aligned with 
the magnet gap. . 


ES. The Reversed Cyclotron. LAURENCE ELLSWORTH 
Dopp, University of California at Los Angeles.—A pro- 
cedure based on theory and practicable experimentally 
would be to operate the cyclotron backward. Such reversal 
would permit the slowing down to some desired speed, 
perhaps zero, of an ionic particle having initially a high 
velocity. The principle is applicable theoretically for the 
stopping in this manner of a high speed cosmic-ray ion. 
Whether this is practicable for example in the case of a 
six-billion-volt cosmic-ray electron appears to be a ques- 
tion of attaining experimentally the ultra-high frequency 
of reversal of potential that would be necessary. 


E6. Portable Demonstration Electron Microscope. 
IGor BENSEN, General Electric Company.—A discussion of 
the design features making possible a compact light-weight 
instrument completely contained in two suitcase units. 


E7. Potential Nuclear Monochromatic Electron Sources. 
M. L. Poo, Ohio State University—Among the some 300 
synthesized radioactive nuclei now available a few have 
been found which emit electrons in a narrow energy range. 
Most nuclei, however, emit the well-known continuous 
beta-ray spectrum. The above-mentioned few special 
nuclei are identified with internally converted nuclear 
gamma-rays. It is desirable that the internal conversion 
coefficient be close to 100 percent. Monochromatic electron 
energy bands between 50 and 300 Kev may be had orig- 
inating from certain nuclei with half-lives up to several 
months. Detail data for the isotopes Pb*™*, TI, Ba™*, and 
Ba will be presented as well as a summary of the field. 
Eventual applications of nuclear sources seem probable 
for certain uses in electron microscopes, electron diffrac- 
tion, and special radio tubes. (Invited Paper—20 minutes.) 


E8. Electron Microscope Study of Ferromagnetic 
Domains. L. Marton, Stanford University.—Ferromag- 
netic domains can be investigated by means of the electron 
microscope by using a replica technique. Such replicas 
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can be prepared by suspending a ferromagnetic colloid 
in the solution of a plastic to be cast on the polished surface 
of a ferromagnetic material. After solidifying the plastic on 
the surface of the magnetized sample it can be stripped in 
the same way as any other metallographic replica. Pre- 
liminary experiments show that gelatine is a good sus- 


pending medium; however, it offers difficulties in stripping. 
The preparation of similar colloids in the usual plastics of 
replica technique offer some difficulties also. A variant of 
the method consists in placing a very thin collodion film 
in contact with the sample and placing a droplet of the 
ferromagnetic colloid on the film. After evaporation of the 
solvent the film is lifted and placed in the microscope. 
Patterns obtained on different materials and at different 
degrees of magnetization are shown. 


E9. Studies of Purified Preparations of the MV Strain 
of Poliomyelitis Virus by Means of the Electron Micro- 
scope. H. S. Lorine, C. E. SCHWERDT, AND L. MARTON, 
Stanford University.—Several highly purified preparations 
of the MV strain of poliomyelitis virus have been examined 
in the Stanford electron microscope. The samples were 
prepared by the procedures previously described! with the 
exception that in some cases thesextracts were allowed to 
stand in the refrigerator for several weeks before the solu- 
tion was ultracentrifuged. In biological tests for activity 
with some preparations 100 percent of the animals de- 
veloped poliomyelitis when injected at concentrations of 
5X 10-" ¢g of nitrogen. The micrographs show the presence 
of a relatively homogeneous material. The majority of the 
particles appear spherical or nearly spherical in shape. They 
are also considerably less well-defined than the reference 
particles on the same micrographs. Estimations of the 
diameters of the presumed virus particles range from 15 to 


20 mu. There is, however, a definite indication of asym- 
metry ranging from 1:2 to 1:5 in some of the micro- 
graphs. The possible explanation of this asymmetry, as 
being due to a chain-like aggregation of the smaller par- 
ticles, must be left open at present. No evidence has been 
obtained in any of the preparations for the highly asym- 
metric or thread-like particles reported for the murine 
strain by Gard and Pedersen. 


1H. S. Loring and C. E. Schwerdt, J. Exper. Med. 75, 395 (1942), 


T1. Note on the Inductive Measurement of Magnetic 
Field Strength* F. T. Rocers, JRr.,f The University of 
Houston. To be read by title—Two of the minor correc- 
tions! to which the standard expression! MI/2AN for 
magnetic field strengths is subject, are listed here in a 
convenient approximate form. For small effects of twisted 
leads to the rotating coil, multiply MI/2A N by 


F,=1—a/2AN, (1) 


where a is the small “‘projected area’ of the (equal) ‘‘loops”’ 
in the leads. For small effects of incorrect orientations of 
the coil, multiply MI/2A N by 


F;,=2 sec ¢/[cos (@—«)+cos 6], (2) 



















where « is the small defect of the rotation angle from 180°, 
¢ is the small angle between the axis of rotation and the 
plane (P) perpendicular to the field, and @ is the small 
angle between the plane of the coil in its initial prerota- 
tional position and P. 


* Contribution No. 90 from The Division of Sciences of The Uni- 


versity of Houston. 


t Now with The Lukas-Harold Corporation, Indianapolis. 
1 See F. T. Rogers, Jr., Phys. Rev. 52, 381 (1937) for the definitions 
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MEETING OF THE NEW ENGLAND SECTION HELD AT BosTON, MASSACHUSETTS ON APRIL 15, 1944 


HE twenty-third regular meeting of the New England Section of the American Physical Society 
was held at Simmons College, Boston, Massachusetts, on Saturday, April 15, 1944. Forty-four 
members were in attendance. The invited papers were as follows: 


The Review of the Papers Given at the January 14-15 Meeting of the American Association 
of Physics Teachers at Columbia University. Given by Louisa BILLINGs AND LILLY LORENTz, 
Smith College, AND MILDRED ALLEN, Mount Holyoke College. 

Physics at Simmons College. James M. Hyatt, Simmons College. 

Echo-Sounding by Flying Bats. DonaLp R. GrirFin, Fatigue Laboratory, Harvard University. 

Opportunities in and Training for Industrial Optics. ALLAN E. PARKER, Worcester Polytechnic 


Institute. 


Some Reflections on the Relation of Mathematics and Physics. R. B. Linpsay, Brown Uni- 


versity. 


In addition the following short research paper was presented. 


Magneto-Current Phenomena in Nickel. Henry A. 
PERKINS, Trinity College—Two experiments on related 
magneto-current effects in a nickel bar may be outlined as 
follows. In one of these a bored-out bar was used and the 
effect on residual magnetism, caused by a current sent 
through the metal during the process of demagnetization, 
was compared with that caused by a current flowing 
through a wire run through the hollow center of the 
cylindrical bar. When the current in the nickel was stopped 
after the magnetizing field had been reduced to zero, there 
was an abrupt decrease in the residual flux, but this de- 
crease was much larger when the sense of the current was in 
opposition to that of the flux than when with it. This sug- 
gests that in one direction the current favors the magnetic 
state more than it does when in the opposite direction. 
When the same current flows along the axis and not 
through the metal, there is about the same change of flux 


MILDRED ALLEN 
Secretary-Treasurer 


caused by breaking it, but a much smaller alteration in this 
change when the current’s direction is reversed. 

To see if any difference due to reversing the current 
exists in the often observed magneto-resistance phenom- 
enon, a slender nickel bar, donated (like the one above) by 
the International Nickel Company, was tested in a Kelvin 
double bridge. Its resistance was measured when unmag- 
netized, and then when magnetized to a strength of 1900 
gauss. The proportional increase AR/R due to magnetiza- 
tion was about 0.6 percent, but no difference could be 
detected when either the field or the current was reversed. 
However, an alternating 60-cycle field, of about the same 
intensity as before, produced a decrease in AR/R of an even 
greater magnitude. This was so marked that it easily 
overcame the increase of resistance due to heating when the 
field was applied, although this increase was fairly large. 
When iron was substituted for nickel, neither of the above 
effects could be observed with any certainty. 
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MEETING OF THE SOUTHEASTERN SECTION HELD AT ATLANTA, GEORGIA ON MAY 5 AND 6, 1944 


HE tenth annual meeting of the South- 

eastern Section of the American Physical 
Society was held in the Atlanta Biltmore Hotel 
and at the Georgia School of Technology, 
Atlanta, Georgia, on Friday and Saturday, May 
5-6, 1944. Approximately sixty-five members and 
guests attended the meeting. Local arrangements 
were made by a committee headed by J. H. 
Howey. 


The regular program consisted of eight papers, 
abstracts for five of which are appended hereto. 
Abstracts of the other papers will be found in the 
August, 1944 issue of the American Journal of 
Physics. 

Additional features of the program were a talk 
by H. T. Briscoe of the Manpower Commission's 
Professional and Technical Division of the 
Bureau of Training, on ‘Manpower problems in 
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the sciences;” an address at the annual dinner by 
Homer L. Dodge, Director of the Office of 
Scientific Personnel; and a review by F. L. Brown 
of some of the papers presented at the recent New 
York meeting of the American Association of 
Physics Teachers. 

At the business meeting the election of the 
following officers for the Section was announced: 
Chairman, O. T. Koppius; Vice Chairman, F. L. 
Brown; Secretary, E. Scott Barr; Treasurer, C. 
B. Crawlev; and Member of the Executive Com- 
mittee, Eric Rodgers. 

The time and place for the 1945 meeting is to 


be announced later. 
E. Scott BARR 


Secretary 


1. Transient Arc Discharge Between a Stationary and a 
Rotating Electrode. James C. Barnes, University of 
Virginia.*—This paper reports a study of the transient arc 
discharge in air and in hydrogen between a stationary and 
a rotating electrode. A highly damped series circuit con- 
sisting of a condenser, spark gap, and resistance was used 
to form the electrical discharge between the fixed electrode 
and the cylindrical surface of a rotor. It has been previously 
found! that the luminous column of the discharge is carried 
around on the rotor surface. The distance this column is 
carried is as long as 12 cm for an applied voltage of 16 kv. 
This distance has been studied in the present work as a 
function of the rotor speed, the rotor material, and the gas 
pressure. Rotational speeds ranged from 7 to 900 r.p.s. 
(peripheral speeds of 0.17 X 10° to 28 X 10° cm/sec.); rotor 
materials included copper, bronze, duraluminum, and steel; 
gas pressures ranged from one centimeter to atmospheric 
pressure. It was found that the arc column elongates in a 
different manner depending on whether the rotor is anode 
or cathode. Arc current was investigated and found to be a 
function of luminous column length. Some calculations of 
the potential gradient in the arc column and the current 
density at the electrode have been made. 


* Present address Kendal! Mills, Paw Creek, North Carolina. 
asin Beams and L. B. Snoddy, Phys. Rev. 55, 504 (1939); 55, 1145 


2. Simple Relaxation Models to Explain the Dependence 
of Fracture Type of Deformation-Speed and Temperature. 
ARTHUR E. Ruark, University of North Carolina.—There 
are materials (such as tars and asphalts) which break with 
a so-called brittle fracture when suddenly deformed, while 
they exhibit considerable plastic deformation when they 
are slowly loaded. Such effects depend also on the experi- 
mentation-temperature. Since the actual mechanisms of 
flow and fracture may be very complex and may differ con- 
siderably from one material to another, it is advantageous 
to have a simple mechanical model which exhibits the 
essential features of such behavior. We shall discuss a 
suitable model, the essential feature being that in the model 
two or more deformation processes leading to final fracture 


AMERICAN PHYSICAL SOCIETY 





proceed at different rates which depend on the temperature. 
[The temperature dependence is assumed to be on the form 
exp (— E/KT), although the exact nature of this depend- 
ence is not essential. ] 


3. The Biological Action of X-Rays as a Function of 
Intensity. A. A. BLess, University of Florida.—lIt is a 
well-known fact that for the same amount of energy ab- 
sorbed the biological action of alpha- and neutron rays is 
several times as great as the action of x-rays. This differ- 
ence in action was attributed by Zircle and by Failla to be 
due to the higher ionization density produced by alpha- 
particles and by neutron ejected protons compared to the 
density of ionization produced by x-rays. Zircle estimates 
the density of ionization produced by alpha-particles to 
be 100 times as great as that produced by x-rays. The object 
of this work was to investigate whether or not x-rays for a 
given dose administered at high intensity produce effects 
different from those administered at low intensity. Wheat 
seedlings were subjected to given doses of x-rays at in- 
tensities varying from 25 to 5000,r.p.m. No intensity effects 
were observed; the effects are a function of the energy 
absorbed and not of the rate at which the energy is ad- 
ministered. It is possible that the action of alpha- and 
neutron rays is different from that of x-rays in character as 
well as in ionization density. 


4. Dynamics of Rise of Water in a Capillary Tube. Wm. 
A. RENSE, Louisiana State University—Water was caused 
to rise in a tube 0.0632 cm in diameter. The rise was 
photographed in stroboscopic light of interval 1/20 sec. 
The data obtained were found to fit the curve: z=4.65 
+ 1.042 log, t, z being the height and ¢ the time. The nature 
of the rise can be predicted by assuming that the rise is 
turbulent (resistance proportional to the velocity squared). 
If v is the speed: dv/dt= —k*v*, k? being a constant. A 
simple integration yields the form of the equation which 
explained the observed points. An abrupt deviation near 
the end of the rise indicates a change from turbulent to 
streamline flow and suggests a dynamic method of finding 
Reynolds number. 


5. Statistical Variations in the Maximum Load and 
Elongation of Copper Wire and Their Relation to Size 
Effect. RopNEY M. SouirES AND ARTHUR E. RUARK, 
University of North Carolina.—The distribution of maxi- 
mum loads and the distribution of elongations are obtained 
for 200 15-cm samples from a single continuous length of 
No. 24 copper wire. The average deviation of the breaking 
loads is only 0.37 percent and for shorter lengths, properly 
chosen, it is only 0.16 percent. These deviations are much 
smaller than those recorded-for larger and presumably less 
homogeneous bodies of copper. The deviations are not due 
to diameter variations, but to variations in the structure or 
composition of the metal. The smallness of the deviations 
implies a corresponding smallness of size effect. It is shown 
how the test data can be used directly to determine the size 
effect associated with change of specimen length. Larger 
bodies of commercial copper show an average deviation of 
the tensile strength amounting to about 1.3 percent. The 
usual value encountered in tests of glass is 12 to 18 percent. 
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MEETING OF THE OHIO SECTION HELD AT COLUMBUS, OHIO ON May 6, 1944 


HE 1944 spring meeting of the Ohio Section of the American Physical Society was held at 

Ohio State University, Columbus, on May 6. About fifty were present. At the business session 

J. M. Davies, W. E. Singer, and Leon E. Smith were elected as Chairman, Vice Chairman, and 

Secretary-Treasurer, respectively, for the ensuing year; and plans were initiated for a meeting of the 

Section in conjunction with the meeting of the A.A.A.S. scheduled for September 11—16 at Cleveland. 
The programme of the scientific sessions follows. 


Contributed Papers 


The Application of Radioactive Tracers to the Chemistry of Minute Concentrations. J. D. 
KURBATOV AND M. H. Kurpatov, Ohio State University. 

Spontaneous Neutron Emission. M. L. Pooi, Ohio State University. 

Radioactive Isotopes of Tantalum. R. V. ZUMSTEIN, Ohio State University. 

A Comparison of Donders’ and Duane’s Tables of the Amplitude of Ocular Accommodation. 
HENRY WILLIAM HorstettTer, Ohio State University. 

An Animated Diagram of the Oscillating L-C Circuit. J. H. Rusu, Denison University. 

A New Double Monochromator for the Infra-Red Region. R. H. NosBLe, Ohio State University. 

Rotation Vibration Energies of the Non-Linear XYZ Molecular Model. Ropert P. ScHUMAN 
AND W. H. SHAFFER, Ohio State University. 

Infra-Red Spectrum of Methyl Fluoride. K. P. Yates, Ohio State University. 


Invited Paper 
The Use of the Netatron in Nuclear Physics. T. J. WANG, Ohio State University. 


Leon E. SMITH 
Secretary 
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